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FREQUENCY OF SEPARATION AND OF INTERLOCKING OF 
CLUSTERS OF GALAXIES* 


By JERZY NEYMAN AND ELIZABETH L. Scot 
STATISTICAL LABORATORY, UNIVERSITY OF CALIFORNIA, BERKELEY 
Communicated by S. Bochner, June 1, 1953 


1. Introduction.—The present paper is a continuation of the study! 
by the present authors, to be referred to as Paper I, in which a probabilistic 
theory of the spatial distribution of galaxies was developed. Prior to 
explaining the postulates of this theory we mention a probabilistic concept 
that is used throughout the present paper. This is the concept of the 
probability generating function. For any non-negative integer-valued 
random variable » the probability generating function is defined for 
t} S 1 as the expectation of ‘“ and denoted by G,(t). The conditional 
probability generating function subject to a hypothesis // is denoted by 
G,(t| H). 

The postulates of the theory are as follows. 

Postulate (1): Galaxies occur only in clusters. 

Postulate (ii): To every Borel measurable set R in space there corre- 
sponds a random variable y(R) that represents the number of cluster 
centers falling within R. The distribution of y(R) depends only on the 
measure, say V(R) = V of R but not on any other properties of this set. 
If R,, Re, ... isan at most denumerable sequence of disjoint Borel measur- 
able sets, then the corresponding random variables y(R,), y(R2), ... are 
completely independent. 

The probability generating function of y(R) is denoted by G,(t, V). 
It was found that, in order to satisfy (77), it is necessary and sufficient 
that for every set R of finite measure V’ 


VA(t) 


with 


ho +- z. h,t* 
k=l 
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where fy > O and h, 2 O, for k = 1, 2, . and 


yd h,. = hy. 
k=1 


Postulate (111): Let R, be a set of positive measure and R», a measurable 
subset of R;. Let, further, a), do, ..., @m be an arbitrary combination of m 
numbers selected out of 1, 2, ..., m 2 m. Consider the case where R; 
is known to contain exactly m cluster centers, say Ci, Co, ..., C,. It is 
postulated that the conditional probability that R, will contain the m 
cluster centers C,, Cy, ..., Ca,, and no others, is independent of the 
particular combination a), dz, ..., @m but depends only on n, m and, 
naturally, on the measures of R; and Ro». 

Postulate (iv): The number y of galaxies belonging to a given cluster 
is a random variable capable of assuming any positive integer value. 

Postulate (v): The coordinates of any galaxy are random variables. 
Given the coordinates of the center of the cluster, the coordinates of a 
galaxy possess a conditional probability density represented by a con- 
tinuous function f(n) of argument 7, the distance of the galaxy from its 
cluster center. Given the coordinates of all the cluster centers, the 
triplets of coordinates of the galaxies are all completely independent. 

The above five postulates determine the general structure of the spatial 
distribution of galaxies studied. Among other results it was shown 
that no generality is lost in adopting the assumption that y(R) is a Poisson 
variable so that 

G,(t| V) = e~*VO-9 (4) 


where \ represents the expected number of cluster centers per unit volume. 

In order to connect the properties of the model of the distribution of 
galaxies in space with what may actually be observed, it was necessary 
to adopt still another postulate. 

Postulate (vi): Given the distance & of a galaxy from the observer and 
given a fixed observational setup, characterized by the limiting apparent 
magnitude m,, there is a function @(m, &) that represents the probability 
that the galaxy in question will be visible on a properly oriented photo- 
graph. When the distances of several galaxies and the limiting apparent 
magnitudes of the plates are fixed, then the visibility of one galaxy is 
independent of that of any other. 

In a subsequent paper,’® to be described as Paper II, joint with C. D. 
Shane, the above general model was specialized by postulating that the 
universe is static, by assuming (4) and that the unknown function f(n) 
can be approximated by, say 


- l . a ine . 
fini a) = (+) 
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where o > Ois an adjustable parameter. Also, a definite form was adopted 
for #(m, &) resulting from the studies of E. Hubble and M. Humason. 
Finally, at the end of the paper it was tentatively assumed that the dis 
tribution of y — | might be approximated by a negative binomial so that 


G,(t) = (1 + B - Bl) ” (6) 


with a > Oand 8 > 0. 

Using the observational data of Shane and Wirtanen® collected at the 
Lick Observatory, tentative estimates were obtained for all the parameters 
intervening in the specialized model. Unfortunately, these estimates 
depend very much upon the values of the constants determining the 
probability @(m, &) which are only vaguely known. As a result, instead 
of single estimates of the parameters, there is a whole region of plausible 
combinations of values of these estimates. 

In a third paper‘ of the same series, to be referred to as Paper III, by 
E. L. Scott, C. D. Shane, and M. D. Swanson, the consequences of adopting 
the numerical values 


A = 1.6 X 10~' cluster centers per cubic parsec 
a = 5.0 X 10° parsecs 
a= 100, 8 = 915 


were closely examined. It was found that the value of o was probably 
overestimated and that, perhaps, o = 2.5 X 10° is closer to reality. On 
the whole, however, it was found that the implications of the specialized 
model of Paper II are in good agreement with the observations. 

It will be noticed that the original model is capable of representing a 
great variety of distributions. On the one extreme, there is the uniform 
(Poisson) distribution of galaxies. This may be obtained from the model 
by assuming that vy = 1 for all clusters and that y follows a Poisson law. 
On the other extreme, the model implies the possibility that galaxies 
occur in compact clusters isolated from each other by great gaps of empty 
space. Finally, the possibilities covered in the model include distributions 
of intermediate types, that may perhaps be described intuitively as repre- 
senting a continuous field of galaxies with smaller and larger local con- 


centrations. The first of these possibilities is flatly rejected by the 


observations,* ® so that it must be taken for granted that the distribution 
of galaxies in space is not statistically uniform. Also, several large clusters 
are now believed identified and are studied per se. However, information 
appears to be lacking regarding the last two possibilities mentioned, 
namely, whether such clusters as exist are separate isolated entities or 
interlock with each other and should be considered as local condensations 
in the generally uninterrupted field of galaxies. 

Stated in this form, the question is phenomenal rather than mathe- 
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matical in character. The purpose of the present paper is to report on the 
authors’ results that include the translation of the above question into 
mathematical terms, the deduction of certain formulae meant to answer 
the mathematical problem thus obtained, and numerical results based on 
the specialized model of Paper II and the estimates of the parameters 
reproduced above. 

2. Penetration of One Cluster by Galaxies from Other Clusters.-Consider 
a particular cluster I’, of galaxies and denote by C; its center. Also let 
vy, be the (random) number of galaxies belonging to f;. These galaxies 
will be numbered in the order of their distance from C\, beginning with the 
most distant galaxy. Let stand for the distance between C; and the 
kth galaxy. 

In accordance with the model of Paper I, the distance is a random 
variable. The symbol S(x«) will be defined for x 2 0 to represent the 
sphere of radius x centered at C). 

Now consider an arbitrary cluster Il, other than I',, and let v2 be the 
(random) number of galaxies it contains. If a galaxy of the cluster I, 
is located within S(,), then we shall say that this galaxy penetrates I’; to 
the depth k. If every one of the » galaxies of T, penetrates I’; to the depth 
k, then we shall say that the cluster I’; is imbedded in T’, to the depth &. 

Our present probabilistic treatment of the phenomenal problem of the 
interlocking of clusters of galaxies consists in computing the probabilities 
that the cluster I’; will be penetrated by galaxies from other clusters and of 
probabilities that other clusters will be imbedded in [';. The proba- 
bilities of this kind are interesting in relation to various hypotheses re- 
garding v; and yv. For example, we may consider that a cluster of, say, 
five galaxies is not very interesting and may wish to confine our con- 
siderations to clusters TI, containing at least a specified number JN, of 
galaxies. Also, it may be interesting to make similar restrictions regard- 
ing the size of the cluster [,. For example, one may ask how frequently 
a cluster I; containing N; or more galaxies will have imbedded a specified 
number of small clusters [,, each with v2 S Ne» galaxies. Alternatively, 
one may be interested to know how frequently a substantial cluster I’; 
(with », 2 MN, galaxies) is penetrated by galaxies from other substantial 


Accordingly, the results given below are always based on the assumption 
that », 2 N, where N, is a number that may be specified arbitrarily. 
Also, the same results include the possibility of subjecting v2 to a restriction 


either ve S Ne or vw 2 No. The letter // is used to denote any such hy- 


pothesis regarding v.. It appears in the symbols denoting the various 
probability generating functions. Also the symbol P(/7) is used to denote 
the probability, easily computed from G,(¢), that the number v of galaxies 
in a cluster will satisfy the hypothesis //. 
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The proofs of the formulae given below are somewhat long and will be 
published elsewhere. 

3. Penetration of V; by Galaxies from the Nearest Cluster.Given anv 
two clusters of galaxies I’ and ['” with centers at C’ and C”, the term 
distance between I’ and I” will be understood to mean the distance be 
tween the points C’ and C”. With this terminology, let '* be the clustcr 
that is nearest to l';. Let v* be the number of galaxies belonging to I'*. 

The subject of our interest in this section is the random variable X (k) 
defined for | S k S » as the number of galaxies of the cluster I'* thet 
penetrate I, to the depth &. The result given below represents the condi- 
tional probability generating function of X(k), given that », 2 N,; 2 | 
and given that »* satisfies the hypothesis 7. For 1S k S N, we have 


Gyw lt] (na = N)H] = Se tet, k, Ni) So'G,(r| Dduldv (7) 


where 


| 
( 
(k—1)!Ply,= N,} 


P(v, k, Nj) 


. ‘ 
l—v a v*—*Piv, = n} 


(S) 
r=1-— (1 — 2) Flx(u), y(v)] 


Fé,n) = SSS Ale — 22+ y? + 27)'"} dx dy dz 
S(n) 
with x(u) and y(v) defined for 0 S u,v < 1 by 
* |—3 log (1 1 
|/—+ log — 1%) 
x(u) = 
\ tarAP(IT) 


te Sit xif(x) de = 0, a 


The functions defined in (8), (10), and (12) will be used in Sections 4 and 5. 

In order to compute (7) numerically, it is sufficient to know the functions 
fandG,. Obviously, if one substitutes V, = k and interprets // as v* 2 1, 
then the conditional probability generating function (7) becomes the 
unconditional probability generating function of X(k). Formula (7) 
contains the answers to such questions as the following: How frequently 
is a cluster I'\, of at least V, = 100 galaxies, penetrated to the dep‘h 
k = 5 by X(5) galaxies all of which belong to '* which is the nearest of 
all other clusters containing at least Ne = 50 galaxies? Of particular 
interest is the question: How frequently is [; not penetrated by any of the 
galaxies of the cluster ['*? The answer is obtained by substituting k = | 
and ¢ = 0 into formula (7). 
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While formula (7) does provide interesting information regarding the 
interlocking of clusters, this information is a little one-sided. The reason 
is that the given cluster I; may be isolated from the nearest cluster I'* 
(so that X(1) = 0) but yet have penetrations from more distant clusters. 
This problem is dealt with in the next sect on. 

4. The Number of Clusters Which Have Some Galaxies Penetrating a 
Given Cluster.-As formerly, let k be a positive integer not exceeding the 
number »; of galaxies belonging to a cluster ['). The letter v. will be used 
to denote the number of galaxies of a cluster I, other than T;. In this 
section we consider Y(k) defined as the number of clusters I’, that have 
at least one galaxy penetrating I; to the depth k. Thus, if ¥(1) = 0 then 
the cluster I’; is completely isolated. 

The following formula gives the conditional probability generating 
function of Y(k), given that », 2 N,; 2 k and given that v satisfies the 


hypothesis /7. 


Gywltl(n2 N)H] = Stet 8l™ bv, k, Ni) dv (13) 


where 


I(v, HW) = 4eP(H) SZ” #11 — Gl — Fé, y(w))| W)} dé. (14) 


The probability that a cluster containing at least N, galaxies will be 
completely isolated from all those clusters in which the number v of 
galaxies satisfies the hypothesis // is obtained from (13) by substituting 
k = landt = 0. 

5. Number of Clusters of Galaxies That Are Imbedded ina Given Cluster. 
In Sections 3 and 4 we considered the problem of penetration of a cluster 
by galaxies belonging to other clusters. In the present section we are 
concerned with a random variable Z(k) characterizing the frequency with 
which a given cluster I’; includes within its bounds other clusters, with 
all their galaxies. Physically, these clusters imbedded in ['; would prob- 
ably be considered as irregularities in the structure of T). 

The random variable Z(k) is defined as the number of clusters 2, other 
than I, that are imbedded in I’; to the depth k. The conditional proba- 
bility generating function of Z(k), given that the cluster I; contains at 
least »;} 2 N; 2 k galaxies and given that the number v2 of galaxies in the 
cluster I’) satisfies the hypothesis //, is given by the formula 


Gew lt] (m 2 N)H] = J e*t- 9%} & (v, k, Ni) dv (15) 
where 
J(v, W) = dn SS” &G,[ Fé, y(v))| HW) dé. (16) 


6. Some Numerical Results—-As mentioned in the Introduction, 
the available estimates of the parameters \, o, a, and 8 involved in the 
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model of Paper II are still very vague and the detailed numerical study 
of the actual distribution of galaxies using the formulae given above must 
be postponed until these constants are known with reasonable precision. 
Nevertheless, it appears interesting to use such estimates as are now avail- 
able in order to gain some idea as to the degree of mutual interlocking of 
clusters of galaxies. With this in mind, the above general formulae were 
specialized by substituting (5) and (6). Finally, using the constants 
A, a, and @ as given in the Introduction and also two alternative values of 
«= 5 X 10* and o = 2.5 X 10° parsecs, the following results were ob- 
tained. All of them refer to the hypothesis that », 2 1 and »,2 1. 
o = 5X 105 o = 25 X 108 

P\X(1) = 0} 0.0098 0.25 

E({|X(1)] 69.5 27.4 

P\Y(1) = 0} 0.0017 0.20 

E[Y(1)] 14. 1.8 


It will be seen that the conclusion regarding the general appearance of 
the distribution of galaxies in space depends very much on the value of 
o. If the larger of the two values is close to reality, then only about 1 
per cent of all clusters are isolated from the nearest cluster and the 
average number of galaxies of the nearest cluster that penetrate the given 
cluster is 69.5. Also, the percentage of all clusters that are not penetrated 
by galaxies from any other cluster is only 0.17 per cent. Finally, the 


average number of clusters that have galaxies penetrating the given 
cluster is 14. It follows that, with the particular estimates of the param- 
eters, a cluster of galaxies as a physical entity is, probably, a rarity. 
Thus, if one attempts to build a dynamic model of the structure of the 
universe, then this model should consider all the galaxies as forming a 
single complicated system, rather than many separated systems each 
represented by a cluster. 

On the other hand, if o = 5 X 10° is an overestimate (as indicated in 
Paper III) and ¢ = 2.5 X 10° is closer to actual fact, then the above 
figures indicate that a substantial proportion of clusters are isolated from 
the others and, probably, could be considered as separate dynamic systems. 
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THE INFLUENCE OF TEMPFRATURE UPON THE 
INACTIVATION OF A BACTERIAL VIRUS BY X-RAYS 


By C. S. BACHOFER, C. F. Euret, S. MAYER, AND E. L. Powers 


DIvIsIOn OF BIOLOGICAL AND MEDICAL RESEARCH, ARGONNE NATIONAL LABORATORY, 
LEMONT, ILLINOIS; AND DEPARTMENT OF BIoLoGy, UNIVERSITY OF NOTRE DAME, 
NoTRE DAME, INDIANA 


Communicated by W. H. Taliaferro, June 25, 1958 


A discussion of the nature of the mechanisms of radiation damage in 
living systems requires distinction between so-called direct effects on 
sensitive sites or “‘targets,’ and indirect effects, by which the biological 
action is accomplished via some intermediate entity in the system such as 
a free radical produced directly by the radiation. One of the methods 
used to distinguish direct from indirect action consists in measuring 
radiosensitivity of the system while suspended in liquid media of different 
composition. In this manner one presumably can estimate the amount 
of direct action by comparing results obtained in pure water with those 


in some organic medium which, having high affinity for very reactive 
radiation-produced substances, should protect the cells from these inter- 
mediate substances. Such an approach has been used by a number of 
authors, a recent instance being Watson's study of radiation effects in the 
bacteriophages of Escherichia coli." * 


It cannot be presumed, however, that changes in the outside environ- 
ment allow discrimination of direct and indirect effects within the cell 
itself, or even within the virus particle, and for distinction on a subcellular 
and subparticle scale other methods have to be employed. 

One method is a study of the effect of temperature on radiosensitivity. 
Above the very low temperature region temperature changes probably 
should not affect ionization and direct inactivation of a vital target, but 
temperature changes should affect the activity of any intermediate sub- 
stance that affects secondarily the sensitive site, either by changing 
diffusion rates or by affecting reactivity itself, or by both methods. 

In this investigation we have studied the effects of widely separated 
temperatures during x-irradiation on survival of the Tl bacteriophage of 
I. coli. These particles appear to be inactivated by single ionizations 
and they can be dried and carried to very low temperatures with little 
inactivation. In this dried state one is dealing with effects within the 
particle itself, and any separation of effect brought about by the tempera- 
ture differences occurs within the virus. 

A preliminary report has appeared previously.’ 

Materials and Methods.Suspensions in broth of Tl phage were dried 
in aluminum sample cups, irradiated in vacuo at controlled temperatures, 
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and assayed for survival on £. coli, strain B. The method of preparation 
involved distribution of 50 A or 100 X aliquots of broth suspensions of 
approximately 5 X 10° particles per ml. as a drop at the bottom cf each 
cup. The cup was an Aloe-Willett cap of 22 mm. outside diameter, 
modified at the bottom by a depression 14 mm. in diameter and 1 mm. 
in depth. The depression defined the maximum diameter of the drop. 
Such preparations were dried from the liquid state in a 160-mm. vacuum 
top desiccator, by evacuation with a Duo Seal Vacuum pump in a cold- 
room at 0-3°C. Evacuation was controlled to proceed in slow stages, 
with 3-4 hrs. elapsing in the first stage before the pressure was allowed to 
drop below 5 mm. Hg. Inactivation and loss associated with boiling in 
early stages was thus avoided. After the first stage, the pump was allowed 
to run continuously for 16-18 hrs., the pressure falling below 1 mm. Hg 
within the first hr. of this interval. This drying procedure resulted in 
less inactivation and less variability between replicate samples than any 
other method employed, including lyophilization. Nevertheless, the 
data from several experiments had to be rejected on account of excessive 
inactivation attributable to the drying process. In the four experiments 
reported here, survival through dehydration and rehydration was higher 
than 95%. The dried preparations were stored at 25°C., and were kept 
at atmospheric pressure except during irradiation and temperature ex- 
posure. 

The irradiation apparatus included a 50-kv. Picker machine with a 
Machlett beryllium window tube operated at 40 ma. without additional 
filtration. A doughnut-shaped chamber with an internal bore of 20 mm. 
diameter was affixed to the head of the tube. At one end of the bore was 
the Be window and at the other end the sample cup, the irradiated sample 
being 10 ecm. from the target of the tube. Evacuation of the chamber 
was effected through a side-arm connected to the bore and monitored by a 
manometer. Pressure was maintained at less than 0.1 mm. Hg throughout 
the exposure interval. Temperature at the bottom of the cup was moni- 
tored by an iron-constantan thermocouple located on its outer surface, 
while the entire sample cup was immersed in an appropriate isothermic 
medium within a Dewar flask. The isothermic media employed were 


liquid nitrogen (77°K.), solid carbon dioxide-acetone (176-180°K.), 
ice water (273°K.), and water (310°K.). 
The dosimetry of the system employed has not yet been determined 


accurately. Variations in doses delivered were accomplished by varying 
the time of exposures, and all doses are expressed in terms of minutes. 
The best estimate of actual dosage is based on the manufacturer's estimate 
of approximately 62,000 r/min. at 10 cm. As noted below the results 
obtained agree very well with those of others.‘ 

Immediately after irradiation, samples were restored to 25°C. and 
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atmospheric pressure. The phage particles were resuspended in broth, 


and dilution and plating were carried out in conventional fashion.’ The 
maximum time lapse between irradiation and plating was 4 hrs., and 
replicate samples from the beginning or end of this interval indicated ab- 
sence of a standing effect when broth is the suspending medium. Effects 
of temperature and pressure changes per se, if present, were too small to 
be detected beyond the survival limits of the controls. In each experi- 
ment an average of 3 
sample cups per point 





1.0 T T 
0.90 


as employed, and titers 
aii was employed, and titers 


were calculated from 
oH plaque counts of 5 plates 
0.60} per cup. 

Results.—The relation- 
0.50 ship between treatment 
and survival at a _ par- 
0.40 ticular temperature (Fig. 
1) is considered exponen- 
tial and the relation, 


N/Ny = e~*r? (1) 


appears to be valid, with 
dose (D) being in min- 
utes of irradiation under 
the conditions described, 
and ky, being the dose 
inactivation constant in 
reciprocal minutes at the 
absolute temperature. 
i 1 ri jl \ The numerical value of 
2 4 6 8 10 ky (+ its standard error) 
RADIATION, BOBARE, (ain) at each temperature was 
FIGURE 1 estimated by a_ least 
The relationship between survival of bacteriophage Squares approximation 
Tl dried from broth and relative dosage of x-rays at for each exponential line 
indicated temperatures. The solid rectangular point is to be as follows: at 
mae 77°K. itis0.115 + 0.007; 
at 176°K. it is 0.128 + 0.015; at 273°K. it is 0.180 + 0.010; and at 
310°K. it is 0.166 + 0.014. 
These values are plotted against temperature in figure 2. The type of 
response is not determined by the points, but assuming a linear relation- 











ship, we calculate that 
ky = 0.000277 + 0.091. (2) 
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When this is substituted into expression (1), the survival of dried T1 
phage is empirically described as 


N/Ny = e D(0.091 +-0.00027 T') (3) 


when treated and measured under the conditions described herein. 

It is important to note that this expression is known to be valid only 
for temperatures below 37°C. Even at 37°C. there is some thermal 
inactivation with time. For instance, at 35°C., Tl phage in our labora- 
tory dried from broth is inactivated by about 10% over 24 hr.; at higher 
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w 
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FIGURE 2 


The relationship between the dose inactivation constant and tem 
perature for dried coliphage T1 


temperatures more rapid loss is certainly to be expected. Preparations 
held at 25°C. for 24 hrs. and at 35°C. for 4 hrs. were indistinguishable from 
others held at 0-5°C., for the same lengths of time. At the two higher 
temperatures, manipulations, as indicated above, were rapid enough to 
obviate the added complication of simultaneous inactivation from two 
treatments, heat in the higher ranges and x-rays, and no consideration 
is being given here to temperature results above 310°K. 

The general empirical expression relating survival with dose and tem- 
perature indicates that at absolute zero there may be a real dose inactiva- 
tion constant (the intercept of the line in Fig. 2). The latter extension 
has not, of course, been tested. Indeed the authors do not consider 
that 77°K. is close enough to predict with certainty that inactivation 
will proceed at absolute zero at the rate indicated, and there is need for 
tests at temperatures much closer to zero than it has been possible to mak: 
to date. However, the intercept indicated is a measure of the tempera- 
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ture-independent radiation effect, and if it exists it alone is a measure of 
sizes or cross-sections of targets that can be believed to be directly affected 
by radiations. 

Discussion..-Adams and Pollard* report no apparent change in in- 
activation volume at temperatures between 203°K. and 310°K., but it 
should be noted that the spread of the points in their low temperature 
region together with the narrow temperature range tested mask an effect 
of the magnitude reported in this paper. The temperature effects on 
radiation sensitivity reported by them occur above 40°C., and being 
strikingly different in magnitude should probably be considered different 
in mechanism from those reported here. 

In view of the small change with temperature, the possibility was con- 
sidered whether the physical size of the target decreases as the particle 
contracts in the lower temperatures to an extent which accounts for the 
lowered efficiency of the x-rays. We are not aware of any measurements 
of change in volume of phage with temperature, but if we choose as an 
outside condition one of the highest values for coefficient of cubical expan- 
sion of organic materials with temperature cited (tar = 6-8 xX 1074 
volume /unit volume/degree C.), over the range of 77 to 310 we calculate 
a volume increase for this material of only 18.6% (233 K 8 K 1074). It 
is unlikely that the virus particle physically expands with increasing tem- 
perature to this extent. Yet relative ‘target’ size of the phage particle 
at 77°K. as compared with that at 310°K. calculated by comparing e~! 


doses at the two temperatures (from equation 3) is 56°) greater (8.9 
min./5.7 min. = 1.56), nearly three times that which can be accounted 
for in terms of cubical expansion. Therefore, even if cubical expansion 
of the particle contributes to the effect, it seems unlikely that it is more 
than one of the minor reasons why a temperature rise results in increased 


radiation efficiency. 

An analysis of the change with temperature in rate of inactivation 
by x-rays in standard chemical kinetic terms shows that the usual test for 
simple kinetics (logarithm of rate against reciprocal of temperature) using 
the Arrhenius formulation does not yield a straight line, nor, indeed, is 
such a result necessarily to be expected. Furthermore, in the absence of 
much more elaborate evidence, a thermodynamic and kinetic explanation 
critically demonstrating a mechanism or mechanisms is impossible. 

However, a tentative model in kinetic terms is possible. Suppose that 
the inactivation constant is composed of two components, one tempera- 
ture independent and one temperature dependent. The constant observed 
at any given temperature should, then, be separable into these two com- 
ponents by subtraction, and the value which when subtracted leaves 
residues which are in a first order Arrhenius relation (k’ = e~°”/*") with 
temperature is, according to this model, the inactivation constant of the 
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temperature independent portion. It was found that a straight line 
relation is approached as the value of the temperature independent portion 
of the constant approaches the value of the total inactivation constant at 
77°K., and the “activation energy” indicated by the residues obtained 
in this way is about 1100 calories. 

The temperature independent part of the inactivation constant may 
be thought of as the rate at which x-rays inactivate the virus particle by 
affecting directly the sensitive sites. But there are several conceivable 
meanings that could be attached to the temperature dependent portion 
of the inactivation constant. For instance, the A// of about 1100 calories 
or some portion of it can represent the activation energy of some reactive 
compounds produced by radiation. The value estimated here is not out 
of line with the estimated values for free radicals. Or, the estimated 
constant may describe, at least in part, the effects of temperature on some 
diffusion mechanism. That such exist in the solid state is indicated by 
Mayburg and Lawrence’ who describe changes in conductivity of poly- 
ethylene after y-irradiation at different temperatures and who calculate 
Q values as low as 3000 cal., representing to them the control by tempera- 
ture of the migration of protons in this solid system. Other similar 
studies have been made on the relation between temperature and the 
development of latent photographic images.* In biological systems a 
diffusion mechanism has been offered to explain the action of liquid nitro- 
gen temperatures in modifying x-ray effects in plant seeds and _ pollen,® 
with effectiveness of a given dose being reduced by SO“) at the low tem- 
peratures. 

The meaning of this interpretation in regard to exponential survival 
curves is that while the ‘ 
may occur very close to but not necessarily within the sensitive site itself. 


‘single hit’’ interpretation is still allowed, the hit 


The intermediate substances formed are very reactive chemically as 
indicated by the low estimate of over-all activation energy, can migrate 
but short distances before reacting and removing themselves from the 
system, and, therefore, have very short lives. The consequence is an 
appearance at ordinary temperatures of temperature independence and 
direct action. 

The data of course do not require the interpretation described above. 
Indeed, the model suggested by Adams and Pollard‘ which demands the 
simultaneous breakage of several critical bonds can conceivably be made 
to fit. Furthermore, the possibility of competitive reactions with different 


energy requirements, recombinations within vital targets, or several 


combinations of these with the diffusion model are not eliminated. Our 
statement of the energy requirement is an over-all description that does 
not separate one possibility from any other. 

Even without a clear physicochemical explanation, these results impose 
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a refinement on all considerations of target size in T1 coliphage deduced 
from x-irradiation data. Using our inactivation value for 310°K., and 
the rated dosage of the 50 kv.p. x-ray machine, we calculate the sensitive 
volume of T1 in our laboratory to be 2 K 10~'* cm.*, which agrees with 
the value reported by Adams and Pollard‘ for the range 203° to 313°K. 
The results here indicate, however, that the value at very low temperatures 
imay be one-half of that indicated by measurements at room temperature, 
and any consideration of the structure of phage based on x-irradiation 
data is subject to this restriction. 


Acknowledgment.--We wish to acknowledge valuable discussions with 
G. A. Sacher and G. R. Hennig of this laboratory. 
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ELECTRON MICROSCOPY OF THE NUCLEIC ACID 
RELEASED FROM INDIVIDUAL BACTERIOPHAGE 
PARTICLES* 


By DEAN FRASER AND ROBLEY C. WILLIAMS 
Virus LABORATORY, UNIVERSITY OF CALIFORNIA, BERKELEY 
Communicated by W. M. Stanley, May 25, 1953 


It has been known for some years! that certain bacterial viruses can be 
“shocked"’ osmotically by rapid dilution into distilled water of a suspension 
of virus in a concentrated salt solution. As a result of this treatment the 
suspension loses most of its virus activity, and becomes extremely viscous. 
When the mixture is subjected to differential centrifugation, instead of the 
typical virus particles normally observable in the electron microscope! 
only flattened virus ghosts are found in the sediment. Herriott® has shown 
that these ghosts are free of the nucleic acid which makes up some 40% of 
the intact virus, and are composed of protein. The ghosts retain the origi- 
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nal ability of the virus to absorb onto, kill, and lyse the host bacterium 
but do not have the ability to reproduce. By means of studies in which 
the protein portion of the virus was labeled with radioactive sulfur and the 
nucleic acid portion with radioactive phosphorus, Hershey and Chase* 
have shown that the sulfur label does not enter the infected host, but re 
mains attached to the cell wall, whereas the phosphorus label enters the 
cell. Since it has also been found that the sulfur-labeled protein can be 
removed from the infected bacterium after absorption without influencing 
reproduction, the presumption is that only the nucleic acid portion is 
responsible for the reorganization of the cell metabolic processes resulting 
in replication of the infecting virus. It cannot be said categorically that 
the nucleic acid moiety is absolutely protein-free—only that it contains no 
sulfur and no detectable protein by various tests—but the most obvious 
conclusion is that the genetic material of the virus is nucleic acid. From 
all appearances and from the above-mentioned work it also seems apparent 
that the protein portion left behind on infection is essentially the ghost 
seen when the phage is shocked osmotically. It can be concluded that the 
material released by shocking is primarily nucleic acid and that this is the 
genetic substance. In this paper this material will be called nucleic acid. 

Discussion and Results.— Because of the great importance of this genetic 
substance it has seemed to us most attractive to examine the nucleic acid 
released from burst bacteriophages in the electron microscope. Our first 
experiments consisted of examining the supernatant fluid left after osmoti- 
cally shocking and centrifuging a sample of purified bacteriophage, essen- 
tially by the methods previously used, but using in the original, hypertonic 
solution a completely volatile solute such as ammonium acetate or glycerol. 
We found that the preparation contained largely elongated filamentous ma 
terial—what one might expect from the known properties of nucleic acid 
and previous electron micrographs.‘ These experiments, however, had the 
disadvantage that there was no way of identifying any given strands with 
any particular virus particle. 

An obvious improvement would be to shock the virus particles im sifu 
on the electron microscope screen. The techniques previously developed 
by Backus and Williams® for spraying minute droplets and by Williams*® 
for freeze-drying biological specimens for the electron microscope seemed 
particularly well-suited for this purpose. Specimens in water sprayed in 
the usual way® and resprayed, before drying, with a concentrated glycerol 
solution showed a considerable proportion of shocked phage particles where 
sets of droplets overlapped. Stranded material was evident in the shocked 


specimens and completely absent in controls which had received a single 


spraying of either the virus suspension or the glycerol solution. Unfortu- 
nately, however, a glycerol-water solution cannot be sublimed at a practical 
rate at temperatures below —60°C. During sublimation the eutectic 
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concentration is inevitably reached and even at —60°C. the frozen material 
partially melts. 

In the course of exploratory experiments on alternative methods we 
found, somewhat surprisingly, that if a sample of virus sprayed from dis- 
tilled water and dried on the electron microscope screen, either air-dried or 


FIGURE 1 
Electron micrograph of a broken T6 bacteriophage particle. The empty membrane 
of the head, and the intact tail, are shown at the bottom of the picture. The fine 
strands of deoxyribonucleic acid are seen spread out in an approximately circular array 
over most of the micrograph. Shadowed with uranium at tan~*'/;. 94,000. 


frozen-dried, is rewetted by a second spraying of distilled water or by con- 
densed water vapor, the particles take on the flattened, eviscerated ap- 


pearance typical of osmotically shocked phages. In our opinion the phe- 


nomenon is not an osmotic effect, but is induced by the forces of surface 
tension which act on the virus particles on the microscope screen as a drop 
dries or forms. 
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This observation made possible a very simple technique for breaking 
open phages. In the current experiments, the virus particles are sprayed 
from water upon a collodion-covered electron microscope screen, are frozen- 
dried in the usual way,* and are then exposed to moist air, which enters 
the sublimation chamber with the specimen at —20°C. The phage 
particles on the screen serve as nuclei for condensation of the supercooled 
water vapor and are broken open by surface tension. The sublimation 
chamber is then cooled to freeze the condensed water droplets and the frozen 
droplets vacuum-sublimed. The virus concentration is so adjusted that 
-ach condensation droplet contains only a few particles. An isolated phage 
ghost is thus easily identified with the released nuclei acid. 

The results of such an experiment are shown in figure 1. At the lower 
left of the picture the typical phage ghost can be seen and above it a 
collection of the nucleic acid which has unraveled almost completely into 
individual fibrils. The strands are so fine that they are discernible and 
measurable only because of the uranium shadows. Their diameter can be 
approximated by measuring the width of the shadows, and is estimated to be 
20 + 5 A. This value agrees excellently with the known diameter of 
elementary deoxyribose nucleic acid fibrils as determined by x-rays and 
other methods.” 4 The number of strands is most difficult to determine 
because of their very small diameter and their baffling interweaving. At 
the ptesent time no significance can be attached to any of the details of 
structure of the fibrils. 

Summary.—A method has been devised whereby bacterial viruses can 
be broken open im situ on an electron microscope screen in such a way that 
the contents of the phage are released and can be identified with a given 


virus particle. The fibrillar strands which can be presumed to be the nu- 


cleic acid or chromosomal material of the virus, have a diameter of 20 A, a 
value which agrees closely with that obtained by other methods for the 
elementary fibrils of deoxyribose nucleic acid. 


Acknowledgment.—The authors take this opportunity to express their 
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work. 
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TRANSFER OF C' BY LUPINE MITOCHONDRIA THROUGH 
REACTIONS OF THE TRICARBOX YLIC ACID CYCLE* 


By D. O. BRUMMOND AND R. H. Burris 
DEPARTMENT OF BIOCHEMISTRY, UNIVERSITY OF WISCONSIN 
Communicated by C. A. Elvehjem, June 25, 1953 


Many observations of the behavior of organic acids in plants are ex- 
plicable in terms of the functioning of a tricarboxylic acid cycle. Gener- 
ally the most abundant organic acids in plant tissues are malic, citric, and 
oxalic acids. Oxalate is relatively inert and appears to be an end-product 
rather than an active metabolite. As malate and citrate show marked 
quantitative responses to changes in the physical environment of the plant 
and behave otherwise as active metabolites, it has been proposed that 
several of their conversions are mediated through a tricarboxylic acid 
cycle. For example, Pucher, Wakeman, and Vickery! have shown that 
excised tobacco leaves in the dark lose malate and accumulate citrate, 
and Zbinovsky and Burris? have found that tobacco leaves infiltrated 
with C''-labeled malate produce labeled citrate. 

Within the past few years, evidence has accumulated in support of the 
existence of the tricarboxylic acid cycle in plants, and the recent observa- 
tions of Millerd et a/.,* 4 and Bonner and Millerd® compel one to conclude 
that the cycle is operative in certain plant tissues. These workers have 
demonstrated the oxidation of the acids of the tricarboxylic acid cycle and 
the oxidative phosphorylation by mitochondria from mung bean seedlings. 
In the experiments reported here, the presence of the cycle is indicated 
by transfer of radiocarbon from pyruvate-2-C'* into the intermediates of 
the tricarboxylic acid cycle. 

Methods and Materials.—Preparations of particles, identified as mito- 
chondria by morphology and affinity for janus green, were obtained from 
the cotyledons of etiolated white lupine seedlings in a manner similar to 
that described by Millerd ef a/.6 The seeds (Lupinus albus) were soaked 
for 12 hours in distilled water prior to incubation in moist sand at 28°C. 
The cotyledons of seven-day-old seedlings were removed, and 70 g. of 
cotyledons, 70 ml. of 0.2 . sucrose solution and enough solid K,CO; to 
bring the final pH of the brei to 6.5-7.0 were ground together with quartz 
sand in a mortar. All operations preceding the measurement of activity 
were performed at 0°C. and solutions and utensils were prechilled to this 
temperature. The brei was strained through cheese cloth and centrifuged 
at 1000 G for 10 minutes to remove nuclei and cell debris. The super- 
natant was then centrifuged at 20,000 G for 10 minutes, and the yellow 
pellet obtained was resuspended in 0.2 A/ sucrose and recentrifuged at 
20,000 G for 15 minutes. This yellow sediment resuspended in 20-25 
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mil. of 0.2 M sucrose was used as the enzyme preparation. The suspension 
contained 0.5-1.0 mg. N/ml. and constituted half the volume of the 
reactants in the manometric flask. All of the experiments were performed 
at 30°C. in conventional Warburg manometric apparatus. 

Pyruvate, oxalacetate, a-ketoglutarate, and oxalosuccinate were sepa- 
rated as their 2,4-dinitrophenylhydrazones on silica gel columns by a 
method similar to that employed by Sanger‘ for the separation of the 2, 4- 
dinitro-N-phenyl derivatives of amino acids. The columns were prepared 
by uniformly mixing 1.4 ml. of an aqueous buffer solution, which was 
equilibrated against the organic solvent used for developing the column, 
with 2.0 g. of silica gel (Mallinckrodt, 100 mesh, analytical reagent silicic 
acid). About 25 ml. of the organic phase was mixed with the silica gel, 
and the resulting slurry was poured into a | cm. diameter column. The 
solvent systems used were 5 per cent butanol, 95 per cent chloroform (ail 


solvents given as volume per cent) equilibrated against a 1.0 AJ potassium 
phosphate buffer of pH 6.8; 15 per cent butanol, 85 per cent chloroform 
equilibrated against a 1.0 .\V potassium formate buffer of pH 3.8; and 7 
per cent butanol, 93 per cent chloroform equilibrated against a 1.0 M 


potassium phosphate buffer of pH 2.0. 

The hydrazones of the keto acids were prepared by mixing an excess 
of a 0.1 per cent solution of 2,4-dinitrophenylhydrazine in 2.0 N HCl 
with the free acids or with the preparation to be determined according to 
the method of LePage.’ The hydrazones then were extracted with ether 
in a Kutscher-Steudel extractor for | hour; the ether was removed under 
vacuum, and the residue was taken up in 0.5-1.0 ml. of the organic de- 
veloping solvent and was placed on the appropriate column. For separa- 
tion of the oxalosuccinate derivative, the 7 per cent butanol, 93 per cent 
chloroform solvent was employed. Derivatives of oxalacetate and a- 
ketoglutarate were separated with the 15 per cent butanol, 85 per cent 
chloroform solvent, and the pyruvate 2,4-dinitrophenylhydrazone was 
recovered with 95 per cent chloroform, 5 per cent butanol as solvent. 
In these systems, the pyruvate and oxalacetate 2,4-dinitrophenylhydra- 
zones were resolved into their geometric isomers. Recoveries of 98-101 
per cent were obtained from known pyruvate, oxalacetate, and a-keto- 
glutarate when about 100 ug. quantities of these acids were determined. 
The quantity of each 2,4-dinitrophenylhydrazone was measured in a 1 cm. 
cell in the solvent in which it was eluted from the column; the quantitative 
observations were made in a Beckman model DU spectrophotometer at 
the wave length of maximum absorption for each derivative (range 340 
to 365 my). 

Experimental Results and Discussion.-As shown in figure |, the addition 
of magnesium ions and ATP (adenosine triphosphate) stimulated the 
oxidation of pyruvate plus malate by mitochondria, but Armour Labora- 
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tories’ “‘liver coenzyme concentrate’ stimulated much more. Figure 2 
indicates that the addition of DPN of greater than 95 per cent purity 
stimulated the oxidation of malate plus pyruvate. Armour Laboratories’ 
which contains DPN, TPN (di- and tri- 
phosphopyridine nucleotides) and CoA (coenzyme A) and which stimu- 
lated the uptake of oxygen, was added routinely in later experiments. 


“liver coenzyme concentrate, 


100) wl. OXYGEN 
ul. OXYGEN UPTAKE 
UPTAKE 











0 ¥ v 7 r T ] 
oO 40 80 120 fe) 20 8640 60 


MINUTES MINUTES 


FIGURE 1 FIGURE 2 
Fig. 1—-The influence of ATP (adenosine triphosphate) plus MgSO,, and of CoA 
(coenzyme A) on the oxidation of pyruvate plus malate by lupine mitochondria. Each 
flask contained mitochondria, 50 wm potassium phosphate buffer pH 7.0, 20 um L 
malate, and 20 wm pyruvate; final volume 2.0 ml. Temperature, 30°C. (17) 1 mg 
“liver coenzyme concentrate’ added. (2) 2 um MgSO, plus 0.2 um ATP. (3) Nothing 
added. Fig. 2-—-The influence of diphosphopyridine nucleotide and coenzyme A on the 
oxidation of pyruvate plus malate by lupine mitochondria. Substrates and buffers 
supplied as in figure 1. (/) 2.5 mg. “liver coenzyme concentrate” added. (2) 
0.25 or 0.05 mg. DPN added. (3) 0.25 mg. “liver coenzyme concentrate” added. 





In contrast to the preparations of Millerd ef a/.,* it was unnecessary to 
grind the lupine tissue in phosphate buffer; added cytochrome-e did not 
influence the uptake of oxygen. The rates of oxidation of pyruvate in the 
presence of L-malate or oxalacetate are shown in figure 3. As oxalacetate 
inhibited oxygen uptake for the first 30 minutes, L-malate was used subse- 


quently as the ‘“‘sparking”’ acid. 





a 
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Manometric measurements indicated that malate and pyruvate in 
combination were oxidized at a more rapid rate by the mitochondrial 
preparations than the sum of their rates singly. These preparations thus 
seemed to provide an enzymatic system suitable for testing incorporation 
of a tracer into intermediates of the tricarboxylic acid cycle. Tests were 
made by adding to a mitochondrial preparation 200 um of pyruvate-2- 
Cl, 200 um of L-malate and 200 um of another acid of the cycle to serve 
as a trapping acid. During an hour’s incubation at 30°C. in a 125-ml. 
Warburg flask, the trapping acid 
equilibrated with the same acid as it 300) pl. OXYGEN 
was formed enzymatically. Then the UP TAKE 
trapping acid and any other acids 
present in sufficient quantity were 
isolated by partition chromotography 
on silica gel’ and were analyzed for 
their content of C'. 

Data presented in table 1 show that 
there is an appreciable incorporation 
of C' from pyruvate-2-C' into every 
acid of the tricarboxylic acid cycle; 





this indicates the participation of 





each acid in the cyclic oxidation of * ¥ 
‘ 80 160 


pyruvate. In the experiment in which 

malate was supplied, there was a net MINUTES 

: : ee FIGURE 3 

accumulation of 21 um of citrate. = 

TI fi ae f tl * rhe effect of malate and oxalacetate 
1e specific activity of the citrate was on the oxidation of pyruvate by lupine 

32.6 counts X min.~' X wm! and mitochrondria. All flasks contained 


that of the original pyruvate 31.5 50 um potassium phosphate buffer pH 
counts X min.-' X wm-'. Thus it 79 | ms. “liver coenzyme concen 


is evident that the citrate arose 


trate’ and mitochondria in a final 
volume of 2.0 ml. (J) 10 wm malate 
directly from pyruvate and malate, plus 20 ym pyruvate. (2) 10 ym 
since the molar specific activities of malate. (3) 10 um oxalacetate plus 20 
citrate and pyruvate were essentially = 4" AOS (4) 10 um oxalacetate or 
the same. The substantial accumu- ~ 4” P¥™VAate —— 
lation of citrate suggests the presence of an active condensing enzyme which 
can effect a condensation of acetyl CoA and oxalacetate as occurs in 
isolated animal tissues.* That the subsequent oxidation of pyruvate 
proceeds via the other acids of the tricarboxylic acid cycle is shown by 
the accumulation of appreciable amounts of C'‘ into every acid of the cycle. 
The variability of the amount of C'* found in individual acids after each 
experiment probably resulted from variability in recovery as well as 
differences in the activity of specific preparations. The accumulation 
of C'* in the respired CO, shows that the oxidation of a portion of the 
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pyruvate is complete. The small percentage accumulation of C' in the 
respired CO, in these experiments can be attributed to the effective trap- 
ping of the C'* by the pools of acids added as sparker and trapping agents. 
This small accumulation of C'O, also constitutes evidence for the asym- 
metrical synthesis and oxidation of citric acid in plant tissues, a process 
previously demonstrated in animal tissue.'? As much C!* should have 
appeared in the CO, as appeared in succinate, fumarate, malate, and 
oxalacetate if the conversion had been symmetrical. 

These results, in conjunction with earlier evidence that mitochondria 
can readily oxidize the intermediates of the tricarboxylic acid cycle, 
strongly support the premise that pyruvate oxidation by lupine mito- 
condria proceeds via the conventional tricarboxylic acid cycle. As plants 
can readily form pyruvate from carbohydrates, it is evident that carbo- 
hydrates in plants can be completely oxidized to CO, and H,O in a manner 
analogous to that of the animal cell. How widespread and active the tri- 


TABLE 1 
THE ACCUMULATION OF C!! IN COs AND ACIDS OF THE TRICARBOXYLIC AcID CYCLE 
FROM PyRuVATE-2-C!! SupPLIED TO PREPARATIONS OF LUPINE MITOCHONDRIA 
PER CENT C!4 PER CENT C!l4 PER CENT C!4 
FOUND IN PER CENT C!4 FOUND IN FOUND IN 


TRAPPING ACID FOUND IN CITRATE MALATE 
TRAPPING ACID ADDED RECOVERED COs FORMED FORMED RECOVERED 


Citrate 6.5 l ; 0.7 
Cis-aconitate l 0.5 30.2 0.6 
Isocitrate l 49.0 5.0 
Oxalosuccinate (30 min. 
reaction ) 
a-Ketoglutarate 
Succinate 
Fumarate 
Malate 
Oxalacetate 
Pyruvate alone 


carboxylic acid cycle is in plants and whether it is functional in mature 


plant tissues remain to be established. 

Summary.-— Methods for the chromatographic separation and determi- 
nation of the 2,4-dinitrophenythydrazones of pyruvate, a-ketoglutarate, 
oxalacetate, and oxalosuccinate are described. Recoveries greater than 
98 per cent for 100 wg. quantities of pyruvate, a-ketoglutarate, and oxal- 
acetate are obtained; the range for each acid is 10-500 ug. 

Mitochondria isolated from the cotyledons of etiolated lupine seedlings 
oxidize pyruvate via the tricarboxylic acid cycle in a manner analogous 
to that of animal mitochondria. Individual experiments in which py- 
ruvate-2-C'! plus unlabeled malate were oxidized by lupine mitochondria 
in the presence of one other acid of the tricarboxylic acid cycle demon- 
strated incorporation of C'* into each acid of the cycle so tested. 
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SUBCELLULAR FRACTIONATION OF MOUSE SPLEEN RADIA 
TION PROTECTION ACTIVIT Y* 
By LEONARD J. CoLe, Maurice C. FIsHLER, AND Victor P. Bonpt 


BIOLOGICAL AND MEDICAL SCIENCES Division, U. S. NAVAL RADIOLOGICAL DEFENSE 
LABORATORY, SAN FRANCISCO, CALIFORNIA 


Communicated by Curt Stern, June 9, 1953 


Jacobson, ef a/.' have demonstrated that the survival of mice exposed 
to lethal doses of x-rays can be increased significantly by implantation of 
normal mouse spleens into the peritoneal cavity of irradiated mice after 
exposure. These findings have led to the concept of a humoral radiation 
recovery factor (or factors) present in, and elaborated by, normal mouse 
spleen tissue.” They were followed by experimental studies from this Lab- 
oratory,* * in which it was shown that the intraperitoneal injection of 
neutral phosphate buffer homogenates of normal mouse spleen into mice 
which had received otherwise lethal whole body x-irradiation, afforded 
marked protection against mortality. These studies provided a reproduc- 
ible experimental procedure for obtaining active, injectable spleen prepa- 
rations in which the vast majority of the cells had been disrupted, and in 
which the level of organization was far removed from that of the original 


spleen tissue. Furthermore, the experimental procedure provided a suit 
able biological assay for testing the protective activity of various spleen 
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fractions and for investigating the nature of the spleen factor. The su- 
pernatant obtained after centrifugation of the spleen homogenate at 1400 
G for 16 minutes elicited no protective effect, while protective activity was 
retained in the sediment. In view of these results, it was considered plaus- 
ible that the spleen factor—if a chemical entity—might be associated 
with, or bound to certain particulate or structural components comprising 
the spleen cells. In the present studies, normal mouse spleen tissue was 
subjected to subcellular fractionation by the Schneider and Hogeboom 
method of differential centrifugation in sucrose medium,® ® and the re- 
spective subcell fractions assayed for radiation protection activity. The 
experimental data indicate that protective activity or mouse spleen homog- 
enates is associated with the spleen nuclei fraction, under the conditions 
of these experiments. 

Material and Methods.The radiation procedure was essentially iden 
tical with that described previously. A Westinghouse Therapy Unit 
was the x-ray source. Radiation factors were: 250 kvp.; 15 ma.; 
0.5 mm. Cu plus | mm. Al filter; skin-to-target distance, 100 cm.; 
dose rate approximately 25 r per minute, measured in air with a Victoreen 
r-meter. In each experiment the control and experimental animals were 
matched with respect to sex, age, and body weight, irradiated simulta- 
neously, and caged together, eight or ten per cage. The irradiated animals 
were LAf, mice of both sexes, 10 to 14 weeks old, weighing 20 to 25 g.; 
the spleen donors were non-irradiated young LAf, mice (approximately 
| week old). 

Young (1-week-old) non-irradiated LAf, mice were sacrificed by cervical 
dislocation. The spleens were removed immediately (using clean but not 
aseptic technique), transferred to a tared chilled, sterilized Ten-Broeck 
glass homogenizer, weighed, and homogenized in 3 ml. of chilled 0.25 M/ 
sucrose solution. During this process the homogenizer was kept immersed 
in a water-bath containing cracked ice. The homogenate was diluted with 
0.25 M sucrose solution to a final concentration of approximately 100 mg. 
spleen per ml. The diluted homogenate was then transferred to a chilled, 
graduated 15-ml. centrifuge tube, and centrifuged in a refrigerated Inter- 
national No. 2 centrifuge (5°C.) at 600 G for 10 minutes. The superna- 
tant was withdrawn carefully with a pipette and transferred to a flask. 
The sediment was transferred to a clean homogenizer, rehomogenized in 
fresh 0.25 M sucrose solution, and made up to the original volume with 0.25 
M sucrose solution. The homogenate was again centrifuged at 600 G for 
10 minutes, the supernatant removed as before, and combined with that 
from the first centrifugation. The sediment from the second centrifuga- 
tion, resuspended in 0.25 M sucrose solution, constitutes the spleen nuclei 
fraction. This fraction was turbid, whitish, and was easily injected 
through a No. 26 needle. The combined supernatants from the above 
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were then subjected to centrifugation in a refrigerated Model E Spinco 
preparatory ultracentrifuge at 8500 G for 10 minutes. The rotor was pre- 
cooled to 5°C., and the vacuum chamber was refrigerated throughout the 
run. The supernatant fluid was carefully removed with a medicine drop- 
per, and the sediment, constituting the mitochondria fraction, was resus- 
pended in fresh sucrose solution. The supernatant from the mitochondria 
fraction was then centrifuged in the Spinco ultracentrifuge at 20,000 G 
for 90 minutes, yielding a sediment, the microsome fraction, which was re- 
suspended in fresh 0.25 Af sucrose. The supernatant from the microsome 
fraction was a pink, crystal-clear, optically homogeneous fluid. It con- 
stituted the soluble fraction. The individual subcellular fractions were 
administered without additional washing. 

In the earlier experiments the medium used for subcellular fractionation 


was an aqueous solution of 0.25 \/ sucrose made up with distilled water. 
Since it was found that the protective activity of spleen homogenates pre- 
pared in this medium was not of a high order, it was decided to use a modi- 


fied, sucrose medium containing added salts. The sucrose-salt solution 
had the following composition: 0.242 / sucrose, 0.0094 1f KH.PO,, 
0.0125 M K,HPO,, 0.0015 7 NaHCOs, and 0.0006 7 adenosine triphos- 
phate (Na ATP). This solution represents a modification of one used by 
Wilbur and Anderson’ for isolating cell nuclei. It is isosmotic, fortified with 
potassium salts, bicarbonate, and adenosine triphosphate, and of suitable 
density to permit subcellular fractionation. 

Thirty minutes after harvesting, the subcell fractions were injected either 
intraperitoneally or intravenously into the irradiated mice, usually be- 
tween 2 and 3 hours after radiation exposure. Aliquots of each fraction 
injected never exceeded 1.0 ml., and the amount injected per mouse is ex- 
pressed in terms of equivalent weight of spleen tissue from which it was 
derived. The control, irradiated mice, received injections of equivalent 
volumes of 0.25 M sucrose solution, or sucrose-salt solution. Survival up 
to 30 days post-irradiation and body weight changes were used as criteria 
for evaluation of protective activity in the spleen subcell fractions. 

Both wet and stained dry preparations of whole spleen homogenate and 
of the particulate fractions were examined under the microscope. Cell 
counts were made using a standard hemacytometer. 

Results.—It is apparent from examination of the data in table | that the 
0.25 M sucrose solution was unsatisfactory as a medium for the prepara- 
tion of active spleen homogenates. Whole spleen homogenates prepared 
in this medium exhibited some radiation protection activity, but to a far 
lesser degree than that reported previously for whole spleen homogenates 
in phosphate buffer at pH 7.2. Of the subcell fractions isolated in this 
medium, there was a suggestion of protective activity in the nuclei frac- 
tion, and to some extent also in the mitochondria. Whole spleen homog- 
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enates prepared in 0.88 A/ sucrose were completely ineffective in radiation 


protection. 

In one experiment, spleen homogenate was prepared in a 0.25 A/ sucrose 
medium to which adenosine triphosphate (Na ATP, 5.7 & 1074 A/) had 
been added. Of ten mice, exposed to 750 r, each of which received a single 
intraperitoneal injection of this homogenate (equivalent to 56 mg. of spleen), 


TABLE 1 


Turrty-Day Morta.ity OF [IRRADIATED MICE INJECTED WITH SUBCELLULAR SPLEEN 
FRACTIONS ISOLATED IN 0.25 M Sucrose MEDIUM 
EQUIVALENT 
RADIA AMOUNT OF 
TION SPLEEN 
DOSE INJECTED, SURVIVAL AT 30 DAYS 
ee FRACTION INJECTED MG NO./ TOTAL PER CENT 
700 Spleen homogenate in 0.25 A/ sucrose 48 6/10 60 
700 Nuclei fraction é 0/10 0 
700 Supernatant from nuclei 0/10 0 
700 0.25 M sucrose 0/10 0 


750 Spleen homogenate in 0.25 M sucrose /15 
750 Nuclei fraction 

750 Mitochondria fraction 

750 Microsome fraction 

750 Soluble supernatant fraction 


750 0.25 M sucrose 


750 Spleen homogenate in 0.25 M sucrose 
750 Spleen homogenate in 0.88 M sucrose 
750 0.25 M sucrose 
750 0.88 M sucrose 
750 Spleen homogenate in 0.25 M sucrose 
750 Supernatant—70 G for 10 min. 42 
750 Residue from above 42 
750 0.25 M sucrose 
750 Spleen homogenate in 0.25 M sucrose’ 56 
750 Nuclei fraction 280 
750 Mitochondria fraction (intravenous ) 250 
750 Mitochondria fraction 220 
750 0.25 M sucrose + ATP 
* Donors, 2 weeks old 
» Contains 6 K 10~4 M adenosine triphosphate (ATP) 
All injections were intraperitoneal, except as otherwise indicated. 


eight were alive 30 days after radiation exposure, whereas none of the ir- 
radiated control mice receiving sucrose plus ATP survived. The spleen 
nuclei fraction, isolated in this medium, exhibited some protective activ- 
ity: three out of nine survived of the irradiated mice each injected with 
the nuclei fraction (equivalent to 280 mg. spleen) the average survival 
time of the non-survivors was 15.5 days, compared to 9.2 days for the su- 
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crose controls. The mitochondria fraction when injected intraperitoneally 
(equivalent to 220 mg. spleen) also afforded a small degree of protection: 
two survivors out of seven; average survival time of non-survivors was 
9.4 days. Intravenous injection of mitochondria did not afford protection 
against death; average survival time was | 1.2 days. 


TABLE 2 
THIRTY-DAY MORTALITY OF IRRADIATED MICE INJECTED WITH SPLEEN SUBCELLULAR 
FRACTIONS ISOLATED IN SUCROSE-SALT MEDIUM 
EQUIVALENT 
AMOUNT OF 
TION SPLEEN 


DOSE, INJECTED, SURVIVAL AT 30 DAYS 
R FRACTION INJECTED MG NO./TOTAL PER CENT 


RADIA 


750 Spleen homogenate in sucrose-salt 56 7/17 41 
750 Nuclei fraction 280 9/9 100 
750 Mitochondria fraction 250 0/8 0 
750 Microsomes + soluble supernatant 140 0/8 0 
750 Sucrose-salt solution rue 0/10 0 
750 Spleen homogenate in sucrose-salt 56 10/10 100 
750 Nuclei fraction 280) 10/10 
750 Mitochondria (intravenous ) 250 0/8 

Nuclei fraction 196 

Nuclei fraction 98 

Sucrose-salt solution 

Spleen homogenate in sucrose-salt 39 

Nuclei fraction 

Nuclei fraction in 1 M NaCl 

Controls (no injection ) 

Spleen homogenate in sucrose-salt 

Nuclei fraction 

Salt-sucrose solution 

Spleen homogenate in sucrose-salt 

Nuclei fraction in 0.5 M NaCl 

Salt-sucrose solution 

0.5 M NaCl solution 
750 Nuclei fraction in 1 M NaCl, diluted 

to 0.75 M NaCl 250) 

750 0.75 M NaCl solution 


All injections were intraperitoneal, except as otherwise indicated 


Subcell Fractionation in Sucrose-Salt-ATP Medium.—The results in- 
dicate that the modified sucrose-salt-ATP solution provided an excellent 
medium for the preparation of protective spleen homogenates, which could 
then be subjected to subcellular fractionation. The survival data are 
summarized in table 2. The degree of post-irradiation protection afforded 
by the injection of spleen homogenates prepared in the sucrose-salt medium 
was greater than that observed with comparable homogenates prepared in 
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0.25 M sucrose medium. This is also illustrated by comparative body 
weight changes of mice exposed to 75( r and treated with spleen homoge- 
nates prepared in both media. (See figure 1.) 

In a series of five separate experiments, marked protection against radi- 
ation mortality was afforded by a single injection of the spleen nuclei 
fraction. In two experiments 100 per cent of the irradiated mice (750 r), 
each injected with nuclei fraction from ten spleens, survived the 30-day 
observation period, whereas all of the irradiated controls (750 r) injected 
with the sucrose-salt solution died. Apparently, the amount of spleen 
substance in the form of the nuclei fraction required for protective activity 
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FIGURE 1 


Body weight changes in x-irradiated mice treated with spleen homogenates prepared 
in sucrose-salt medium and in 0.25 M sucrose 


is greater than that in the form of whole spleen homogenate. However, 
some activity was observed with the nuclei fraction derived from as little 
as 3.5 spleens (98 mg.) was administered to irradiated mice. 

No evidence of radiation protection activity was observed in the mito- 
chondria fraction obtained by differential centrifugation in the sucrose- 
salt medium (table 2), in contrast with a small increase in the number of 
survivors in the group of irradiated mice injected with the mitochondria 
fraction isolated in 0.25 M sucrose medium. Since it is known that mito- 
chondria become aggregated and agglutinated in the presence of salts, 
and as a result are sedimented with the nuclei,’ it seems likely that the nu- 
clei fraction isolated in the sucrose-salt medium contains a considerable 
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admixture of aggregated mitochondria. These considerations suggest 
that a small degree of radiation protection activity may well be associated 
with the mitochondria fraction, but that the apparent lack of activity in 
the mitochondria isolated in sucrose-salt medium may be ascribable to the 
smaller amount of mitochondria injected. 

Desoxyribonucleic acids, in association or combination with protein, 
comprise the major portion of the cell nucleus. Furthermore, Schneider" 
has shown that all of the desoxyribonucleic acid in liver, kidney, and tu- 
mors is recovered in the nuclei fraction. Desoxyribonucleic acid consti- 
tutes as much as 48.5 per cent by weight of the dry, fat-free salmon sperm 
nuclei!! and indeed, among the various tissues of the animal body, the spleen 
is one of the richest in desoxyribonucleic acid content.'*? On the basis of 
the present finding that the radiation protection activity of mouse spleen 
homogenates is associated with the cell nuclei fraction, an attempt was 
made to determine whether crude nucleic acid extracts from the nuclei 
fraction would exhibit radiation protection activity. The procedure of 
Mirsky and Pollister'* involving the extraction of nucleoprotein by 1 A 
NaCl was employed. The spleen nuclei fraction was first prepared in 
sucrose-salt medium as described above; one-half of the nuclei suspen- 
sion was injected into a group of eight mice approximately | hour after 
exposure to 750 r. Solid NaCl, to a final concentration of 1 1/ NaCl, was 
added to the remaining half of the spleen nuclei suspension in a flask im- 
mersed in an ice-water bath. Immediately upon addition of the salt, the 
nuclei suspension became very viscous and quite resistant to flow. Since 
it could not be readily subjected to centrifugation, the viscous mass was 
administered directly to a second group of irradiated mice (750 r). The 
material was administered by intraperitoneal injection through a No. 23 
needle, after packing it into the barrel of the hypodermic syringe with a 
glass rod. The data in table 2 show that the viscous | 17 NaCl suspen- 
sion of the nuclei fraction was devoid of radiation protection activity, as 
was the spleen nuclei fraction suspended in 0.5 17 NaCl solution, 

Inactivation of Spleen Homogenates by in Vitro Exposure to X-Rays and 
to Ultrasonic Vibrations. Additional experimental data pertinent to the 
present discussion show the effect of x-radiation, and of ultrasonic vibra- 
tions on the radiation protection activity of mouse spleen homogenates. 
The homogenate was prepared and divided into two equal portions. One 
portion was exposed to the physical agent while the other portion, although 
kept under otherwise identical conditions, was not exposed. The treated 
and untreated homogenates were then respectively injected into two groups 
of irradiated animals. The source of the ultrasonic radiations was a 


crystal Ultra-Sonorator (Crystal Research Laboratories). The spleen 


homogenate, contained in a Pyrex test-tube, was immersed in the oil bath 
of the Ultra-Sonorator during the exposure. In one experiment, the 
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homogenate was exposed continuously for 5: minutes at 1010 v. and 160 
ma. In the second experiment, the ice-chilled homogenate was exposed 
at 750 v. and 100 ma. for a total time of 4 minutes; the exposure was al- 
ternated | minute in the Ultra-Sonorator and 1 minute in an ice-water 
bath. It seems clear from the data presented in table 3 that both x-ray 
exposure (750 r), and ultrasonic vibrations, in vitro, inactivate the protec- 
tive activity of mouse spleen homogenates. Relative body weight changes 
in x-irradiated mice treated with normal spleen homogenate, and with 
spleen homogenate exposed to 750 r are shown in figure 2. 

Microscopic Observations.—An attempt was made to ascertain the num- 
ber of whole cells, broken cells, and cell nuclei in whole spleen homogenates 
and in nuclei fraction. The attainment of this objective is far more dif- 


TABLE 3 
INACTIVATION OF SPLEEN HOMOGENATE RADIATION PROTECTION ACTIVITY BY IN VITRO 
EXPOSURE TO X-RAYS, OR TO ULTRASONIC VIBRATIONS 


RADIATION 
DOSE SURVIVAL AT 30 DAYS 


R MATERIAL INJECTED NO./TOTAL PER CENT 
650 Spleen homogenate (32 mg. )* 10/10 100 
650 Spleen homogenate (32 mg.) exposed to 

ultrasonic vibrations 8/17 47 
650 Phosphate buffer 8/13 
700 Spleen homogenate (36 mg.) 14/15 
700 Spleen homogenate (36 mg.) exposed to 

ultrasonic vibrations 3/15 
700 Phosphate buffer 1/10 
750 Spleen homogenate (56 mg.) 7/10 
750 Spleen homogenate (56 mg.) exposed to 

THO r 0/10 


750 Phosphate buffer 0/10 


@ Equivalent amount of spleen. 


ficult than would appear at the outset, since adequate microscopic criteria 
for distinguishing intact spleen cells from partially disrupted ones and from 
nuclei, are not readily established. In order to obtain base line criteria 
for the normal appearance and size of intact spleen cells, imprints of young 
mouse spleens were prepared. The freshly cut surface of the spleen was 
touched lightly on a glass slide; the slide was then air-dried and stained 
with Wright-Giemsa stain. Lymphocytes of all sizes were seen, as well 
as reticulum cells, megakaryocytes, normoblasts, polymorphonuclear 
leukocytes, blast cells, and an occasional eosinophil. The vast majority 
of the cells ranged in diameter from 6 » to 13 yu. A few macrophages 
19-20 uw in diameter were seen. The cell nuclei stained dark purple and 
comprised the largest proportion of the cell diameter. The cytoplasm ap- 


peared as pale blue rims around the nuclei. 
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Microscopic examination of the stained preparations of phosphate buffer 
spleen homogenate revealed that the vast majority of cells in this prepara- 
tion are obviously disintegrated or disrupted. A few apparently intact 
cells resembling small lymphocytes were seen; the remainder of the prepa- 
ration consisted of a great number of deeply staining, free nuclei, a num- 
ber of partially disintegrated cells in which portions of the cell surface 
had been torn away, and granular cell debris. Cell counts indicated the 
presence of approximately 4.3 X 10° cellular bodies, i1.e., whole cells plus 
disrupted cells plus nuclei, in sucrose-salt homogenates derived from 100 
mg. of spleen. Of these bodies, approximately 3.6 X 10° appeared to be 
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FIGURE 2 


Effect of in vitro exposure of spleen homogenates to x-rays (750 r) on radiation pro 
tection activity 


intact cells. On the basis of the value of 22.9 & 10° cells per milligram of 
normal mouse spleen recently reported by Mizen and Petermann,'* it was 
calculated that the single homogenization involved in the preparation of 
whole spleen homogenates had disrupted approximately SO per cent of 
the spleen cells beyond recognition as cells or cell bodies. Of the remaining 
20 per cent recognizable as cell bodies, approximately one-fifth appeared 
to be whole cells. Thus, the spleen homogenate contains about 2 per cent 
of the original number of cells in seemingly whole form. It will be re- 
called that the whole homogenate derived from one mouse spleen will 
afford marked protection against radiation mortality. Assuming that the 
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remaining whole cells after homegenization are representative of the spleen 
population, implantation of one-fiftieth of one spleen would almost cer- 
tainly not provide the same order of radiation protection (if any) as that 
afforded by spleen homogenate derived from one spleen and containing 
the same number of whole cells as that contained in one-fiftieth of one 
spleen. The additional protection, then, must have come from the dis- 


rupted cells. 

Microscopically, the nuclei fraction consisted largely of nuclei which 
were, in many instances, disrupted. The number of seemingly whole cells 
in this fraction was very low, and indeed, many microscopic fields were en- 
tirely devoid of whole cells. Visual comparison of the smear of the whole 
homogenate and of the nuclei fraction indicates that the procedure for 
nuclei isolation, involving, as it does, two separate homogenization steps in 
sucrose-salt medium, results in a great reduction of the number of appar- 
ently intact cells as compared with the original homogenate. Smears were 
also prepared from normal unhomogenized spleen fragments, and treated 
identically with the smears of spleen homogenates and the nuclei fraction. 
Practically all the cells in such smears appeared to be intact, indicating 
that the method of smear preparation per se did not disrupt the normal ap- 
pearance of the cells. 

Discussion.— The present finding of the association of radiation protec- 
tion activity of mouse spleen homogenates with the cell nuclei fraction 
bring to the fore the problem of the radiosensitivity of the cell nucleus, 
and the role of nucleic acids and nucleoproteins as determinants in the 
biological effects of x-rays. Certainly, the earliest damaging effects of 
ionizing radiations, such as inhibition of cell division, are most readily 
observed in the nucleus. The well-known induction of chromosome breaks 
and mutations by ionizing radiations provide additional examples of the 
great vulnerability of the cell nucleus to radiation damage. An excellent 
review of the problem of radiation sensitivity of the nucleus has been writ- 
ten by Sparrow." 

The experimental evidence is available to support the supposition that 
the primary disturbance——on the biochemical level—elicited by ionizing ra- 
diations may be an aberration in the chain of metabolic events leading to the 
biosynthesis or utilization of nucleoproteins. Thus, Hevesy"® observed 
that x-irradiation inhibits the incorporation of C' carboxy-labeled ace- 
tate into desoxyribonucleic acid (DNA) purines. He stated that the find- 
ing ‘‘that the rate of formation of a main constitutent of cell nuclei is re- 
duced appreciably after irradiation helps to explain the blocking effect of 
irradiation on cell division and the formation of anomalous mitotic prod- 
ucts in the irradiated tissue.’’ Skipper and Mitchell’? demonstrated a simi- 
lar inhibition of incorporation of radioactive (C'*) bicarbonate and formate 
into viscera nucleic acids and nucleic acid purines and mice, 6 hours after ex- 
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posure to x-rays (950 r). Furthermore, the work of Abrams'* and of 
Holmes'’ shows that the rates of biosynthesis of nucleic acids are reduced 
markedly by x-rays, even though the synthesis of proteins is relatively un- 
affected. Also pertinent to this discussion is the observation of x-ray- 
induced depolymerization of thymonucleohistone and of Na thymonu- 
cleate.” *! Limperos and Mosher”? attempted to demonstrate this ef- 
fect in vivo. Rats receiving 1000 r whole-body irradiation, were sacrificed 
either immediately or 24 hours after exposure, and the thymus DNA was 
isolated and analyzed. The results indicated that a splitting off of either 
purine or pyrimidine bases had occurred in the DNA isolated from the rats 
immediately sacrificed; the thymus DNA isolated from animals sacri- 
ficed 24 hours after irradiation was almost completely depolymerized, as 
revealed by the structural viscosity, streaming birefringence, and sedimen- 
tation constants of the DNA solutions. Also of considerable significance 
to the present discussion is the experimental demonstration of nucleic 
acid depolymerization by chemically produced hydroxyl radicals,** and 
the evidence for an oxygen effect in the action of x-radiation on nucleic 


acid in vitro.”* 
The present finding of radiation protection activity in the spleen nuclei 


fraction and the destruction of this activity by treatment of the nuclei 
fraction with 1 / NaCl, or by x-irradiation of spleen homogenate in vitro 
appears to be in line with above observations when viewed in the light of 
the following hypothesis: It is postulated (a) that the cell nuclei of nor- 
mal mouse spleen contain a specific macromolecular nucleoprotein complex 
which is required for the cell growth and division, regeneration, or matu- 
ration of critical hematopoietic tissue (bone marrow, spleen); (6) that 
the biological activity of this nucleoprotein is inhibited by ionizing radi- 
ations. The inhibition may result from interference with the biosynthesis 
of this compound, or as a consequence of a physico-chemical alteration of 
the nucleic acid moiety which may result in depolymerization. Accord- 
ing to the above hypothesis, the exposure of mice to an LD, dose of total 
body x-irradiation, would result in a cessation or inhibition of the new- 
formation of the specific nucleoprotein complex. A deficiency in the supply 
or biological activity of this complex would then lead to a defect in the 
production or new-formation of leukocytes and erythrocytes and, since 
the biological half-life of the circulating leukocytes is known to be much 
shorter than that of the erythrocytes, the defect in hematopoiesis at the 
site of the formation would be manifested in the characteristic leukopenia 
of the acute radiation syndrome—the essential features of which can indeed 
be largely explained on the basis of cell damage leading to deficiencies that 
are first apparent in areas of rapid cell turnover, e.g., hematopoietic and 
intestinal epithelial tissues. 

The biological activity of the postulated nucleoprotein may, therefore, 
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be the limiting factor in the recovery of mice exposed to LDj total body 
x-rays; and the administration of the spleen nuclei fraction (obtained 
from normal, non-irradiated mice) to irradiated mice would provide a new 
source of the intact, non-degraded, biologically active nucleoprotein com- 
plex, able to stimulate or elicit hematopoietic recovery in the irradiated 
host mouse, and thus prevent death. 

The present observations that exposure of spleen homogenates to x- 
radiation (750 r) and to ultrasonic vibrations, in vitro, results in loss of 
radiation protection activity are consistent with the above hypothesis, 
since both these physical measures are known to induce degradation and 
depolymerization of nucleic acid macromolecules.” 

In view of the functional similarities betweeen the spleen and bone 
marrow in the mouse,” and the fact that injection of mouse bone marrow 
suspensions into irradiated mice affords them marked protection against 
death,” it seems likely that the protective factors in spleen and bone mar- 
row are closely related in nature, or even identical. This conclusion finds 
some support in the observation that the sensitivity of splenectomized 
mice to total body x-irradiation does not differ appreciably from that of 
intact, non-splenectomized mice.** Presumably the spleen and bone mar- 
row factors can replace each other. It is significant, too, that spleen and 
bone marrow, both tissues of high cell turnover, appear to be the major 
sources of the radiation protection factor in mice. In this connection it 
is of interest that Kelly and Jones*® were able to demonstrate an alteration 
in the turnover of desoxypentose nucleic acids in the liver, spleen and kid- 
neys of mice bearing growing mammary carcinoma transplants. These 
workers suggested “‘that the presence of a rapidly dividing tissue mass may 
influence the rate of cell division (as evidenced by the nucleic acid turnover 
rate) in the other body tissues.” 

Finally, it seems pertinent to note at this point that thiourea inhibits 
the depolymerization of DNA by x-radiation and by chemically generated 
free radicals, whereas ascorbic acid (another reducing agent) does not have 
this action.*’ Since thiourea protects rats against lethal x-ray doses when 


given prior to irradiation™ while ascorbic acid administration confers no 
radiation protection,*! it would appear plausible that preirradiation pro- 
tection by compounds such as thiourea or cysteine may be related to their 
ability to inhibit DNA depolymerization by x-rays. Thus, radiation pro- 


tection by these compounds would involve “protection’’ of the postulated 
spleen nucleoprotein factor against depolymerization and resultant loss 
of biological activity. Cronkite, ef a/., have previously suggested that 
the accelerated recovery of the bone marrow and of splenic myelopoiesis in 
glutathione-pretreated irradiated mice ‘“‘may be due to protection of hu- 
moral factors which control orderly hemopoietic regeneration.” 

The experimental data presented in this paper strongly suggest that the 
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post-irradiation protection factor present in the spleens of young LAf, 


mice is non-cellular in nature and is associated with the cell nuclei. Stud- 
ies aimed at the isolation and further identification of the factor are now 
in progress. 

Summary.—The subcellular particulate fractions of mouse spleen were 
isolated by differential centrifugation of spleen homogenates in sucrose 
media. The experimental data indicate that post-irradiation protection 
activity of mouse spleen homogenates is associated with the cell nuclei 
fraction. The highest degree of protective activity was exhibited by the 
nuclei fraction isolated in a modified, isosmotic sucrose-salt medium for- 
tified with adenosine triphosphate. In two of the experiments 100 per 
cent of the irradiated mice (750 r), each injected with the nuclei fraction 
derived from 280 mg of spleen tissue, were alive 30 days after exposure, 
whereas all of the irradiated controls (750 r) died. No activity was ob- 
served in the mitochondria, microsome, or soluble supernatant fractions 
isolated in this medium. A low level of protective activity was observed 
in whole spleen homogenates and in the spleen nuclei fraction prepared in 
0.25 M sucrose medium. Crude nucleic acid extracts prepared by treat- 
ing the nuclei fraction with | 17 NaCl were inactive. The data strongly 
suggest that the radiation protection spleen factor is non-cellular in nature. 
It is proposed that the factor may be a nucleoprotein. 
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THE MECHANISM OF PROPIONIC ACID FORMATION BY 
SUCCINATE DECARBOXYLATION. I. THE ACTIVATION 
OF SUCCINATE 


By H. R. WHITELEY* 
HOPKINS MARINE STATION, PACIFIC GROVE, CALIFORNIA 
Communicated by C. B. van Niel, June 5, 1953 


The studies of Johns'~* and Delwiche* with suspensions of anaerobic 


micrococei (Micrococcus lactilyticus or Veillonella gazogenes) and Propiont- 
bacterium spp. showed that these bacteria produced propionate by the de- 
carboxylation of succinate. This decarboxylation can also be accomplished 
by cell-free extracts of M. /actilyticus and the requirement for specific 
cofactors in this reaction has been briefly described.’ The present paper 
contains detailed evidence for the participation of these cofactors in the 
activation of succinate and for the nature of the active product. 
Methods.—M. lactilyticus (strain 221) was grown in tap water with 
2-3%% sodium lactate, 20% peptone, 1%) yeast extract, salts, and 0.0015 
thiamin. The latter was added routinely to all media, except when other- 
wise noted, because some batches of yeast extract, apparently deficient in 
thiamine, yielded small cell crops and extracts ef such cells had low decar 
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boxylase activity on a unit protein basis. When the deficient media were 
supplemented with thiamin, greater cell crops resulted, and extracts pre- 
pared from these cells had a high activity. 

Cultures were incubated at 37°C. for 16-24 hours; the cells were then 
harvested, washed once in tap water, ground with alumina,® the cell paste 
extracted with distilled water, or with a 0.01°% solution of a reducing 
agent, and centrifuged at approximately 12,000 X g for 20 minutes. Ex- 
tracts so prepared, with cysteine as the reducing agent, could be stored for 
several days in the frozen state without great loss in activity. 

Details of manometric procedure and of methods for the determination 
of volatile acids have been published elsewhere.’ The protein content of 
the extracts was estimated by the method of Stadtman, et al.* 

Total hydroxamic acid was determined according to Lipmann and 
Tuttle;® identification and quantitative estimation of the component 
acids was accomplished by paper chromatography.'® For the separation 
of succinhydroxamic acid from other hydroxamic acids a 10:2 mixture of 
phenol-water was used. 


TABLE 1 
THE Errect OF PH ON THE DECARBOXYLATION OF SUCCINATE 


ue. COs/10 MG. EXTRACT 
PROTEIN/HOUR 


190 
543 
560 
520 
512 
359 
199 


Pabst Co. preparations of adenosine-tri-phosphate (ATP) and coenzyme 
A(CoA, approximately 200 units/mg.) were used. Concentrations of 
CoA are expressed in terms of units as defined by Kaplan and Lipmann."! 

CoA and ATP were adsorbed from the extracts by treatment with Dow- 
ex-1 and Norit A;* 98% of the CoA, and 90-95% of the ATP content 
were thus removed. The former value was established by CoA assay 
with pigeon-liver enzyme,'! the latter by measurement of absorption at 
260 mu. 

Dialysis was performed at 5°C. against a 0.001% cysteine solution at 
pH 5.5; if continued too long, the extracts are permanently inactivated. 
Consequently the optimal duration of dialysis was established by deter- 
mining the activation of the extracts at frequent intervals. 

Results.-Conditions for testing the decarboxylation of succinate: All 


experiments were conducted under anaerobic conditions inasmuch as in air 


succinate is not decarboxylated but oxidized to acetate. The optimum 
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substrate concentration for carbon dioxide production is approximately 
20 uM /inl. with 20 mg. of extract protein. The highest rate of decar- 
boxylation occurred at pH 5.5 (table 1), and all experiments in which car- 
bon dioxide production was measured were performed at this pH. With 
freshly prepared extracts the rate of decarboxylation was approximately 
the same in various buffer solutions (acetate, propionate, benzoate, man- 


240 


A 
B 


MICROLITERS COo 








4 

60 

MINUTES 

FIGURE 1 
Decarboxylation of succinate by extracts of M. lactilyticus. A: Original extract, LO 
uM succinate. B: Supernatant fluid after centrifugation at 60,000 xX g, 10 uM 
succinate. C: Sediment after centrifugation at 60,000 X g, 10 uM succinate; also 

endogenous, original extract 


delate, and phthalate); citrate buffer was strongly inhibitory, and phos- 
phate caused a slight decrease in the rate of decarboxylation. Since MgCle 
partially relieved the inhibition by citrate, but not by phosphate, the 
action of the former was attributed to the binding of metal ions. The 
addition of MgCl and reducing agents increased the rate of succinate 
decarboxylation, the latter probably by maintaining the cofactors in the 
reduced (active) state. Cysteine, glutathione, thioglycollate, and 1,3- 
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dimercaptopropanol zt a final concentration of approximately 0.01 M 
were equally effective; NaS was toxic. NaF was added to inhibit an 
ATP-ase present in the crude extracts. Thus, carbon dioxide production 
from succinate by freshly prepared bacterial extracts was measured at 
30°C., in No», with the following reaction mixture: | ml. extract contain 
ing 10-20 mg. protein, 10 uM reducing agent, 10 uM MgCl, 100 uM NaF, 
90 uM succinate, and 50uM buffer, the choice of buffer depending on the 
nature of the experiment; total volume 2.2 ml. 

Figure | shows the production of carbon dioxide from succinate by an ex- 
tract of M. /actilyticus. Uf an extract prepared as described is further cen 
trifuged at approximately 60,000 * g for 40 minutes, a small amount of 
sediment is found. This sediment is not capable of decarboxylating suc- 
cinate, but the supernatant fluid is, and at a rate comparable to that de 
termined for the original extract. This shows that the enzymes involved 
in the decarboxylation of succinate are soluble. 

Suspensions of intact cells decompose succinate to equimolar amounts 
of carbon dioxide and propionate;' * this is true also for extracts. The 


decomposition of succinate is therefore a straight decarboxylation. f 


Cofactor requirements: Freshly prepared extracts, tested under the above 
conditions, decarboxylate succinate rapidly without the addition of any 
cofactors (Fig. 1). As described previously,’ an increase in the rate of car- 
bon dioxide production occasionally results from the addition of cocar- 
boxylase. This effect is restricted to extracts prepared from cells grown 
in a thiamine-deficient medium, or to extracts subjected to prolonged dial- 
ysis. 

Experiments with cell suspensions of biotin-deficient Propionibacterium 
pentosaceum™ showed that biotin increased the rate of succinate decar- 
boxylation. Since M. lactilyticus does not require this cofactor for growth,'® 
biotin-deficient cells could not be obtained. However, the addition of 
biotin to untreated or ‘‘aged”’ cell suspensions,'* or to freshly prepared, 
dialyzed, or depleted extracts did not increase the rate of succinate de- 
carboxylation. Adenylic acid and a factor prepared by autoclaving glu- 
cose with H.SO, have been reported to be involved in biotin-mediated re- 
actions; '® their addition, singly or in combination with biotin, did 
not affect the decarboxylation of succinate by the above-described prepara 
tions of M. lactilyticus. 

Freshly prepared extracts respond to the addition of ATP and CoA, al 
though a greater increase in decarboxylase activity upon the addition of 
these cofactors can be demonstrated with extracts that have been dialyzed 
or stored in the frozen state. The greatest effect is observed with prepara- 
tions from which CoA and ATP have first been removed by treatment with 
adsorbents. Such extracts have low activity, but almost complete 
reactivation can be accomplished by adding both ATP and CoA, in the 
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presence of cocarboxylase,® provided the pretreatment of the extract has 
not resulted in significant inactivation of the enzymes concerned. 

If traces of ATP remain in the extract after adsorption, the addition of 
CoA alone is sufficient to cause considerable carbon dioxide production.° 
However, a more complete removal yields extracts that do not respond to 
the addition of either cofactor alone; both ATP and CoA are required for 
succinate decarboxylation. The rate of carbon dioxide production by 
such extracts increases, though not proportionally, with increasing amounts 
of ATP. The experimental resuits, summarized in table 2, show that the 
addition of even small amounts of ATP causes the decomposition of large 
amounts of succinate in the presence of CoA. 


TABLE 2 
Tue Errecr or ATP AND COA ON THE DECARBOXYLATION OF SUCCINATE BY DEPLETED 
EXTRACTS 
no CoA AppDED 18 uNntrs CoA AappED 
uM ATP added 0 10.0 1.0 2.5 5.0 10.0 
ul. CO2/10 mg. extract 
protein /hour 46 40 226 234 338 380 


TABLE 3 
THe Errecr or ATP ON THE FORMATION OF SUCCINHYDROXAMIC ACID BY A DEPLETED 
EXTRACT 
uM HypDROXAMIC AciID/10 
uM ATP appep MG. EXTRACT PROTEIN/HOUR 
0 0.2 
5.0 L. 
10.0 3. 
15.0 4 
20.0 4.2 
10.0 3 
10.0 0.0% 


Reaction mixture contained: 1000 4M hydroxylamine, 500 uM TRIS buffer pH 7.3, 
50 uM succinate, 10 uM glutathione, 10 uM MgCh, 100 uM NaF; incubated at 30°C. 
under pure No. 

“ Hydroxylamine added at end of incubation period. 

» Reaction mixture heated for 5 minutes at 100°C. (pH 5.0) at end of incubation period 


and hydroxylamine added. 


Activation of succinate: The reactivation of depleted extracts by ATP 
and CoA indicates an initial conversion of succinate to succinyl-CoA or 
succinyl phosphate prior to decarboxylation. This inference is supported 
by the fact that a hydroxamic acid, identified as succinhydroxamic acid by 
chromatography, is formed in extracts to which hydroxylamine has been 
added. ‘The following observations suggest that, under certain conditions, 
the substance reacting with hydroxylamine is succinyl phosphate rather 


than succinyl-CoA. 
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The addition of hydroxylamine at the start of the incubation period com- 
pletely suppresses carbon dioxide evolution; the amount of succinhydrox- 
amic acid formed under these conditions does not depend on the quantity 
of CoA present,’ but is solely determined by the amount of ATP (table 3). 
The addition of hydroxylamine at the end yields only traces of hydroxamic 
acid if sufficient CoA is present; in this case, decarboxylation of succinate 
occurs. With extracts supplied with ATP, but unable to decarboxylate 
succinate owing to the virtually complete removal of CoA, nearly the same 
amount of succinhydroxamic acid is obtained when hydroxylamine is added 
at the end, as when it was present from the start. Since the quantity 
of activated succinate produced in the absence of hydroxylamine by depleted 

TABLE 4 
THE EFFECT OF VARIOUS SUBSTANCES ON CARBON DIOXIDE PRODUCTION AND SUCCIN- 
HYDROXAMIC ACID FORMATION 


“> OF CONTROL 
ADDITION® . HYDROXAMIC ACID 


Potassium arsenate § l 
Sodium chloride 63 
Lithium chloride 70 
Potassium phosphate 130 
@ Final concentration 0.01 M. 
TABLE 5 
THE EFFect OF ARSENATE ON THE FORMATION OF SUCCINHYDROXAMIC ACID BY A DE- 
PLETED EXTRACT 


uM ARSENATE ADDED 
UNITS COA ADDED 20.0 50.0 100.0 


0 §: 1.56 1.43 1.33 
9 0.28 1.20 0.00 
18 0.10 0.00 
36 1. 50 0.05 
Values given as uM hydroxamic acid/10 mg. extract protein/hour; 104M ATP in all 
reaction mixtures; experimental conditions listed under table 3; hydroxylamine added 
at end of incubation period. 
@ Hydroxylamine added at start of incubation period 


extracts depends on the ATP concentration, and far exceeds the trace 
amounts of CoA present, it appears likely that succinyl phosphate rather 
than succinyl-CoA is the compound giving rise to succinhydroxamie acid. 

If the reaction mixture resulting from the incubation of a depleted ex- 
tract, succinate, and ATP is heated for 5 minutes at 100°C. (pH 5.0), the 
addition of hydroxylamine no longer vields hydroxamic acid (table 3). 
Hence, the precursor of the latter is heat labile, a characteristic of acyl- 
phosphates, but not of acyl-CoA compounds,” although Sanadi and 
Littlefield’ recently found that CoA esters may also be heat labile in crude 


preparations. However, in an aliquot of the \/. /actilyticus extract used 


for testing the stability of the hydroxamic acid precursor, acetyl-CoA 
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was not decomposed by heat, thus supporting the proposal that the former 
is succinyl phosphate, rather than succinyl-CoA. 

Supplementary evidence of a more indirect nature for the formation of 
a phosphate derivative comes from experiments with arsenate. The ad- 
dition of arsenate completely inhibits carbon dioxide evolution and the 
production of succinhydroxamic acid (table 4), provided CoA is present. 
With increasing amounts of CoA, the arsenate inhibition of succinhydrox- 
ainic acid formation also increases (table 5). This effect can be accounted 
for by assuming a reversible, enzyme-mediated, reaction between succinyl 
phosphate and CoA with the formation of phosphate and succinyl-CoA. 
Substitution of arsenate for phosphate in the back reaction yields CoA and 
succinyl arsenate, followed by an instantaneous hydrolysis of the latter. 
This is analogous to the mechanism proposed for the disappearance of 
acetyl phosphate in solutions with phosphotransacetylase, arsenate, and 
CoA." 

Experiments on the effect of Na*, Lit, and PO*,, on the activity of ex- 
tracts also support the view that succinyl phosphate is produced (table 
1). These ions differentially inhibit carbon dioxide production and the 
formation of succinhydroxamic acid. Na* and Li* inhibit both, but the 
addition of PO, causes an increase in succinhydroxamic acid production 
and a marked decrease in carbon dioxide evolution. It is known that 
Na* and Lit have an inhibitory effect on phosphotransacetylase;'* the 
action of PO*, is explicable by postulating the reversal of a reaction in which 
phosphate is an end-product, such as the above-mentioned reaction be- 
tween succinyl phosphate and CoA. 

The results of these experiments therefore show that succinate is ini- 
tially activated in cell-free extracts of M. lactilyticus. The ATP-dependent 
accumulation of a substance capable of forming succinhydroxamic acid 
suggests that under certain conditions succinyl phosphate is the ‘‘active”’ 
form of succinate; experiments with inhibitors imply that succinyl-CoA 
may be formed. Evidence has recently been obtained for the formation of 
the both compounds by heart muscle preparations during the oxidation of 
alpha ketoglutarate.*” *! 

Succinyl phosphate may be produced in the bacterial extract either from 
succinate and ATP, or by the initial formation of succinyl-CoA from succin- 
ate, ATP and CoA, followed by trans-succinylation. 
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couragement of Dr. C. B. van Niel in the course of this work. 
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+t The production of equivalent amounts of carbon dioxide and propionate from 
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succinate can be explained most simply as the result of succinate decarboxylation and 
this interpretation is used in the present papers. It is, however, recognized that ex- 
periments with purified enzymes may show that the process iayolves successive oxida- 
tions of four-carbon compounds and reductions of three-carbon compounds, 
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THE MECHANISM OF PROPIONIC ACID FORMATION BY 
SUCCINATE DECARBOXYLATION. II. THE FORMATION 
AND DECARBOXYLATION OF SUCCINYL-CoA 
By H. R. WHITELEY* 

Hopkins Marine Station, Pacific Grove, California 


Communicated by C. B. van Niel, June 5, 1953 


The effects of arsenate, phosphate, Na*, and Li*, described in the pre- 
ceding paper,' suggest that succinyl-CoA is enzymatically produced from 
succinyl phosphate and CoA in extracts of Micrococcus lactilyticus. The 
CoA requirement for succinate decarboxylation thus indicated that suc- 
cinyl-CoA might be the actual substrate for the decarboxylase. Evidence 
is presented here for the decarboxylation of succinyl-CoA and the mecha- 


nism of its regeneration in such extracts. 
Methods.—-Sulfanilamide acylation was tested by the method of Kaplan 


and Lipmann,’ sulfanilamide being determined according to Bratton and 
Marshall.’ 
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Pyrophosphate formation was investigated under the conditions described 
by Lipmann, ef al.;‘ pyrophosphate was estimated by acid hydrolysis,’ and 
enzymatically with yeast pyrophosphates,® orthophosphate by the method 
of Sumner.’ Pyrophosphatase-free extracts were obtained by treatment of 
the crude extracts with the calcium phosphate gel which does not absorb 
pyrophosphates,® and elution of the gel with buffer at pH 5.7. 

Succinoxidase was purified by the method of Keilin and Hartree.’ Ace- 
tyl phosphate was synthesized from isopropenyl acetate and phosphoric 
acid;'’ suecinyl-CoA and propionyl-CoA, according to Simon and Shemin.'! 

Results. Decarboxylation of succinyl-CoA: Extracts depleted of ATP 
and CoA by treatments with adsorbents are unable to decarboxylate 
succinate.! However, the addition of succinyl-CoA to such extracts re- 
sults in a rapid evolution of carbon dioxide (table 1, exp. 1). The fact 
that succinyl-CoA decarboxylation proceeds at a much greater rate than 
that of succinate in the presence of free CoA precludes the possibility of 
attributing the results to a preliminary decomposition of succinyl-CoA by 
a deacylase. 

TABLE 1 
CaRBON DrioxipE PRopucTION BY DEPLETED EXTRACTS FROM SuccinyL-CoA AND 
SUCCINATE WITH VARIOUS ADDITIONS 
EXP. | EXP. 2 EXP. 3 
10 uM 50 uM 50 uM 
ADDITIONS SUCCINATE SUCCINATE SUCCINATE 
None 36 61 20 
10 uM ATP 35 70 22 
18 units CoA 58 87 38 
360 units CoA? 82 110 a 
10 uM ATP, 18 units CoA 143 180 141 
10 uM suceinyl-CoA‘ 2325 
1 uM propionyl-CoA “ae 244 
4 uM acetyl phosphate 
1 uM acetyl phosphate, 18 units CoA 
Values given in wl/10 mg. extract protein/hour. 
“ Amount of CoA used in preparing propionyl-CoA. 
» Succinate omitted. 
Bicarbonate-free preparation. 


The decarboxylation of succinyl-CoA may be represented either asa 
process yielding directly carbon dioxide, propionate, and free CoA, or as 
a reaction in which one of the end-products is released as a CoA compound. 
The latter possibility was examined by incubating succinate, ATP, CoA, 
and an amount of hydroxylamine insufficient to cause a complete inhibi- 
tion of succinate decarboxylation, with an extract capable of activating 


succinate but virtually incapable of activating free propionate. Chroma- 


tography of the reaction mixture disclosed the presence of both succin- and 
propionhydroxamic acids, and comparison of the quantity of the latter 
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formed from succinate with that produced in a similar reaction mixture 
containing free propionate instead of succinate showed that the amount 
produced from succinate was by far the greater (table 2). Of special sig- 
nificance is the close correspondence between the quantities of propion- 
hydroxamice acid recovered and carbon dioxide liberated during the decom- 
position of succinate. Under the experimental conditions used, the action 
of CoA transphorase, described below, cannot be detected. These results 
strongly suggest that propionyl-CoA is a direct product of the decarboxyla 
tion of succinyl-CoA. 
TABLE 2 
QUANTITIES OF PROPIONHYDROXAMIC ACID FORMED BY A DEPLETED EXTRACT BY THE 
ACTIVATION OF PROPIONATE AND THE DECARBOXYLATION OF SUCCINATE 


5 uM PROPIONATE 50 wM SUCCINATE AS SUBSTRATE 
AS SUBSTRATE uM PROPION 
INCUBATION PERIOD, uM PROPION HY DROXAMIC 
MIN HYDROXAMIC ACID ACID uM COs 


30 0.038 0.55 0.58 
45 0.06 0.77 0.82 
60 0.20 0.88 0.93 


Reaction mixture contained: 800 4M hydroxylamine, 50 uM phosphate buffer pH 
7.3, 10 uM glutathione, 10 uM MgCl, 100 uM NaF, 10 uM ATP, 36 units CoA; incu 
bated at 30°C. under pure Ne. Values expressed in u4M/10 mg. extract protein/hour. 


TABLE 3 
ACYLATION OF SULFANILAMIDE BY A DEPLETED EXTRACT 
SULFANILAMIDE 
ACYLATED, 
uoG./10 MG. EXTRACT 
uM ATP appEp UNITS COA ADDED PROTEIN/HOUR 


0 0 
0 18 
10 0 
10 18 l 
10 36 I 
20 18 2 


Reaction mixture contained: 504M TRIS buffer pH 8.0, 0.2 ml. pigeon liver prepara 
tion, 50 uM succinate, 100 wg. sulfanilamide, 10 wM glutathione, 100 wM NaF, 10 «6M 
MgCl; incubated at 380°C. under pure No. 


Reversibility of the decarboxylation reaction was tested by incubating 
an extract with propionate and bicarbonate. The formation of succinate 
was revealed by the reduction of 2,6-dichlorophenol indophenol under the 
influence of partially purified succinoxidase. Addition of malonate in- 
hibited dye reduction, as did omission of bicarbonate from the mixture. 
These observations eliminate oxidation of propionate or other substances 
present in the extract as the cause of the reduction. On the other hand, 
addition of ATP and CoA greatly increased the reduction; this was not 
the case when succinate was used as substrate. 
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Formation of CoA compounds in extracts: Because the evidence for 
suceinyl-CoA production by phosphotrans-succinylase! is indirect, an 
attempt was made to determine whether acylation of sulfanilamide, con- 
sidered specific for acyl-CoA derivatives, can occur in extracts provided 
with ATP, CoA, succinate, sulfanilamide, and ‘‘acylating enzyme.’’? 

The results of these experiments show unequivocally that sulfanilamide 
acylation occurs (table 3), though only to a small extent. If succinyl-CoA 
is the acylating agent, this is readily understandable since this substance 
is also rapidly decarboxylated; the acylating enzyme must therefore com- 
pete with decarboxylase for its substrate. In agreement with this expla- 
nation is the fact that dialyzed extracts, with reduced decarboxylase ac- 
tivity, cause a larger sulfanilamide acylation. Moreover, in a reaction mix- 
ture incubated at pH 8.0, considerably above the pH optimum for decar- 
boxylation, more sulfanilamide disappears, although this may be partly 
the result of more favorable conditions for the action of the acylating 
enzyme. 

Both ATP and CoA are essential for sulfanilamide acylation; a small 
amount of CoA suffices for maximum effect, but a substantial increase in 
sulfanilamide disappearance is caused by augmenting the ATP supply. 
Regeneration of CoA as a consequence of the transfer of the acyl moiety of 
an acyl-CoA compound to sulfanilamide readily accounts for the former; 
the ATP effect is discussed in more detail below. 

Mechanism of succinyl-CoA formation: Small amounts of ATP and CoA 
in extracts of M. lactilyticus permit the decomposition of large quantities 
of succinate to propionate and carbon dioxide.' If this process involved 
the decarboxylation of succinyl-CoA with the liberation of free CoA, the 
need for only catalytic amounts of the coenzyme would be understandable. 
However, the situation with respect to ATP is different; succinyl-CoA 
formation exclusively by way of succinyl phosphate, or by some other 
mechanism involving both ATP and CoA, should require a quantity of 
ATP equivalent to that of succinate decomposed, since it is difficult to 
conceive of a regeneration of ATP from ADP, or AMP, and inorganic phos- 
phate coupled with a simple decarboxylation. 

Because traces of ATP suffice for an extensive succinate decomposition, 
a mechanism of succinyl-CoA formation independent of ATP appeared 
probable. Such a mechanism would be provided by a decarboxylation of 
succinyl-CoA with the liberation of an acyl-CoA compound which could 
react with free succinate to regenerate succinyl-CoA through the media- 
tion of a CoA transphorase. A reaction of this type has been postulated to 
account for the activation of propionate'* and formate’ by extracts of 
Clostridium kiluyveri in the presence of acetyl phosphate and CoA. 

Formation of propionyl-CoA by the decarboxylation of succinyl-CoA 
seemed probable; hence the occurrence of a transphorase reaction between 
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propionyl-CoA and succinate was investigated. Depleted extracts, able 
to decarboxylate succinyl-CoA but not free succinate, were incubated with 
succinate and a small amount of propionyl-CoA. As seen from Table 1, 
exp. 2, rapid decarboxylation occurred, and the rate of carbon dioxide 
evolution exceeded that from succinate in extracts supplemented with ATP 
and CoA. 

In view of these results it seemed possible that other CoA compounds, 
e.g., acetyl-CoA, would likewise serve as CoA donor. Acetyl-CoA was 
not added as such, but generated in the extract from acetyl phosphate and 
CoA through the action of phosphotransacetylase whose presence had been 
established by tests on the arsenolysis of acetyl phosphate and acetyla- 
tion of sulfanilamide. Table |, exp. 3, shows that depleted extracts can 
decarboxylate succinate if supplemented with acetyl phosphate and CoA. 

The following experiment provides more direct evidence for suceinyl- 
CoA formation froin succinate, acetyl phosphate, and CoA. Depleted 
extracts, incubated with succinate, acetyl phosphate, and CoA, yield three 
hydroxamic acids upon the addition of hydroxylamine; these have been 
identified as acet-, propion-, and succinhydroxamic acids. The latter is 
formed by the above reaction mixture only if CoA is added, and since in 
the presence of both acetyl phosphate and CoA, succinate is decarboxylated, 
the formation of activated propionate could be expected. It should be 
noted that if hydroxylamine is added at the start of the incubation period, 
CoA transphorase action cannot be demonstrated and only acethydroxamic 
acid is found in the above reaction mixture. 

Incubation of depleted extracts with acetyl phosphate, CoA, and propi 
onate, followed by the addition of hydroxylamine, yields propionhydrox 
amic acid, as well as acethydroxamie acid, showing that CoA transfer can 
also occur between acetyl-CoA and propionate. Activated formate is 


produced by substituting formate for propionate in the above reaction 


mixture. The quantity of succin- and propionhydroxamic acids found 
in these experiments is far greater than the available amount of CoA, 
suggesting that phosphate derivatives are formed from the CoA compounds 
by the action of phosphotrans-succinylase and -propionylase. 

The effect of propionyl-CoA and acetyl-CoA on the decarboxylation of 
succinate, and the evidence for the formation of activated succinate from 
succinate and acetyl-CoA indicate that generation of succinyl-CoA by 
CoA transfer can occur. Hence a decarboxylation of succinyl-CoA with 
the production of an acyl-CoA compound provides a simple mechanism 
for the continued formation of succinyl-CoA independent of ATP. The 
experiment summarized in table 2 permits the conclusion that propionyl 
CoA is this compound, rather than a one-carbon CoA compound." 

In the absence of acyl-CoA compounds, the initial formation of succinyl] 
CoA must, of course, proceed by other means. That succinate can be 
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decarboxylated in the presence of CoA and ATP shows that this can be 
accomplished through a reaction in which ATP takes part, either through 
the primary formation of succinyl phosphate,! or by a different type of 
reaction involving both ATP and CoA. Such a reaction could result in 
the concomitant formation of ADP and inorganic phosphate, as described 
for the activation of succinate by heart preparations,” or of AMP and pyro- 
phosphate, as described for the activation of acetate by yeast preparations.‘ 
Since pyrophosphatase-free extracts of 7. lactilyticus, incubated with suc- 
cinate, ATP, CoA, and hydroxylamine‘ showed an increase in inorganic 
phosphate and gave negative tests for pyrophosphate, it must be concluded 
that an appreciable quantity of succinyl-CoA is not formed by the latter 
activation mechanism. 

The results of the experiments to date indicate that the decomposition 
of succinate could be schematically represented by the following diagram: 


Succinate 


ATP + CoA 


' 
' 
! 
' 
' 


Succinyl-CoA + Propionate 
Tt 7 
; 
Propionyl-CoA + Succinate 
+ 
CO» 

Theoretically the formation of an infinitely small amount of succinyl- 
CoA suffices to initiate the decarboxylation. In the presence of hydrox- 
ylamine, or of sulfanilamide and an “‘acylating enzyme,’ the acyl-CoA 
compounds of the cycle are partly drained off. Thus the cycle is inter- 
rupted, and larger quantities of ATP are required for the primary synthesis 
of succinyl-CoA. 

The diagram indicates that succinate formation by a reversal of the de- 
carboxylation reaction would require carbon dioxide and propionyl-CoA. 
However, this synthesis of succinate may be considered to be analogous to 
that of acetoacetate from acetate. It has been established that the latter 
requires two “‘active’’ acetate units, in the form of acetyl-CoA."* Hence 
it appears reasonable to suggest that the synthesis of succinate, too, may 
be possible only if both propionyl-CoA and an activated one-carbon unit 
are available. This view receives support from the experiments of Del- 
wiche, ef al.'* The effect of ATP and CoA on the synthesis of succinate by 
extracts of M. lactilyticus might be attributed te their involvement in the 
production of both these ‘‘active’’ building blocks. 
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GENERALIZED MATHEMATICAL RELATIONSHIPS FOR 
POLYPEPTIDE CHAIN HELICES. THE COORDINATES OF 
THE WU. HELIX* 

By BARBARA W. Low AND H. J. GRENVILLE-WELLS 
UNIVERSITY LABORATORY OF PHySICAL CHEMISTRY RELATED TO MEDICINE AND PUBLIC 
HEALTH, HARVARD UNIVERSITY, BOSTON, MASSACHUSETTS; AND LABORATORY FOR 
INSULATION RESEARCH, MASSACHUSETTS INSTITUTE OF TECHNOLOGY, CAMBRIDGE, 
MASSACHUSETTS 
Communicated by J. T. Edsall, June 4, 1953 


The formulation of polypeptide chain configurations is concerned 
principally with intra-chain packing and _ interactions between 
atoms in neighboring residue units of the skeletal backbone chain 
(—C:-CO-NH-—),. The folded or coiled configurations are held together 
by intra-chain (CO-HN) hydrogen bonds between residues. In helical 
configurations the residues are structurally equivalent; they all have 
identical intra-chain environments (exclusive of amino acid residue varia- 
tions). Such configurations, therefore, effectively provide maximal 
intra-chain hydrogen bonding. 

Detailed studies depend on the use of precise dimensions for the poly- 
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peptide chain. Pauling and Corey have predicted as the fundamental 
unit a chain residue with dimensions and configuration derived from 
detailed studies of crystal structure analyses of amides, amino acids, 
and peptides.' On the same basis they have also defined appropriate 
limits for the hydrogen-bonded oxygen to nitrogen distance and for the 
angle between the NH and NO vectors. 

Only two helical®: * configurations, the a and y helices, fulfill completely 
all of the requirements specified by Pauling and Corey. A third helix, 





Cfra"a" 


FIGURE 1 


Orientation of the chain residues in the ‘‘a’’ and ‘‘y” series. Deri- 
vation of the x and z¢ coordinates for the atoms of the first residue. 
The x and z axial directions in each series are perpendicular and paral- 


lel, respectively, to the helical axis. 


the m helix,! involves so small a deformation (4°) in the intra-chain a 
carbon bond angle that it also should probably be considered to fulfill 
their requirements. Deformations of this order have been observed in 
amino acid structures.® 

The mw helix was discovered while working with three-dimensional scale 
models which had dimensions slightly different from the Corey-Donohue! 
dimensions first used by Pauling and Corey. In an attempt to determine 
the best set of atomic coordinates for this structure, an analytical method 
was developed. The procedure establishes certain relationships between 
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different polypeptide chain helices which appear to have been confused in 
the past. The method, is rapid and quite general; 1t may also be useful 
in the study of non-peptide helical structures where there exists a choice 
among different values for the independent variables. 

In the three helices, a, 7, and y, which satisfy the Pauling-Corey re- 
quirements, each residue (~-C-CO+-NH-+C—-) is coplanar and lies in one of 
a zone of planes parallel to and at a common distance from the helical 
axis. The helices may be described by specifying the carbonyl oxygen 
to which each NH group is hydrogen bonded. Thus, in the m helix there 
is a hydrogen bond between one amide group and the fourth amide group 
beyond it along the chain. This description, however, is inadequate 
since it overlooks the two fundamentally different ways of coiling a 
helical structure. 

In figure 1 the fundamental dimensions for a_ single residue 

*Cy-CiO;- Ni Hy Coe have been drawn to scale for the most recently 
revised set of values® and inscribed in a circle of diameter C,;--C». The 
residue lies in a vertical plane and if one end (C;) is kept fixed, the residue 
may be tilted either up or down. If P is the pitch of the final coil, the 
vertical translation per complete turn, and m is the number of residues 
per complete turn, then in a given helix this vertical translation per residue 
is P/n. If the residue is tilted upward so that C,—C) (solid line) is its 
horizontal projection (normal to the helical axis), then the CO and NH 
groups are almost vertical and parallel to the helical axis. A P/n value 
of 1.14 A was used to construct this figure and over the range P/n = 
0.9 — 1.5 A deviations from the vertical are not very large. Tan @ = 
P/nL where L = C;—C3. The a and 7 helices are both constructed in 
this way and we have considered all such helices as belonging to the 
“a” series. 

When the residue is tilted downward so that C; -C3 (broken line) is 
horizontal then the CO and NH groups are both sharply inclined to the 
helical axis. The C, c and N,-—-C., bonds are practically horizontal 


a * bd a hs . . . 
and parallel to C;—-Cy, over a wide range. The y helix is constructed in 


this way and we have considered all such helices as belonging to the ‘“y”’ 
series. 

The residue-unit configurations differ markedly between the two series. 
Direct transformations between members of different series is impossible. 
Thus, the early suggestion’ that the y-helix might represent supercontracted 
keratin and supercontraction corresponds to a direct transformation 
a — y¥ is theoretically invalid without reference to the comparative in 
stability of either configuration. 

Figure 2 (a) shows a perspective drawing of the mw helix with 4.4 turns 
and a pitch P = 5.02 A and hence a unit translation per residue P/n = 
1.14 A. Each succeeding residue, C\(CO»-NH)Co, C.(CO-NH)C;, ete., 
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lies in a vertical plane parallel to the helical axis and as the plan view, 
figure 2 (b), shows the angle between the two planes is the helical angle 
¥ = 369/n. The angle 180 — y is the horizontal component of the tetra- 
hedral angle N,C.C,’ and so the tetrahedral angle depends on both the 
number of turns (m) in the configuration and the vertical unit residue 
translation. Both the plan view (in which the nitrogen positions have 
been omitted at the cross-over points) and the drawing, illustrate that the 
CO and NH groups are almost vertical and the angle between the NH and 
NO vectors is small. 

The m helix shown has a right-handed screw axis. The mirror image 
left-handed helix corresponds to the mirror image of figure 1 with C,-—C, 
plotted from right to left. Looking at the helix upside down corresponds 
to a horizontal projection C}—C» (not drawn) parallel to C,—C3. 

As a helix is wound up, a residue may theoretically form hydrogen 
bonds with the mth amide group beyond it. Figure 2 shows that for the x 
helix this corresponds to the chemical sequence, O 

| 

C—(NH-C(HR)+CO), N 
read in either direction. For the ‘‘a’”’ series the sequence is therefore 
i Pre 


I 
C—(NH-C(HR)-CO), N with 3n + 4 atoms in the hydrogen-bonded 
loop so formed. Bragg, Kendrew and Perutz* have identified helices by 
symbols which indicate the number of turns in a unit repeat (screw sym- 
metry) with subscript giving the number of atoms in the hydrogen-bonded 
loop. The m helixt would thus be written 4.416. We have avoided this 
identification where possible, for it leads to ambiguities; the helices 
3.6,3 and 3.7), are both the a helix,® whereas 3.64 is a helix of the ‘‘y’’ series. 

Variations in the number of turns per unit repeat are correlated with 
changes in the ‘‘tetrahedral’’ angle and hydrogen bond length and angle. 
Minor revisions are therefore likely to be made from time to time. The 
right-handed helix wz in the drawing is shown with two of the 8 carbons 
of the amino acid residues drawn in (8C2, BC,) at the correct position for 
naturally occurring L amino acids." 

The ‘‘y"’ series residue unit in figure | would lead to left-handed helices. 
In winding up these units so that a residue may form hydrogen 
bonds with the mth amide group beyond it, the chemical sequence 


Fig. 2 (a) —Perspective drawing of the right-handed ap helix. The link-end carbon 
The 8 carbon and hydrogen atom positions for L- 


atoms are numbered consecutively. 
Plan view of mp helix. 


amino acids are shown attached to Cs(a@) and Cy(@). Fig. 2 (b) 


The hydrogen positions are shown with broken lines. Nitrogen positions are omitted 


at the two cross-over points. 
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NITROGEN 


CARBON 


OXYGEN 


HYDROGEN 





FIGURE 2 (a), ABOVE; 2(6), BELOW 
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H ies © is obtained. There are 3(n — 1) + 5 


N (C(HR).CO.NH),—:C(HR) C 
atoms in the hydrogen-bonded loop so formed. 

Atomic coordinates for both series of helices are derived as follows. 
If the origin of the rectangular coordinate system is taken at C, with the 

P ‘ .* ‘ - . 

x axis along C,;—-Cy and the z axis parallel to the axis of the helix, the 
coordinates of atoms in the first residue are calculated from equations of 
the kind: 


0) 

C,O; cos (@ + 4) 

CC; cos (6 + 62) (1) 

C,C, cos ¢ = L (projected length of the residue) 

0 

C,O; sin (@ + 4;) 

CC} sin (@ + 4) 

CC. sin @ = P/n (2) 
In the ‘“‘y’’ series the z axis direction is changed. The same expressions 
hold except that the signs of the angles 6,, 6, 6; and 4, are reversed. For 
all atoms in the first residue y = 0 (the residues are both planar and 


vertical). 
‘rT 4 4 «/ s 
Phe lengths C,O,, CiC), and angles 6, #2, etc., are constants for a given 


residue model and may be calculated directly. The x and z coordinates 
for the atoms of the first residue unit are functions of the vertical rise per 
residue, P/n. In figure 3 (a) and (b) the atomic coordinates are plotted 
for the ‘‘a’’ and “‘y”’ series of helices, respectively, over the range of values 
of P/n from 0.95 to 2.0. They were prepared using Pauling and Corey’s 
most recent set of fundamental residue dimensions. For example, the 
broken vertical line gives values for @ helix with P/n = 1.48. Similar 
expressions for the Corey-Donohue configuration are shown graphically 
in figures 4 (a) and (6). There are some marked differences near the 
limits of the range explored even though the dimensions of the two models 


are closely similar. 

Planes containing successive linked residues are equidistant from the 
helical axis. The a carbon atoms C)CoC,, ete., lie, therefore, on the surface 
of a cylinder about the helical axis of radius R = L/(2 sin y/2). In 
horizontal projection they lie, therefore, on the circumference of a circle 
(Fig. 2 (6)). The coordinates of the rth link-end position are: 


xc, = L(O + 1 + cos ¥ + cos 2y ... + cos (r — 2)¥) (3a) 


Fig. 3 (a)—Atomic coordinates for the ‘‘a@’’ series atom positions plotted as a 


function of P?/n. (Pauling-Corey dimensions.*) 
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yor = L(O+ 04 sin py + sin 2y .,. + sin (r — 2)yp) (3b) 
= (r—1)P/n (3c) 
The coordinates of the point P, in the rth link is given by: 
Xpp = Xer + Xp,scos (r — 1)y¥ 
Yer = Yor + Xp,ssin (r — 1)y 
sp, = (r — 1)P/n + 2p, 
The equations (3a) and (4b) may be rewritten: 
Xcr = Leu, (5a) 
Vc, = (5b) 
where u, and v, are the unitary coordinates of the link-end positions. 
u, = (1 + cosy + cos 2y ... + cos (r — 2)p) (6a) 
v, = (sin yy + sin 2y ... + sin (r — 2)y) (6b) 


In figure 5 the unitary coordinates for successive a carbon atoms are 
recorded for r values of 2 to 6 over the range n = 3.0 to 5.5. 
The equations (4a) and (4b) may be rewritten: 


xp, = L.u, + A;*Xp, (7a) 
yp, = Lev, + B,+xp, (7b) 


where A, = cos (r — 1)W and B, = sin (ry — 1)y. Values of A, and B, 
forr = 2 to 6 are given in figure 6 plotted against 7 the number of residues 
per complete turn. 

Atomic coordinates of the atoms N,, Co, c are required in order to 
determine the ‘‘tetrahedral’’ angle. The direction cosines of the two 
vectors N,;C, and C.C} may be readily evaluated. The ‘tetrahedral’ 
angle N,C.C) = x is given by the usual cosine expression. This simplifies 
to 


j(Aoexe,) (L — xn,) + to,°(P nm — &y,)I 
COh XY Moe ee ee ee ee (S) 
| Ni Coe CoC, f 


Similarly the angle + between the N,—-O,; and N,—-H, vectors is given by 


the equation: 


Fig. 3 (b)-Atomic coordinates for the ‘‘y"’ series atom positions plotted as a func- 
tion of P/n. (Pauling-Corey dimensions.*) 
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CoS T = 
(X13, — X0,)(XNr — XH,) + yn-(ynr — YH,) +(2nr — 20,)(2Nr — 2H,) (9) 
N,H,-N,O, 


where 


N, -— Oy = V (xn, — Xo,)? + (ynr)? + (2nr — Zo)? (10) 


If the helical axis is taken as the origin of the coordinate system (the 
most convenient form for calculation of the radial distribution function) 
cylindrical coordinates* (p, 6) of the atoms may be readily evaluated. 


Thus for example 


pn, = WV (xy)? + R? — an, (11a) 


‘ ae (= y *) (11b) 
PN, 

Coordinates for the 6 carbon atom may be easily derived once the N,C:C; 

angle has been evaluated. 

Atomic coordinates for two variations of the m helix are given in table 1. 
When the helix is constructed with 4.3 residues per turn the strain in the 
tetrahedral angle is minimized, but the O,N;H; angle is larger than that 
shown in figure 2. Series of atomic coordinates have been calculated using 
(a) the original scale model dimensions," (6) the Corey-Donohue! dimen- 
sions, and (c) the Pauling-Corey dimensions.® From these, values of the 
tetrahedral angle and hydrogen bond dimensions have been derived. For 
values of P/n = 1.14 to 1.15 and n = 4.3 to 4.4 (1.e., P = 4.9 to 5.06), 
the following range of values was obtained. 


-oO 


ZN,C,C; 113.5° — 115.5 
O,.-N, 2.75A~2.83A 
ZONsHs 3.5° =» 18° 


No attempt has been made to develop differential functions; however, 
the x and z coordinates in the graphs, figures 3 and 4, are the differen- 


tials of each other. 


Acosx = —Asinx 
x 
For a given residue configuration the pitch P and the number of residues 
per turn, m, may be varied independently. Once a set of values has been 
chosen it is obvious what type of adjustment is needed to improve the 
solution. A graphical method!? has also been devised which permits 
extremely rapid trial and error determinations of the three critical variables. 
Four helices of the ‘‘a’’ class have been described. ‘‘a’’,-: corresponds 
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to the 2; or a, configuration proposed by Huggins,'* Zahn,'* and Mizu- 
shima.'© Examination of figure 1 shows that in this helix the O,N»H» 






angle will be Jarge for the values of P/n required. ‘‘a’’,<2 corresponds 
to the 3.0j helix suggested by Taylor'® and Huggins. ‘‘a@’’,.3 is the a 







rABLE 1 






ATOMIC COORDINATES, AND VALUES FOR N.C;C,, Ns—O;, AND O.NH;; x, yandcin A 











y 2 







x 










n = 43, P = 4.90, P/nm = 1.14 
Ci 0 0 0 














C, 1.18 0) 0.98 
) QO, 0.99 0 2. 20 
i Ni 2.36 0 0.39 
H, 2.5] 0 ~0.60 

Cy 3.62 0 1.14 
' Helical axis 1.81 2.02 


C, 










ron 1.18 0 0.98 
O; 0.98 0 2 91 
Ni 2.36 0 0.40 
H, 2.51 0 —0.59 
, C; 3.62 0 1.15 
Helical axis 1.81 2 09 
ZNiC.C, 115 
N;—O, 2 SOA 





ZO NSH; 











C, 















_£ 1.23 
O; 1.10 0 2.15 
Ni 2.38 0 0.27 
| H, 2.48 0 —0.73 
Cc 3.68 0 0.95 
Helical axis 1.84 2 66 
| ZN,C roy 122 
: N;—O: 2.93 A 
' ZO\NH 17° 








helix and ‘‘a’’,-4 is the mw helix. ‘“‘a’’,<5 has been studied by us. The 
horizontal component of the N,C.C, angle is 180-360/(5 + 6) where 6 is 
fractional, and is, therefore, itself of the order of 110°. In the ‘‘a”’ series 
the vertical component of this angle is considerable. For n = 5.2, P = 
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1.95, and the Pauling-Corey model, figure 1, the N,C.C; angle is 122° 
(table 1). This deformation is large. However, as the range of values for 
the helix shows, slight changes in the model and in the values for P and n, 
may lead to considerable differences in the tetrahedral angle derived. 
Hydrogen bond length and angle are reasonable for the structure. All 
the members of the group “‘a’’,.. > “‘a’’,=5 are acceptable if the restrictive 
criteria are somewhat relaxed. Transformation between different con- 
figurations may be accomplished by simultaneously decreasing or increas- 
ing the wind up angle, and either reducing or increasing the vertical tilt 
per unit residue. 

Four helices of the ‘‘y’’ series have been described. In this series the 
angle N,C.C;3 is almost equal to its horizontal component ISO — yp. “‘y’’»=2 
(2s)'* and ‘‘y’’,<3 (the 11 atom ring)" are therefore completely excluded 
for the coplanar configuration. In “‘y’’,.4 this angle is still rather small. 
For the variant with 4.3'* residues per turn the horizontal component is 
96.3°. ‘“y'',=5 is the y helix. ‘“‘y’’,<. and higher members of this series 
give ‘‘tetrahedral’’ angles which are much too large. 

The relative stabilities of polypeptide chain helices have been evaluated 
in terms of structural features which are not considered here. Pauling 
and Corey” have estimated, for the a and y helices, the strain involved in 
possible “unfavorable orientation’’ about the N,C. and Cc bonds, and 
incomplete van der Waals’ stabilization.’ Robinson and Ambrose” and 
more recently Donohue'* have both made further estimates of the quali- 
tative or semiquantitative agreement between certain helices and a wide 
range of restrictive criteria. The fundamental residue dimensions used 
are not specified in detail for either study. Both groups of investigators 
have used models different from the Pauling-Corey residue unit, and 
considerable latitude, including the introduction of slight distortions from 
the planar amide configuration, has been permitted in adjusting the 
residue dimensions to improve the agreement with their theoretical re- 
quirements. 

The graphs reproduced in this paper permit the rapid variation of 
helical configurations based on a particular model, but must be redrawn 
for each new residue model. The amount of variation to be expected 
is illustrated by comparison of Figs. 3 and 4. These graphs are for a 
residue lying wholly in one vertical plane parallel to the helical axis. 
While this restriction can be easily relaxed in any particular case, once 
values of P/n and n have been chosen, systematic allowance for lack of 
coplanarity would involve the construction of families of graphs. This 
elaboration did not appear justified at this time. 
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CONDUCTANCE AND VISCOSITY OF CONCENTRATED 
AQUEOUS SALT SOLUTIONS AT 50.5° 


By M. Joun Ricr, Jr.,* f AND CHARLES A. KRAUS 
METCALF RESEARCH LABORATORY, BROWN UNIVERSITY 


Communicated June 12, 1953. 


Our knowledge of the electrical conductance and the viscosity of con- 
centrated salt solutions is very limited. Campbell and Kartzmark' have 
measured the conductance and the viscosity of solutions of silver nitrate and 
ammonium nitrate in water at 95° up to saturation. Seward? has measured 
the conductance and viscosity of solutions of tetrabutylammonium picrate 
in butanol at 93° over the complete range from pure solvent to pure fused 
salt; Strong* has measured the conductance of some quaternary ammo- 
nium salts in benzene at 25° up to fairly high concentration but viscosity 
data are for the most part lacking. Moessen and Kraus‘ have measured 
the conductance of solutions of tetrabutylammonium bromide and tri- 
methylammonium chloride in bromine at 25° to fairly high concentrations 
but there are no viscosity data for these solutions. The same is true of 
potassium iodide in iodine at 140°C. which was measured by Lewis and 
Wheeler.’ 

Cesium formate has been reported to be soluble to the extent of one mole 
of salt to one-half mole of water at 50°.6 However, this datum is in error; 
the limit of solubility is one mole of salt to approximately two moles of 
water. However, potassium formate is soluble to the extent of approxi- 
mately one mole of salt to 1.1 mole of water at 50°. The conductance and 
viscosity of solutions of these two salts have been measured at 50.5° 
to near saturation. Solutions of the potassium salt were measured to 
lower concentration in order to approximate the value of Apo. 

In the preparation of cesium formate, cesium chloride was converted to 
cesium nitrate by repeated treatment with nitric acid. After repeated 
recrystallizations, the nitrate was converted to carbonate by treatment 
with oxalic acid in the presence of a small amount of water and subsequent 
calcination. The carbonate was converted to forinate by neutralization 
with formic acid and evaporating to dryness. The salt was recrystallized 
by dissolving it in 96% alcohol and adding ether to precipitate the desired 
proportion of salt. The final product was dried in vacuo at 150° for 4-5 
hours, m.p. 264°C. (corr.). The nitrate, when tested spectroscopically, 
showed no traces of metallic elements other than cesium. 

Potassium formate was prepared by neutralizing potassium carbonate 
with formic acid and evaporating to dryness. The salt was recrystallized 
from 96% alcohol on addition of ether. The product was dried in vacuo at 
150° for 8 to 10 hours. This material was found to contain 1.5% water. 
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The water could be completely removed by regrinding the salt and heating 
in vacuo at 150° for 10 hours. When necessary, operations were carried 
out in a dry-box. 

The solutions were made up in the cells in which their conductance was 
measured. Three cells of different constants were employed. Two of 
these were of the Erlenmeyer type with the electrode chamber attached to 
the outside of the flask. These cells were similar to those described by 
Daggett, Bair, and Kraus,’ although much smaller. The third cell con- 
sisted of two cylindrical tubes which were joined at the bottom by a 55-mm. 
length of smaller tubing of approximately 5 mm. diameter. The cylindrical 
cells were closed by ground glass caps which carried tubes at the ends of 
which were sealed platinum electrodes that normally projected into the 
solution. The cell constant of this cell was 38.430 and was independent of 
the precise setting of the caps or the volume of solution in the cell. This 
cell was provided with stopcocks so that the solution could be forced from 
one arm to the other by air pressure for the purpose of mixing. The other 
two cells had constants of 0.3081 and 10.669. 

Cell No. 1 was of approximately 300 cc. capacity and was used for meas- 
uring the most dilute solutions. A known weight of salt was introduced and 
weighed quantities of water were added successively. When the cell was 
full, the greater portion of the solution was withdrawn into a weight pipette 
and a new series of additions of water was begun. Cell No. 2 had a capacity 
of approximately 800 cc. and was used for intermediate concentrations. 
Measurements could be made with a minimum volume of 50 cc. of solution. 
The desired concentration range could be covered by merely adding known 
weights of water to a known weight of salt. Cell No. 3 had a capacity of 
approximately 100 cc. and measurements could be made with a very small 
volume of solution. The desired concentration range could be covered by 
merely adding successive weights of water. 

To determine the densities, known solutions were made up by adding 
known weights of water to known weights of salt in a large weighing bottle. 
Samples of these solutions were transferred to a pyecnometer which had been 
calibrated at 50.5°. Since the density of the solutions varies slightly from 
linearity, data are given for densities at a series of concentrations. 

The density of dry cesium formate was determined by weighing the de- 
hydrated salt in a pycnometer under hexane. From the weight of hexane 
and its density and the weight of salt, the density of the salt is readily com- 
puted. Assuming the density of hexane to be 0.6550, the density of the 
salt was found to be 2.99. The density of potassium formate was taken to 
be 1.91.5 

Viscosities were measured with a modified Ostwald viscometer. The 
viscometer was provided with a ground glass cap and solutions were made 
up in the viscometer. The instrument was calibrated with water and the 
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time of efflux was determined as a function of volume of water. The vis- 
cometir with its contents was weighed on a sensitive balance. Concen- 
tration could be changed by addition of water and further dilution could 
be made by withdrawing a known weight of solution and then making fur- 
ther addition of water. 

Densities of solutions in water at 50.5° are given for cesium formate and 
potassium formate in tables 1A and 1B, respectively. Molal concentra- 
tions appear in column 1, molar concentrations in column 2, and densities 
in column 3. 


TABLE 1 
DENSITY OF FORMATE SOLUTIONS AT 50.5°C. 


(A) CESIUM FORMATE— ————-—-(B) POTASSIUM FORMATE 
MOLAR MOLAR 
CONCENTRATION, C DENSITY CONCENTRATION, C DENSITY 


0.00 0.9878 0.00 0.9878 
1.193 1.142 3.967 1.160 
2.414 1.303 5.643 1.232 
4.892 1.613 9.452 1.392 
9.068 2.147 


TABLE 2 
CONDUCTANCE OF FORMATE SOLUTIONS AT 50.5°C. 


——=(A) CESIUM FORMATE ——~ -—-———(B) POTASSIUM FORMATE — 
EQUIVALENT EQUIVALENT 
MOLAR CONC COND, MOLAR CONC, COND, 


10.13 24.80 6. 503 63.03 
10.06 24.84 3.492 96.17 
9.034 18 10.45 30.82 
6.913 33 8.602 44.24 
4.629 78 15.52 9.780 
0.8104 1 13.66 14 
0.5430 5.8 . 5999 142. 
0.3055 5 . 2722 154 
3.231 98. 84 . 15438 161.6 
2.437 5 .404 110 
1.785 9 .325 126 
1.400 07471 168. 
.01973 176.6 
01088 179 


Conductance data for the cesium and potassium salts are presented in 
tables 2A and 2B, respectively. Molar concentrations are given in column 
1 and equivalent conductances in column 2. 

Viscosities for solutions of cesium formate and potassium formate are 
given in tables 3A and 3B, respectively. Molar concentrations are given 
in column 1, relative viscosities in column 2, and absolute viscosities in 
centipoises in column 3. 
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In tables 4A and 4B are presented values of the conductance-viscosity 
product for the cesium and potassium salts, respectively, at rounded molar 
concentrations. In column 2 are given values of the ratios of moles of 
salt per mole of water. Values of equivalent conductance and of viscosity 


TABLE 3 
VISCOSITY OF FORMATE SOLUTIONS AT 50.5°C. 
(A) CESIUM FORMATE (B) POTASSIUM FORMATE 
MOLAR CONC. n X 10? porses MOLAR CONC 7 X 10? porses 
9.605 2.612 13.95 5.173 
8.203 1.809 8.224 1.463 
6.834 1.338 6.090 028 
3.989 ). 8500 4.199 8277 
808 ). 7313 3.125 7296 
819 1.117 15.64 . 22% 
784 0.6533 12.26 3.403 
7496 0. 5847 10.38 216 
1.516 6245 


TABLE 4 
CoNDUCTANCE-VISCOSITY PRODUCT OF FoRMATES AT 50.5°C. 
MOLES H2O PER n X 102 
MOLAR CONC MOLE SALT CONDUCTANCE POISES An-PRODUCT 
(A) Cesium Formate 
00 192.9 5449 051 
25 oF 157.2 . 560 879 
00 01893 134 .610 .817 
00 ).04080 115.£ 670 774 
00 .09141 88 840 746 
00 1981 52 a7 72: 
00 .8018 32. 22 726 
00 3704 25.0 95 
(B) Potassium Formate 

00 ee 187.9 . 5449 
50 oxy 145.4 565 
00 01889 132.8 593 
00 06223 102.3 725 
00 1150 78.8 
00 1802 58 .¢ 
00 2634 41 .( 


905 .713 


3 165 683 
) 680 689 

00 5176 16.8 4.09 687 
14.00 6146 13.6 5.19 706 
15.50 9.9 89 781 


appear in columns 3 and 4, respectively, and values of the conductance- 
viscosity product in column 5. 

In figure 1 are shown plots of the equivalent conductance of cesium and 
potassium formates as a function of the square root of molar concentration. 
At low concentrations (not shown on figure), the conductance of the ce- 
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sium salt parallels and lies several A-units above that of the potassium salt. 
This is due to the fact that the conductance of the cesium ion is several 
A-units higher than that of the potassium ion. According to Noyes and 
Falk, the conductance of the cesium ion at 18°C. is 3.5 A-units greater than 
that of the potassium ion. These early measurements are not of high pre- 
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cision and the precise difference in the conductance of the cesium and po- 
tassium ions at 50.5° remains uncertain. However, it is safe to assume that 
the conductance of the cesium ion is the greater. We have extrapolated 
the conductance of potassium formate to zero concentration, using the 
A—<J/C plot. We obtained the value 187.9 at 50.5°. 





VoL. 39, 1953 CHEMISTRY: RICE AND KRAUS 807 


Using the data of Harned and Owen," for the conductance of the for- 
mate ion at 25° and those of Auerbach and Zeglin'! for the conductance of 
sodium formate at 18°, we obtain a vulue of 76 for the conductance of the 
formate ion on extrapolating linearly to 50.5°. From the equations of Har- 
ned and Owen we obtain a value of 112.1 for the conductance of the potas- 
sium ion at 50.5°. Thus, we have the value of 188.1 for the limiting con- 
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ductance of potassium formate at 50.5° in good agreement with our extra- 
polated value. 

At higher concentrations, the conductance curves of the two salts inter- 
sect at a concentration of approximately 2.5 molar. At still higher con- 
centrations the curves diverge increasingly. 

At the highest concentrations measured, the conductance of cesium and 
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potassium formates are 24.80 and 9.38, respectively. At a concentration 
of 19 molar, the conductances are 25.0 and 33.7, respectively. 

While the conductance of electrolytes at low concentrations is governed 
by the interaction between their ions, in concentrated solutions conduct- 
ance is mainly dependent on the viscosity of their solutions. In figure 2 
are shown values of the absolute viscosity of solutions of the two formates 
as a function of molar concentration. As may be seen from the figure, 
viscosity increases markedly with concentration in the case of both salts, 
being somewhat greater for the cesium salt. For a 10-M solution, the 
viscosity of the cesium salt is 5.4 times that of the pure solvent. For the 
potassium salt it is 3.9 times that of the solvent at 10 molar and 14.8 
times that of the solvent at 15.5 molar concentration. Up to a concentra- 
tion of C = 12, the measured viscosities in centipoises may be reproduced 
reasonably well by the equation : 


n = 0.545 + 0.0755C — 0.0135C? + 0.008C? 
for cesium formate, and by 
n = 0.545 + 0.0473C — 0.004C? + 0.00107C? 


for potassium formate. 

To a first order of approximation, the conductance of a completely dis- 
sociated electrolyte at different concentrations is proportional to the re- 
ciprocal of the viscosity of the solutions. Accordingly, we may compen- 
sate for the viscosity effect by multiplying the conductance by the viscos- 
ity. If the electrolyte is completely dissociated into its ions and the nature 
of the ions undergoes no change, the product should remain constant. 

In figure 3 are shown values of the conductance-viscosity product for the 
two formates as a function of concentration to near saturation. The prod- 
uct for the cesium salt is greater than that of the potassium salt over the 
entire concentration range. At low concentrations, the product for the 
cesium salt is greater than that of the potassium salt by an amount that 
corresponds to the greater conductance of the cesium ion. At lower con- 
centrations, when the viscosity effect is small, the conductance falls off 
sharply with increasing concentration in accordance with theory. As con- 
centration increases, the two curves diverge as a result of the greater vis- 
cosity of the solutions of the cesium salt. Ata concentration of 5 molar, the 
curves flatten out over a considerable concentration range. The product 
for cesium formate appears constant between 6 and 9 molar; for potas- 
sium formate it is constant between 7 and 13 molar. Over these ranges of 
concentration the number of moles of water per mole of salt decreases from 
6 to 3 for the cesium salt and from 5 to 2 for the potassium salt. 

It is apparent that at a concentration in the neighborhood of 4 to 5 
molar, where the solutions contain approximately 10 moles of water per 
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mole of salt, the character of the conductance process undergoes marked 
change. At these higher concentrations, the conductance-viscosity prod- 
uct changes but little while the number of water molecules becomes com- 
parable with the number of ions in solution. The hydration of the ions 
must be changing in this region. Apparently, the increase in mobility due 
to decreasing hydration just compensates for the increase in viscosity and 
the conductance-viscosity product remains constant. It will be noted 
that at the highest concentrations the curves for both salts have a marked 
upward trend. 

Campbell and Kartzmark have measured the conductance and viscos- 
ity of concentrated solutions of silver nitrate and ammonium nitrate in 
water at 95°. Values of the conductance-viscosity product for these 
salts along with our own for the formates are shown in figure 4, moles of 
salt per mole of water being plotted as abscissae. At lower concentrations, 
the product for the two nitrates falls off more sharply than it does for the 
formates. This is due to the lower values of dielectric constant and vis- 
cosity of water at the higher temperature, which results in larger values of 
the constants a and 6 of the Onsager equation. At higher concentrations 
the curves for the two nitrates differ from each other as they do also from 


those of the formates. 

For ammonium nitrate the product decreases continuously up to the 
highest concentration measured, two moles of water per mole of salt. The 
product for silver nitrate passes through a pronounced minimum at a con- 


centration which corresponds approximately to 12 moles of water per mole 
of salt. At concentrations above the minimum, the curve rises and ap- 
pears to flatten out but its precise form at the higher concentrations can- 
not be determined because of the scattering of the experimental values. 

In going from 22 to 9 moles of water per mole of salt, the conductance- 
viscosity product for silver nitrate decreases 0.8%; in going from 9 to 
6 moles of water per mole of salt the viscosity increases 277% while the con- 
ductance-viscosity product increases 9%. To the best of our knowledge, 
silver nitrate is completely dissociated into its ions. As the number of 
water molecules per mole of salt decreases, the hydration of the ions must 
decrease. However, over the concentration range from 22 to 9 moles of 
water per mole of salt, it is difficult to account for the larger effects that are 
observed on the basis of hydration alone. One cannot but wonder whether 
the postulates that underlie the theory of dilute electrolyte solutions re- 
main applicable at high concentrations. The product for ammonium 
nitrate at higher concentrations is much lower than that of the other salts. 
One is led to suspect that ammonium nitrate is not completely dissociated 
into its ions at high concentrations. There is a possibility of hydrogen 
bonding between the ammonium and the nitrate ions. 

From the density data, the volume change accompanying the process of 
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solution of the two formates at different concentrations has been computed. 
Values are shown graphically in figure 5. For cesium formate there is a 
volume contraction of 11.6 cc. at a concentration of 1 molar; it decreases 
gradually to a value of about 8.5 cc. at 5 molar and remains sensibly con- 
stant at higher concentrations. The volume contraction for potassium 
formate at | molar is 8.7 cc.; with increasing concentration the contraction 
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decreases sharply and levels off at a value of approximately 2.3 cc. at a con- 
centration of 5 molar. A 5 molar solution contains approximately 8 
moles of water per mole of salt. All properties of these solutions indicate 
a change in the nature of these solutions in the neighborhood of 10 or 8 
moles of water per mole of salt. It is difficult to escape the conclusion that 
the laws governing concentrated salt solutions differ fundamentally from 
those of dilute solutions. It is not to be expected that the properties of 
concentrated solutions may be accounted for by merely modifying some- 
what the laws that apply at low concentrations. 


* The investigation here reported was supported by ONR under Contract N7onr 
35809. 

+ This paper is based on a portion of a thesis presented by M. John Rice, Jr., in partial 
fulfillment of the requirements for the degree of Doctor of Philosophy in the Graduate 
School of Brown University, June, 1953. 
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DEVELOPMENT AND GENETICS OF A SEX-INFLUENCED 
TRAIT IN THE LIVERS OF MICE* 


By C. H. DANFORTH AND ELIZABETH CENTER 
DEPARTMENT OF ANATOMY, STANFORD UNIVERSITY SCHOOL OF MEDICINE 
Communicated June 8, 1953 


It is generally recognized that there are many traits which may be in- 
fluenced by the sex of the individual in which they occur. If one of the 
principal genes necessary for its appearance is located in the Y-chromosome 
a trait is designated as sex-limited, if in the X-chromosome as sex-linked, 
and if in addition a major autosomal gene is also necessary the trait is com- 
monly referred to as sex-influenced. However clear-cut these categories 
may appear in their verbal formulation, it is probably true that one could 
arrange a finely graded series from traits that are sex-limited in the strictest 
sense to those whose expression has little or no relation to sex. We wish at 
this time to report a condition in the mouse liver which might serve as one 
member of such a series, and to call attention to features in its development 
which explain why, on genetic grounds alone, this particular trait might 
sometimes appear to be almost completely sex-limited, and at other times 
not. 

Morphology of the Trait.—The liver of the mouse, like that of many re- 
lated rodents, is commonly described! as consisting of four primary divi- 
sions, a right lateral, a medial, a left lateral and a posterior, each of which, 
except for the left lateral, is in turn subdivided into two major lobes. In 
the present connection we are concerned with the right lateral division, 
whose smaller posterio-lateral, and larger anterio-medial components may 
conveniently be designated as lobes 1 and 2, respectively. The manifesta- 
tions in which we are especially interested involve only lobe 1, the dis- 
tinctive features of which may be mentioned briefly. 

In the females this lobe has the appearance of a somewhat distorted tri- 
angular pyramid whose base is penetrated by the inferior vena cava and 
whose apex is bent backward over the ventro-medial aspect of the right 
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kidney. The free part of the lobe presents three well-defined surfaces and 
three sharp margins, all of which converge to a pointed tip. One surface 
is in contact with the duodenum, one with the right kidney and suprarenal 
gland, and one with lobe 2 of the liver and with the peritoneum of the lateral 
body wall. In some males lobe | is similar to that of the female, in others 
it shows distinct deviations. 

When lobe | of the male differs from lobe 1 of the female it may do so in 
any one of several ways. It may show a slight notch in either the duodenal 
or the renal margin of its anterior surface; the two notches may both be 
present and connected by a shallow groove; the groove may be deeper 
causing the apex of the lobe to appear more like an appendage of the basal 
portion; or, finally, the basal part may itself become pointed and over- 
ride the morphological tip reducing it to the status of an inferior lobe. 
There are various minor variations and intergradations, but most speci- 
mens fit fairly well into one or another of the six categories set up in table 
1, the data for which were derived from some 3500 autopsies on mice rep- 
resenting both sexes and several different strains. 


TABLE 1 


DISTRIBUTION OF VARIATION 3 ACCORDING TO SEX 


LOBE 1 DUODENAL RENAL SHALLOW DEEP INFERIOR 
SEX INTACT NOTCH NOTCH GROOVE GROOVE LOBULE TOTALS 


Males 490 117 170 175 450 492 1894 
Females 1598 1 4 1 0 17 1621 


Table | serves to indicate the distribution of variations in a rather large 
aggregate, and especially to show differences in incidence between the two 
sexes. It may be mentioned here that ancestors of some of the mice used 
had been treated with nitrogen mustard or other possible mutagens, but 
where no differences were detected between such mice and their controls 
data from both groups have been included. Excluded from the table are 
records from earlier autopsies where full details were not recorded, a few 
doubtful cases, and instances where the configuration of the liver was sec- 
ondarily modified by partial or complete situs viscerum inversus, or by other 
anomalies. Summarizing the figures, it appears that among these animals 
deviations from an intact lobe 1 occurred in about 74 per cent of the males 
and in less than 1.5 per cent of the females. The true difference was pos- 
sibly slightly greater than this since a few of the variations recorded for 
females were of a minor or doubtful grade. 

When the same figures are broken down according to strain, differences 
among the males are shown to be considerable, those among the females 
relatively insignificant. Of five groups in which more than 100 individuals 
of each sex were examined, the incidence of intact lobes ranged from 97.5 to 
100 per cent for the females, and from less than one to over 71 per cent for 
the males. S-br mice gave the lowest and C3H the highest rates in both 
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sexes. Figures for the four genetically best defined of these groups are 
4 presented in table 2, from which the almost completely consistent data 
The C3H and C57-black animals en- 












from females have been omitted. 








TABLE 2 





DISTRIBUTION OF VARIATIONS ACCORDING TO STRAIN (MALES ONLY) 


LOBE | DUODENAL RENAL SHALLOW DEEP INFERIOR 
STRAIN INTACT NOTCH NOTCH GROOVE GROOVE LOBULE rOTALS 














C3H 82 1 11 2 13 6 115 
S-J 145 27 78 81 69 77 477 
C57-black 68 22 30 22 108 55 305 
S-br 2 11 9 11 125 154 312 






TABLE 3 





Direct GENETIC Tests (MALE YOUNG ONLY) 











puO SHAL- IN- 
MATINGS LOBE | DENAL RENAL LOW DEEP FERIOR 
PEMALES MALES INTACT NOTCH NOTCH GROOVE GROOVE LOBULE 













TOTALS 


















Swiss (4) Swiss 19 3 3 0 2 9 36 
2. S-br (8, 12) S-br (11) 0 0 0 0 5 8 13 
3. S-br Swiss 2 0 2 l 2 0 rf 
4. Swiss (1) S-br l 2 3 9 20 5 40 
5. F, and later 8 5 2 3 11 9 38 
6. S-J (7) S-]J $1 9 4 2 2 2 60 
7. S-J (6) S-br 0 1 1 0 2 3 7 
8. S-br (2,12) S-J 0 3 l 2 7 5 18 
9. F, and later 4 l 1 2 11 8 27 
10. C3H (11) C3H (12) 15 0 5 1 2 0 2: 
11. C3H (10) S-br (2) 0 0 0 1 8 10 19 

S-br (2,8) C3H (10) 0 5 3 4 18 1 31 









Figures in parentheses indicate other matings in which the same mice were involved, 





tered in the table were of the standard strains propagated under these des- 
ignations by the Roscoe B. Jackson Memorial Laboratory and were derived, 
respectively, from more than twenty and more than thirty successive gen- 
erations of brother-sister matings, followed by a further period of close 
inbreeding in which, however, only part of the matings were between actual 
sibs. The S-J animals came from a somewhat inbred but rather hetero- 
geneous ‘‘general-purpose’’ stock, and the S-br from a fairly homogeneous 
strain which has been maintained by the Department of Anatomy at 
Stanford for a number of years, its most distinctive characteristics being 
brown fur and a moderately short tail. 

Genetics of Variations in Lobe 1.—-Keferring to table 2, where only males 
are listed, it appears that when a distinction is made between intact lobes 
and those that are more or less subdivided, the former show frequencies of 
about 71, 30, 21, and 0.06 per cent in the C3H, S-J, C57-black, and S-br 
Corresponding figures for the fe- 
















males examined in the present series. 
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males are 100, 98, 99, and 97.5 per cent. Since all the mice were raised 
in the same laboratory and under similar conditions, these distributions 
indicate at least a partially genetic background for the observed variations. 
Breeding tests were carried out by making the matings first and deter- 
mining the characteristics of parents and offspring later. About 600 
young were obtained from appropriate matings of this sort. Since none 
of the females in this series showed any subdivision of lobe 1, only males 
are listed in table 3, which summarizes the results from direct genetic tests. 
Entries under mating | are of indifferent value for present purposes since 
the records failed in part to differentiate between offspring of males with 
and those without intact lobes. The items are included only for compari- 
son with those from matings 3 and 4. In all matings except | the exact 
condition of each parent is known. All females involved had an intact 
lobe 1, all S-br males a deep groove or inferior lobule, all S-J males an intact 
lobe 1, and the single male used in matings 10 and 12 had a renal notch. 
Matings 3 and 4 are the only ones in the series where any of the imme- 
diate offspring of the S-br parent of either sex showed an intact lobe 1. 
Inbred descendants of these matings showed the expected wide distribu- 
tion of the possible manifestations. Combining figures from the S-J and 
C3H lines, which seem to be the most similar both genetically and mor- 
phologically, we find that when inbred these two strains (matings 6 and 
10) produced a total of 83 male young of which 56 had an intact lobe 1, 
whereas when outcrossed to S-br (matings 7, 8, 11, 12), there were 75 


male young none of which had an intact lobe 1. Reciprocal crosses 
showed no significant differences in distribution of the traits, which indi- 
‘ates that the tendency to a subdivision of lobe | is transmitted about 
equally by both sexes, even though its expression is limited almost exclu- 


sively to males. 

As to the various grades of lobulation, it may be observed that a vertical 
line drawn through table 3 between the columns headed ‘‘renal notch’’ and 
“shallow groove’’ would divide each group asymmetrically, with the larger 
fraction either to the right or left of the line. This might suggest a natural 
separation at this point between two alternative but somewhat overlapping 
traits. In this connection it may be noted that the C3H male used in 
mating 10 had a renal notch but behaved in the tests like S-J or other C3H 
males with intact lobes. The data, however, are not sufficient to differ- 
entiate on genetic grounds between the possibilities that there are two 
alternative traits with incomplete dominance and variable expressivity, 
or that the several manifestations are individually dependent on varying 
combinations of multiple or cumulative genes. 

Embryology of Variations in Lobe 1.—The ontogeny of lobe | throws an 
entirely new light on the problem and provides a complete change of focus. 
It appears from a study of embryos in progressive stages of development 





VoL. 39, 1953 GENETICS: DANFORTH AND CENTER 815 


that the variations in lobe | which have been described are the direct, and 
not indirect, morphogenic effects of the gonads themselves. In preparing 
the smaller embryos for study our procedure has been to remove each one 
from its membranes, immerse it in water for about an hour and then in 
Bouin's solution for from 6 to 20 hours, depending on its size, after which it 
was washed in weak alcohol and dissected under the low power of a bin- 


ocular microscope. In these dissections the routine first step was to re- 


move the skin and subjacent tissues from the right side of the back down to 
the level of the kidney, duodenum, and liver. 

In embryos at about 12 to 13 days after fertilization, and of approxi- 
mately 10 mm. in length,” * lobe | is the most rudimentary of the seven 
major portions of the liver. At this stage it is very small and compressed 
on three sides by the much larger lobe 2, the right kidney and the duo- 
denum, while the relatively undifferentiated right gonad and its accompany- 
ing ducts cover it and project into its dorsal surface. This gives the lobe 
as viewed from above somewhat the appearance of a boat or dugout, with 
a central depression and two raised margins embracing the ventral portion 
of the gonad. Development of embryos from this point on is rapid,‘~* 
and this is especially true of the region under consideration. At about this 
time the ovary, while still small and cylindrical, begins to move obliquely 
across the dorso-caudal pole of the kidney, early freeing lobe | in the female 
to undergo a considerable dorsal expansion, in the course of which the 
groove that had been occupied by the gonad generally becomes completely 
smoothed out leaving no trace of its former existence. The male gonad, 
by contrast enlarges considerably before beginning its descent, and this 
has the effect of still further deepening and widening the groove as well as 
somewhat retarding the development of lobe 1 as a whole. When the 
testis finally does descend, it is not immediately lifted off lobe 1, as is the 
ovary, but continues to impinge on its dorsal surface until it has completely 
passed over the duodenum and on into the pelvis. By the time the testis 
has lost contact with lobe 1, the two sides of the latter have often developed 
appreciably, accentuating the depression between them and causing it to 
remain as a permanent notch or groove. 

Differential early development on the two sides of the groove apparently 
determines the precise configuration which the lobe will ultimately assume, 
and all forms of adult lobulation may be found taking shape in embryos 
of around 15 to 16 days. That in some adult male livers lobe 1 shows these 
traits and in others does not is probably due to the length of time during 
which the mechanical impact of the testis persists. That this period 
varies is demonstrated by the fact that sometimes in members of the same 
litter one right testis will be found deep in the groove, another just leaving 
it, and a third in the pelvis well beyond the liver. It is, of course, impos- 
sible to follow individual development in preserved specimens, but exami- 





816 GENETICS: DANFORTH AND CENTER Proc. N. A. S. 


nation of a considerable mass of embryonic material makes it seem very 
probable that those rudiments which are early freed from the impeding 
influence of the testis are the ones which develop as intact lobes, while 
those in which the testis remains until the liver has bulged out around it 
are the ones which in the adult show varying degrees of lobulation. It 
would thus appear that the immediate factor determining whether or not 
lobe 1 will be intact is a temporal one, involving the time and speed at 
which the testis loses contact with the liver. 

Discussion.—-It has seemed worth while to record these rather trivial 
variations in the mouse liver because of the bearing they may have on 
questions of wider interest. In the first place, it is instructive to find that 
a trait which in one strain might appear to be almost completely sex- 
limited can be shown in another strain not to be so in any strict genetic 
sense. In this instance, the genes immediately involved are clearly auto- 
somal and probably effective through their regulation of developmental 
rate at a crucial moment in embryonic life. For their effect to be made ap- 
parent at all, however, there must be the added impediment of a gonad, 
the presence of which in its turn is determined by sex chromosomes. Ex- 
istence of the autosomal part of the complex would probably never be sus- 
pected from the examination of female material even though its manifes- 
tations are clearly recognizable in the male. 

It is probably significant that more than thirty generations of intensive 
inbreeding in the C57-black line has failed to fix a trait that is obviously 
hereditary (albeit with varying penetrance and expressivity), while less 
intensive selection has come much closer to doing so in the S-br line. The 
latter has been maintained for some ten years with no special breeding pro- 
gram other than the periodic reduction of the stock to three or four selected 
individuals. One might wonder if some strains of mice which persist with 
fair vitality despite intensive inbreeding are in reality being continually 
selected for vigor which fortuitously happens to be associated with hetero- 
zygosity in certain unsuspected genes, as professor P. W. Gregory’ believes 
to be the case in some herds of cattle. 

On the embryological side, if the inferences drawn here are substantiated 
by such experimental tests as it would seem feasible to make, the develop- 
ment of the lobe | in the mouse liver will afford one more illustration of the 
importance of ‘developmental rate in structural expression” often empha- 
sized by Professor Stockard® and later embryologists. 

Summary.—Certain lobulations or incisions in the more caudal lobe 
(lobe 1) of the right lateral division of the mouse liver are found to occur 
with a rather high frequency in males, only rarely in females. A genetic 
and embryological study of these manifestations reveals that they are the 
result of impressions made by the gonads in early embryonic life, and that 
differences in grade of expression are apparently mediated through auto- 
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somal genes rather than those in the sex chromosomes. The difference in 
incidence between the sexes is due to the fact that the ovary is early re- 
moved from conta:t with the liver, while the testis commonly maintains 
such contact with resulting pressure on lobe | during the period of its most 
rapid development. The evidence indicates that grades of expression, 
which in the male range from slight to marked deviations from the female 
type, are conditioned by the time at which the gonad in its descent loses 
contact with the liver, the final expression being in effect an index of the 
developmental rate at a period around the fifteenth or sixteenth day of 
intra-uterine life. 

* The observations reported in this paper are incidental to an investigation of possible 
mutagenic effects of nitrogen mustards, supported by a research grant from the National 
Institutes of Health, Public Health Service. 
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MUTATION RATES AND LETHAL GENE FREQUENCIES 
IN POPULATIONS OF DROSOPHILA MELANOGASTER 


By R. P. Levine* AnD P. T. Ives 
DEPARTMENT OF BroLoGy, AMHERST COLLEGE, AMHERST, MASSACHUSETTS 
Communicated by Karl Sax, June 8, 1953 


It is now generally accepted that gene mutation is the basis of evolu- 
tionary genetic change. Though many data are available concerning the 
relative frequencies of mutations in populations of several species of 
Drosophila,' there is a paucity of information regarding spontaneous mu- 
tation rates in natural populations of these species. The present article is 
a report of preliminary investigations of sex-linked lethal mutation rates 
in natural and artificial cage populations of D. melanogaster and the fre- 
quency of second chromosome lethals in the artificial populations. 

The Experimental Procedure.—The samples of natural populations ana- 
lyzed were obtained from Blacksburg, Va.; Canonsburg, Pa.; Austin, 
Tex.; Chapel Hill, N. C.; and Amherst, Mass. 

All samples, with the exception of Chapel Hill, were handled in the fol- 
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lowing manner. As soon as the flies were obtained from the field, from 
40 to 100 females were isolated singly in 25 K 95-mm. vials of a standard 
culture medium and placed at 25°C. One F, male, no more than three 
days old, was taken from each of the cultures. He was mated for two 
days to two females carrying the X-chromosome complex known as Muller- 
5. The females were allowed to oviposit for a week. The F, were shaken 
in each case to a new culture where they remained for two days to insure 
maximum insemmination. Then they were anesthetized and put into in- 
dividual cultures, one pair per vial, for the standard Muller-5 sex-linked 
mutation rate test. Each family of test cultures from a single F; male was 
numbered. There were about fifty tests in most of the families, meaning 
tests of 50 sperm for the F; male. The numbering of each family is essen- 
tial since the occurrence of more than one mutation within a family might 
result from a single mutation in the F, male early in spermatogenesis, thus 
affecting more than one sperm. When two or more mutations occurred 
within a family, they were localized by cross-over tests with the X-chro- 
mosome markers, yellow, cut, raspberry, forked (y ct ras f) in order to ascer- 
tain how many independent mutations they represented. 

The test generation was scored for lethal mutations by means of a bin- 
ocular microscope. The presence of three or more normal males together in 
the field was considered suff cient to score the test as non-mutant. The 
field showed a maximum of about thirty flies. All other cases were anes- 
thetized for scoring. In cases of less than 10 per cent normal males, retest 
matings were made, and the final scoring was based on a minimum count 
of 100 flies in two generations. Tests showing no normal males were scored 
as lethals. Those showing up to 10 per cent normal males were scored as 
semilethals. Obvious phenotypic mutations were also scored. The mu- 
tation rate was based on the total of the three kinds of mutations which 
occurred independently in each family. All lethals from the same fam- 
ily which appeared to be at the same locus in cross-over tests were consid- 
ered to be the result of a single mutation. Such cases were rare in the 
study. 

Only a small number of females were available from Chapel Hill. They 
were allowed to lay eggs in half-pint bottle cultures, from which many 
young F; males were mated individually to Muller-5 females as above. 

The sex-linked lethal mutation rate is being followed in two cage popula- 
tions of D. melanogaster. The cages used and the technique of their main- 
tenance are described by Wallace.? The first population, known as , 
was derived from a collection of flies from South Amherst, Mass., in 
September of 1952. Fifty fertilized females, sisters of the males tested in 
the Amherst sample, were used to start the cage on October 3, 1952. 
The cage was kept at 25°C. Egg samples were removed from the cage at 
approximately monthly intervals from the time of the emergence of the 





VoL, 39, 1953 GENETICS: LEVINE AND IVES 819 


first generation in the cage. Resulting young males were mated singly to 
Muller-5 females and tested as above. 

The second cage population, known as O, was derived from an approxi- 
mately homozygous Oregon-R strain which had been inbred for 95 gen- 
erations by brother-sister matings; the last 40 generations in this labora- 
tory. The 50 females which were used to start the cage and their brothers 
which were used for an initial determination of mutation rate were obtained 
from a half-pint laboratory stock culture. The cage population began on 
November 17, 1952, and was treated in the same manner as the 7 popu 
lation. The two cages were kept in the same incubator. 

In order to determine the second chromosome lethal gene frequencies in 
the M and O populations, adult males were removed at approximately 


TABLE 1 
MuTATION RATES IN NATURAL POPULATIONS 
NO. OF 
COLLECTION NO. OF NO, OF SEMI MUTATION 
LOCALITIES TIME TESTS LETHALS LETHALS RATE 
Blacksburg ( Va.) 

1 May ’52 3768 x 
2 Oct. ’52 2550 x 
Canonsburg ( Pa.) July ’52 2571 x 
Austin (Tex. ) Nov. ’52 3323 j x 
Chapel Hill (N. C.) Nov. ’52 3456 x 
Amherst ( Mass.) Sept. ’52 3191 x 


TABLE 2 
MUTATION RATES IN THE M (AMHERST) CAGE 
TIME 
(pays NO. OF 


PROM START NO. OF NO. OF SEMI- 
SAMPLE OF CAGE) TESTS LETHALS LETHALS MUTATION RATE 


1 0 3191 3 2 1.56 X 1078 
2 33 3326 7 3 3.01 X 107% 
3 67 3829 12 1 3.40 XK 1073 
4 98 3310 9 0 2.42 X 1078 


monthly intervals and the lethal frequency was determined in the manner 
described by Ives.* Lethals and semilethals (those yielding less than 17 
per cent normal flies) have been grouped together in the tables. 

Mutation Rates——The mutation rates determined for the five natural 
populations studied are given in table 1. The variability in the rates from 
population to population is not significant. (P is 0.47 as determined by 
chi-square.) The populations studied represent collections of flies from 
geographically diverse habitats, but it is by no means certain that the larval 
habitat of the populations differed significantly in an ecological sense dur- 
ing the time of greatest reproductive activity. 

Recorded in table 2 are the mutation rates measured for the / popula- 
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tion. Sample i is the mutation rate at the time the cage was established. 
The remaining values were obtained from egg samples taken from the cage 
at the times indicated in the table. A P value of 0.70 for samples 2, 3, and 
4 indicates only random variation. However, upon comparing the pooled 
data of the M population with the pooled data of the natural populations 
one obtains a P value of 0.016 indicating a significant difference between 
the M cage population and the natural populations with respect to muta- 
tion rate. 

Table 3 shows the mutation rates of the O population. Sample | is 


TABLE 3 
MUTATION RATES IN THE O (OREGON-R) CAGE 
TIME 

(DAYS NO, OF 

FROM START NO. OF NO. OF SEMI 
SAMPLE OF CAGE) THKSTS LETHALS LETHALS MUTATION RATE 
| 0 3278 2 2 22 % 1075 
2 33 3412 0 47 XK 1073 
63 2835 : 0 a0 *% 3075 
xX 10-3 


3 
4 97 3112 0 2.25 


TABLE 4 
SECOND CHROMOSOME LETHAL GENE FREQUENCIES 
TIME 
(DAYS FROM START NO. OF PER CENT 
OF CAGE) TESTS LETHALS OF LETHALS 
M Population 
190 y 37 .§ 
105 2: 22 
200 30 
133 26 
125 28 


55 34. 


O Population 
231 
205 
161 
162 


the mutation rate of the inbred Oregon-R stock at the time that the cage 
population was established. A P value of 0.85 for samples 2, 3, and 4 
suggests only random variation. Since the past history and the artificial 
condition of the O population is quite different from that of the natural 
populations, there is nothing instructive in a comparison of the mutation 
rates of the O population and the five natural populations. However, the 
mutation rates of samples 2, 3, and 4 from both the / and O populations 
make a valuable comparison between flies with different past histories 
subjected to similar cage environments. Such a comparison by a chi- 
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square test gives a P value of 0.26. The differences between the two popu- 
lations are not significant, due to the almost linear increase in the muta- 
tion rate in the O population to a point very close to the mutation rate of 
the / population in the last sample. 

Lethal Gene Frequencies.—-The second chromosome lethal gene frequen- 
cies determined for both the / and O populations are given in table 4. 
It will be noted that the / population had an initial chromosome lethal 
gene frequency of 37.9 per cent. This represents the lethal frequency in 
the sample of the Amherst natural population from which F;, flies were 
used to establish the cage population. Subsequent samples showed a 
decrease in frequency in the first generation in the cage (sample 2) and 
what appear to be random fluctuations thereafter. 

The O population, which at its inception had a lethal gene frequency of 
1.3 per cent, has shown a steady increase in frequency in subsequent 
samplings. 

Discussion.-It has been claimed by Berg* ° that sex-linked lethal mu- 
tation rates in different natural populations of D. melanogaster vary from 
population to population according to ecological differences. She offers 
the hypothesis that the mutation rate is subject to natural selection and that 
rates are lower in populations existing under specialized ecological situa- 
tions in which the range of adaptedness of the individual is narrow. On 
the other hand, populations found in more generalized environments and 
thus under conditions in which there would be wide range of adaptedness 
would be expected to have a higher mutation rate. The action of natural 
selection would not be expected to be as stringent here when compared to 
the situation existing in the specialized environments. Though her data 
are not suggestive, they are based on a relatively small number of tests and 
experimental procedures are often obscure. However, the results of Dob- 
zhansky, Spassky, and Spassky*® obtained from an analysis of lethal mu- 
tation rates in D. willistoni and D. prosaltans are of interest in 
this respect. The former species is common throughout most of tropical 
America while the latter is a rare species, restricted to a small number of 
localities. It has been shown that the mutation rate is significantly higher 
in the widely distributed D. willistont when compared with D. prosaltans 
which appears to require specialized environments for its maintenance. 
Thus, their results will tend to support the Berg hypothesis. They point 
out further that, ‘if mutation rates are adaptively modified by ecological 
conditions, the divergence among species may be expected to be greater 
than amony populations of a single species.’’ The results presented here 
for a number of possibly ecologically diverse populations of D. melanogaster 
bear this out. However, the mutation rates in artificial populations pre- 
sent a somewhat different picture. Here, as was the case with the V/ 
population, the mutation rate increased significantly above the character- 
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istic for the natural populations. In addition, the O population appears 
to exhibit a steady increase in mutation rate above that found to be 
characteristic of the inbred Oregon-R stock. 

These results suggest that mutation rate can be modified in populations 
permitted to develop under completely new environmental conditions. 
The representatives of the Amherst population in the population cage are 
subjected to an entirely new spectrum of selective agents; competition is 
keen for food as well as for egg laying and pupation sites. It is possible, 
therefore, to conceive of the increased mutation rate as being the conse- 
quence of a population’s response to natural selection under the new con- 
ditions afforded by the cage. This would fit at least in part, the hypoth- 
esis of Berg. The drop in the second chromosome lethal gene frequency 
the initial 37.9 to 22.9 per cent after the first cage generation had devel- 
oped was not, however, due to cage conditions. It was due, in all proba- 
bility, to elimination of lethal chromosomes caused by the vigorous brother- 
sister mating of flies which were put into the cage at the start. 

The O population showed a steady increase in mutation rate. The 
ancestors of the cage population have been subjected to the rather special- 
ized conditions of laboratory culture for 95 generations. Genetic varia- 
bility was kept at a minimum level by inbreeding and any factors tending 
to increase mutation rate would also be kept at a minimum mutation rate 
level by this process. However, when permitted to outcross in the popu- 
lation cage, the lethal gene frequency begins to rise and undoubtedly other 
genetic variability also increased. Similar increases in lethal gene fre- 
quency have been obtained for a number of cage populations by Wallace.’ 
Accompanying the increase in lethal frequency is an apparent gradual rise 
in mutation rate. This increase can be explained either on the basis of a 
relaxation of selection following the Berg hypothesis or by the action of 
mutator gene complexes developing under cage conditions of increased pop- 
ulation size. Thus, the increase in mutation rate is perhaps a reflection 
of the greater diversity of the cage environment as compared with the con- 
ditions obtaining in single vials of brother-sister single pair matings. 

The findings of no significant differences between mutation rates of sam- 
ples of flies taken from different populations of D. melanogaster is not in ex- 
act accord with the finding of Demerec.* However, most of his samples 
showed a low mutation rate and only a few were high. In addition, his 
samples were derived from stocks which had been kept under laboratory 
conditions for several generations. It may be that tests of additional 
areas, using our technique, will show some to be significantly different. 
None of the wild samples came from areas of really small breeding units as 
indicated by their frequency of second chromosome lethals (data to be pub- 
lished by Ives). 

It is possible that the apparent difference in mutation between D. 
prosaltans and D. willistoni is due as much to differences in mutator alleles 
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in populations of different sizes as it is to ecological differences between 
their habitats. Having a much larger population area and a larger breeding 
population size, D. willistont would be expected to have a higher concen- 
tration of mutator alleles similar to the mutator high described by Ives? 
and hence a higher mutation rate than D. prosaltans. This is what was ob- 
served to be the case. However, the difference in mutator complexes and 
in population size would both depend in turn on the different ecological 
requirements of the two species. They may, therefore, be mechanisms by 
which ecological differences can express themselves in the mutation rate. 

Further light on the nature of the changes in the frequency of lethals in 
the 7 cage should come from allelism tests of the lethals recovered in the 
various samples. Since only 200 second chromosomes, at the most, con- 
stituted the gene pool, there should be a higher amount of allelism in sam- 
ples of lethals taken from the population cage than in the initial sample 
from the much larger natural population. Preliminary results indicate 
that this expectation will be realized. 

It is obvious that much more extensive work needs to be done to eluci- 
date the problems arising in this preliminary study of the dynamics of 
natural and artificial populations of D. melanogaster. In summary, it 
may be stated that preliminary studies of mutation rates in flies derived 
from five geographically diverse natural populations of D. melanogaster 
showed a comparative uniformity. An artificial population derived from 
flies collected in one area and maintained in a cage showed an increase in 
mutation rate over that found in natural populations. A cage population 
developed from a stock with a long history of close inbreeding also showed 
a tendency for increased mutation rate. It is suggested that these mu- 
tation rate increases reflect the effects of the cage environment. 
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ANTI-f AND THE “NEW” Rh ANTIGEN IT DEFINES 


By RutuH SANGER, R. R. RAck, R. E. ROSENFIELD, P. VOGEL, AND NATALIE 
GIBBEL 
MEDICAL RESEARCH COUNCIL, BLOOD Group RESEARCH UNIT, THE LISTER INSTITUTE, 
LONDON, ENGLAND; THE Mount Sinat HospitraL, NEw YorK, NEW YorRK; AND 
BUREAU OF LABORATORIES, DEPARTMENT OF HEALTH, NEW YORK, NEW YORK 


Communicated by L. C. Dunn, June 19, 1953 


A “new” Rh antibody, anti-f, has been found in the serum of a man, 
Mr. Cl., aged 30 who suffers from hemophilia. A brief summary of the 
investigation of this antibody has recently been published.! 

Since childhood Mr. Cl. has received about 35 transfusions of blood, 
the last of which, in December, 1951, was followed by a rigor. He was 
referred to us, for the first time, in March, 1952, by the Hemophilia Founda- 
tion, Inc. At that time we found three blood-group antibodies in his 
serum: anti-B, anti-A, and anti-S. In addition there was at least one 
other unidentified antibody. In the hope of clarifying the reactions of 
this antibody, 0.25 ml. of O, NsNs, cdE/cde, kk blood was injected on two 
occasions: following this the presence of anti-N was disclosed and, in 
addition, the ‘“‘new”’ antibody, anti-/. 

The blood groups of Mr. Cl. are: A,, MsMs, P+, CDe/cDE(R,R:), 
Lu(a—), kk, Le(a + 6 —), Fy(a +), Jk(a + 6 +). His serum can be 
absorbed free of anti-B, anti-N, anti-S, and anti-A leaving anti-f alone; 
it is the anti-f left in this absorbed serum that is the subject of this paper. 
A detailed account of the Rh blood groups prior to the recognition of anti-f 
may be found in the text of Race and Sanger.’ 


THE REACTIONS OF ANTI-f. 


Methods.—The anti-f in the serum of Mr. Cl. was found to react most 
clearly with an albumin technique: a volume (0.01 ml.) of a 4% suspension 
of red cells in 20% bovine albumin was added to an equal volume of the 
serum. The mixtures were left at 20°C. for from one to two hours, then 
centrifuged at about 1,000 r. p. m. for 90 seconds. The cell deposit was 
then pipetted on to a slide; if agglutination was not visible to the naked 
eye it was looked for with the low power of the microscope. 

More than half of the red cell samples were tested against titrations of 
anti-f, the diluent being 20°, albumin; the rest were tested against the 
undiluted serum. The dilutions were made in bulk and then distributed 
to the appropriate tubes. 

Scores were given to the reactions: agglutination visible to the naked 


eye scored 10, lesser positive reactions scored 8, 5, 3, or 2. It soon became 
clear that there was a marked dosage difference, for cde/cde, cDe/cde and 
cD“e/cde score more highly than do other positive samples: 
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65 cde/ cde, cDe/cde and cD“e/cde samples; mean score = 32.5 
99 other positive samples; mean score = 16.0 
An indication of the variability of the method is given by the following 
scores for a cde/cde person (J.E.N.) tested on 12 occasions spread over 5 
months. 


35, 44, 30, 30, 33, 36, 38, 35, 35, 38, 33, 34; mean = 35.1 + 3.8 
and for a CDe/cde person (M.E.M.) tested on 11 occasions: 
6, 5, 18, 12, 10, 14,:16,. 32; 7,10, 3; mean = 9.9 + 4.0 


Results: The Usual Reactions..-Table 1 shows the results of testing 
with anti-f samples of blood from 328 unrelated white persons and from 
125 of their relatives. From these results it is clear that the antigen 
recognized by anti-f is part of the Rh complex of antigens; this is also 
shown by the linkage tests described below. Contingency tables show 
that there is no evidence of an association between the f antigen and any 
other of the known blood-group antigens, sex, or the ability to taste phenyl 


thiocarbamide. 

Many of the ‘‘unrelated’’ samples in table | were selected, as can be seen 
from the rarity of their Rh groups. It can be calculated that anti-f will 
agglutinate approximately 66° of samples of blood from unselected white 


people. 

The reactions of anti-f have been confirmed by absorption tests with 
red cells representing most of the critical genotypes. A volume of serum, 
0.3 ml., was absorbed with an equal volume of packed cells for 15 minutes, 
at 20°C.; four successive absorptions were done. The following red 
cells failed to absorb anti-f: CDe/CDe (6 samples), cDE/cDE, CDe/cDE, 
Cde/cdE (2 samples). The following red cells did absorb anti-f: cde/cde, 
cDe/cde, CDe/cde (2 samples), cDE/cde, cdi /cde, Cde/cde, CDE/cDe, 
and Cdki/cde; as expected cde/cde and cDe/cde were the most efficient 
bloods for absorption. Anti-f could be eluted from the various absorbing 
cells and the eluates when tested behaved alike, thus giving direct evidence 
that the positive reactions of the anti-f serum with red cells of different 
genotypes are due to one antibody. 

Table 1 becomes clear when it is realized that anti-f is reacting with a 
product of the chromosomes cde, cDe, and cD“e. These chromosomes 
may consequently be distinguished by the addition of f. If the absence of 
f is indicated by F the chromosomes are: 


cdef, <Def, cD“ef, 
CDeF, cDEF, C“DeF, CdeF, cdEF, CD“eF, cD“EF, 
CDEF, CD“EF, C*deF, CdEF, cDE“F, C“DeF, and C“dEF. 
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This shows the usual arrangement of f ir. Caucasian chromosomes, and 
the evidence is given in table 2. 

Less Common Reactions.—For reasons given below we expect certain 
less common combinations to exist—namely cDeF, cdeF, CDef and cDEf. 


TABLE 1 
THE REACTIONS OF ANTI-f WITH THE RED CELLS oF 328 UNRELATED, BUT OFTEN 
SELECTED, WHITE PEOPLE AND 125 oF THEIR RELATIVES 


POSITIVE . NEGATIVE—-— ‘ 
GENOTYPE UNRELATED RELATIVES UNRELATED RELATIVES 


CDe/cde 84 26 
cde/cde 61 28 
cDE/cde 30 18 
cDe/cde 19 4 
cd E/cde 9 0 
Cde/cde 2 
C“De/cde 0 
CD“e/cde 

CDe/cDe 

cDE/cDe 

CdE/cde 

Cvde/cde 

CDE/cDe 

CDE/cde 

CD“E/cde 1 

cDe/cde 3 

C“De/cde 1 

CdE/cde 1 2 
CDe/CDe 0 0 
CDe/cDE 0 0 
cDE/cDE 0 0 
C¥De/CDe 0 0 
CY’De/cDE () 0 
CPe/Cde 0 0 
cDE/cdE 0 0 
CDe/cdE 0 9) 
cDE/Cde 0 0 
CDe/CDE 0 0 
Cde/cdE 0 0 
CDe/CdE 0 6 
CDe/CD*e 0 0 
cD“e/cdE 0 0 
Cde/Cde 0 0 
cDE“/cde 0 0 
—D-—/-—-D- 0 0 


tb 


— bdo oO 
uo 


~ 


—— bo bt ort 


= © 


We have found an example of cdeF/cdef. The red cells of a member of the 
Staff of the Lister Institute, although cde/cde, persisted in giving a single 
dose score with anti-f: his scores were 9, 13, 8, 10, 12, and 7 as against 
an average for all other cde/cde of 32.5 (the next lowest score for any 
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unrelated cde/cde was 21). The absence of f from one of his cde chromo- 
comes is confirmed by family studies. His sister is also cdeF/cdef for she 
scores only 14. Their father is cDE“F/cdeF since his red cells are not 
agglutinated by anti-f nor will they absorb it (family 19, table 3). 


TABLE 2 
EVIDENCE FOR THE EXISTENCE OF THE FOLLOWING CHROMOSOMES 


SAMPLES POSITIVE WITH ANTI-f 
SINGLE DOSE DOUBLE DOSE 


CDe/cde, cDE/cde 
Cde/cde, cdE/cde 
C“De/cde, CdE/cde 
CD“e/cde, CD“E/cde 
C¥de/cde, CDE/cde 
CDe/cDe, cDE/cDe 
CDE/cDe 


CHROMO- 
SOMES 


cdef 


BLOOD SAMPLES NEGATIVE WITH 
ANTI-/ 


cde/cde 
cDe ‘cde 
cD“e/cde 


cDe/cde 


cDe/cde 
CDe/CDe, CDe/cDE, C“De/CDe 
CDe/CdE, CDe/CDE, CDe/Cde 
CDe/cdE, CDe/CD*e 
cDE/cDE, CDe/cDE, C¥De/cDE 
cDE/Cde, cDE/cdE 
CDe/cdE, cDE/cdE, Cde/cdE 
cD*E/cdE 


CDe/cde, CDe/cDe 


cDEF cDE/cde, cDE/cDe 


cdEF cd E/cde 


CdeF 


C’DeF 
CDEF 
CdEF 


CDe/Cde, cDE/Cde, Cde/cdE 
Cde/Cde 

C’De/CDe, C“’De/cDE 
CDe/CDE 

CDe/CdE 


Cde/cde 


C’De/cde 
CDE/cDe, CDE/cde 
CdE/cde 


C“dEF 
CvdeF 
CD“eF 
cD“EF 
CD“EF 
C“DeF 
cDE“F 
—D-—-—-t 


C“dE/cde 
Cvde/cde 
CDe/CD*e CD“e/cde 
cD“/cdE 
CD“E/cde 
C“De/cde 
cDE“/cde* 
—D-—/-D- 
* The genotype is cDE“F/cdeF, see text. 


t The chromosome — D— does not carry f. Until anti-F is found there must be some 


doubt whether it is —D-—— or —D—F. 


We thought we had identified the chromosome cDeF when a sample of 


blood apparently cDe/cde persisted in giving a single dose score with anti-f: 
ultimately we found that this sample was in fact C“De/cde for the cells 
were agglutinated by some anti-C sera but not by others. 

Family investigations.—Table 3 shows the results of testing 47 families. 
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TABLE 3 
THe RESULTS OF TEesTING 47 FAMILIES, OR PARTS THEREOF, WITH ANTI-f 


FAMILY 
NO FATHER MOTHER 


CHILDREN 





1 - —_ 
CDe/CDe CDe/cDE CDe/CDe CDe/CDe 


CDe/CDe CDe/CDe | CDe/CDe CDe/CDe 


CDe/CDe CDe/CDe | CDe/CDe CDe/CDe CDe/CDe CDe/CDe 


cDE/Cde CDE/cdE Cde/cdE 


CDe/cdE D“E/cde I cD“E/cdE cdE/cde 


CDe/cDE Cde/cde | CDe/Cde 

+ - | + + 
cDE/cde CDe/Cde Cde/cde CDe/cDE | Cde/cde 

or 
Cde/cDE 

< 4 | . 

CDe/cDE cdE/cde cDE/cdE cDE/cdE 
+ i + + + - 


CDe/CD“e| cDE/cde | CDe/cde CD"e/cde CDe/cde CDe/cDE 
_ i + + + | = 
CDe/cDE | cDE/cde CDe/cde CDe/cde | CDe/cDE 


+ -_ iH 
cDE/cde CDe/cDE | CDe/cDE cDE/cDE | 
+ _ | + + | _ ~_ 
CDe/cde CDe/cDE | cDE/cde cDE/cde CDe/cDE | CDe/CDe 





4 | 4 | 7 
CDe/CDe cdE/cde | CDe/cde | CDe/cdE 
+ ase | rs - 
CDe/cde | CDe/CDe | CDe/cde CDe/CDe 
+ + + 
cDE/cDE | cde/cde cDE/cde | cDE/cde 


CDe/CDE | cde/cde | CDE/cde CDE/cde 


Cde/cdE cde/cde | Cde/cde | Cde/cde 
+ | + 
cDE/cDe | cDE/CDE| CDE/cDe 
4 | + + 
cDE“/cde CDe/cde cde/cde cde/cde 
4 + i ~ + 
CDe/cde CDe/cde i} CDe/CDe cde/cde 
t oe + . 
CDe/cde CDe/cde CDe/cde CDe/CDe 
+ + - 
cDE/cde CDe/cde | CDe/cDE 
+ + { _ | _ 
CDe/cde CDe/cde | CDe/CDe CDe/CDe 
+ | + - _ - 
CDe/cde | cDe/cde CDe/CDe CDe/CDe | CDe/CDe CDe/CDe 
or | | | 
CDe/cDe 


Notes.—-The fathers in families 1, 28, and 29 are brothers. The mothers in families 
10 and 11 are sisters and so are the mothers in families 15 and 24. The father in family 
22 is the first child in family 10, and the father in family 45 is the child in family 18. 
The mother in family 47 is the third child in family 46. 





VoL. 39, 1953 GENETICS: SANGER, ET AL 


TABLE 3 (Continued) 
FAMILY 
NO. | FATHER MOTHER 


CHILDREN 


cDE/cde cdE/cde | ceDE/cdE eDE/cdE 


+ | _ + + 
CDe/cde cDE/cde | CDe/cDE cde/cde cde/cde 
+ + + 
cDE/cde CD“E/cde CD“E/cde cDE/cde cde/cde CD“E/cde 
+ + | + + 
CDe/cde cde/cde i cde/cde cde/cde 
+ + + + + 
cDE/cde cde/cde cDE/cde cDE/cde «DE/cde 
+ + + + 
CDe/cde cde/cde | cde/cde CDe/cde 
+ + | + + + 
cDE/cde | cDE/cde | cDE/cde cDE/cde cDE/cde 
+ + | + + + 
CDe/cde CDe/cde | cde/cde CDe/cde CDe/cde 
+ + + + + 
ceDE/cde CDe/cde cde/cde CDe/cde cDE/cde 
+ + | + + + + + 
CDe/cde cde/cde CDe/cde CDe/cde cde/cde cde/cde CDe/cde 
+ + 1 + 
CDe/cde cDE/cDe |} cDe/cde 
or or 
CDe/cDe cDE/cde 
+ + i + + 
CDe/cDe cde/cde cDe/cde eDe/cde 
+ + + + + 
CDe/cde cde/cde | cede/cde cde/cde CDe/cde 
+ + + + + + + + 
cde/cde CDe/cde | cde/cde cde/cde cde/cde CDe/cde CDe/cde cde/cde 
+ + i 
CDe/cde CDe/cDe | cDe/cde 
or or | 
CDe/cDe CDe/cde 
+ + | + 
cDE/cde cde/cde | cDE/cde 
+ + + + + + + 
cde/cde cde/cde | cde/cde cde cde/cde cde/cde cde/cde 
. + + + 
CD“e/cde cde/cde | cde/cde cde 


+ + | + + + 
cDE/cde cde/cde | cDE/cde cde cDE/cde 
1 + + + 
cde/cde CDe/cde cde/cde cde/cde 
+ + + + + 
CDE/cDe cde/cde CDE/cde cDe/cde cDe/cde 
+ | + 
CY¥dE/cde | CDe/cde CDe/cde C¥dE/cde CDe/cde 
+ + 
C”¥’dE/cde CYdE/cde 


Notes (Continued).—Twins are bracketed, all are monozygous save the second pair in 
family 12. 

The Rh groups of the following families have previously been published: 1, 4, 5, 
18, 21, 27, 28, 29, 45, 46, and 47. They should not therefore be included in any compila- 
tion of genetical studies of the Rh groups 
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In 11 of the families not all the members have beer: tested with anti-f 
but their Rh groups are given from our records because they reveal the 
true Rh genotype of the members who have been tested. 

Without any of the foregoing considerations, association of anti-f with 
the Rh system would be apparent from the families. Ten of them, families 
7,8, 9, 10, 11, 12, 13, 14, 23, and 24, are of the type to show linkage of the 
f antigen with the RA locus, and this they do with no example of crossing- 
over. Applying the u statistics of Fisher,* * as elaborated by Finney,° 
linkage at the 1 in 100 level of probability is established between f and 
Rh, for X(\), which is 36, is greater than 2.33 V x(k), which is 18.9. 


THE RELATION OF THE f ANTIGEN TO THE RA SySTEM 


The antigen recognized by anti-f clearly belongs to the Rh system, but 
it cannot be identified with any of the known members of the 3 series of 
allelic antigens: among the genotypes which are negative with anti-f 
(table 1) are to be found all the known RA antigens (save c’ which we have 
not tested) and, furthermore, all the common antigens are to be found 
there in the homozygous state, DD, dd, CC, cc, EE, and ee. 


TABLE 4 
Tue Cross-OvER THEORY OF FISHER EXTENDED TO INCLUDE f AND F 
The Common Chromosomes Are Known to Be CDeF, cdef, and cDEF: the Table Shows 
the Expected Chromosomes of the Second Order of Frequency. 
COMMON HETEROZYGOTES WITH THEIR APPROXIMATE FREQUENCIES 
33% 12% 11% 
D ‘d 


Recombinations after a Single Cross-Over 
dCeF DCEF dcEF 
Deef DeeF Deef 
dceF dceF 
DCef DcEf 


A fourth allelic series of antigens.-Anti-f, and the antigen it defines, 
will easily fit in with what is already known of the Rh system only if we 
postulate a fourth series of allelic antigens. The antigen f shows pref- 
erences for its partner antigens: it prefers cde, cDe, and cD“e, though it 
may desert at least cde. This is typical of the behavior of the other Rh 
antigens. The crossing-over idea of Fisher, which takes into account 
these preferences, not only accommodates f but also gains support from it. 

Table 4 shows the crossing-over idea extended to include f and its pre- 
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sumptive allele, F. The right order appears to be DCEF. The order 
FDCE would require all cDe chromosomes to lack f, the order DFCE 
would require the majority of cDe to lack f and the order DCFE would 
require a high proportion of cDe to lack f: in fact, all cDe chromosomes so 
far tested have f. 

Thus having allowed the chromosome cDe to dictate the position of 
F in relation to DCE, certain consequences follow, which may be seen 
from table 4. Cde should lack /, which it does; cdE should lack f, which 
it does; CDE should lack f, which it does; and CdE should lack f, which 
it does. (CdF, according to the theory, arises from either of two hetero- 
zygotes, cDE/Cde or CDE /cde, both involving a second order chromosome. ) 
An occasional cde should lack f, and we have found one. A minority of 
cDe should lack f, but we have not found an example among 21 unrelated 
cDe chromosomes. An occasional CDe and cDE should have f, but these 
we have not yet encountered; a possible difficulty in their recognition will 
be mentioned below. 

Other Possible Interpretations.—-The possibility had to be considered 
that f is a new kind of blood-group antigen, a blended antigen, made up 
of ¢c and e and produced only when their genes are together on the same 
chromosome. This was against all our experience. It was, for example, 
difficult to believe that c and e could have a fundamentally different effect 
when present in Cde/cdE and in CdE/cde; it was even more difficult to 
believe that Mr. Cl., who is CDe/cDE, could have made an antibody 
against a blend of two antigens both components of which he himself 
possessed. The possibility was ruled out when a chromosome cde lacking 
f was found (family 19, table 3). 

Two other possibilities, much more difficult to exclude, are: (a) that f 
is an allele of / and e, say e’, and that the common chromosomes are 
cde’, cDe’, CDe, Cde, cDE, cdE, ete.; and (6) that f is an allele of C, ¢ 
and C”, say c’, and that the common chromosomes are c/de, c/De, cDE, 
cdE, CDe, C’De, Cde, etc. They are difficult to exclude because there 
is evidence, to be given below, that most of our anti-e and some of our 
anti-c sera also contain anti-f. 

We feel that the true explanation is unlikely to lie in either of these 
two possibilities for both require a number of very unusual qualifications 
of what we already know about Rh antisera. Further discussion over 
these two possibilities can be postponed since both should be excluded 
when the existence of chromosomes cDE and CDe ‘‘positive’’ with anti-f 
is proved, or when anti-F is found. 


THE CONSEQUENCES OF ANTI-f 


Phenotype Distinctions.-Anti-f does not promise to be of any great 
practical importance in routine Rh blood-group work. Certainly it splits 
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off CDE/cde and CDE/cDe (anti-f positive) from CDe/cDE, CDe/cdE 
and Cde/cDE (anti-f negative); it also separates CdE/cde (anti-f positive) 
from Cde/cdE (anti-f negative). These are hardly day-to-day problems 
though all these distinctions have been made in practice. Indeed the 
third apparent Cde/cdE tested with anti-f was positive and was thereby 
shown to be the very rare CdE/cde, and this was subsequently confirmed 
by the investigation of the family.® 

A curious fact is that this one serum, anti-f, can distinguish between 
DD and Dd as competently as can the four sera anti-C, anti-c, anti-/, 
and anti-e combined. However it should be said that the use of the latter 
four antisera has a slight advantage in permitting a better estimate to 
be made of the error involved in the guess at the genotype; for example 
the phenotype CDe/CDe is less likely to be the genotype CDe/Cde than is 
the phenotype cDE/cDE to be cDE/cdE. Nevertheless, if anti-f becomes 
more generally available, this application will be of some importance in 
testing the blood of fathers of children with hemolytic disease of the 
newborn. 

Failure to Make Anti-f in Volunteers.-Some time ago, in collaboration 
with Dr. B. H. Kirman, 34 CDe/cDE volunteers were given four intra- 
venous injections, at intervals of three months, of 1 ml. of cde/cde blood. 
We hoped to make anti-d but in this we were disappointed; on looking 
back we should have been doubly disappointed, for our tests showed that 
anti-f was not made either. 

The Need to Look for Anti-F.—lf anti-F exists, as we presume it does, 
it will be difficult to identify with the test cells commonly available. It 
can be deduced from the left-hand column of table 2 that anti-/ will give 
the reactions of a mixture of anti-C + anti-. Anti-F could be made by 
cde/cde or by cDe/cde persons; if made by the former it would probably 
be mixed with anti-D. The most hopeful way to find anti-F is to test 
sera thought to contain anti-C + £, or anti-D + C + KE, with cells of the 
genotype cdeF/cdef. So far we have tested five such sera but without 
success. 

Antt-f in Anti-e and Anti-c Sera.—The sample of blood of the genotype 
cdeF/cdef was recognized early in the work on anti-f. The e and c antigens 
in this blood were thereupon titrated with anti-e and anti-c sera and the 
results showed that most of these sera behaved in saline as if they con- 
tained some anti-f, while in albumin their behavior was that expected 
of anti-e and anti-c. That is to say in saline the cdeF/cdef cells reacted 
much less strongly than did cdef/cdef cells, while in albumin the reactions 
were about equal. The effect was most marked with our anti-e sera: 
one of them, Br., absorbed 4 times with CDe/cDE cells left anti-f prac- 
tically free from anti-e. 

This evidence of two antibodies in the anti-e serum Br. explained a 
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puzzle of some years’ standing. From 1946 to the end of 1946 the Blood 
Group Research Unit was doing comparative titrations with this serum. 
There was an interval of seven months, from October, 1948, to May, 1949, 
when no titrations were done. It was noticed, but not understood, that 
before this interval CDe/cDE cells scored more than cDE/cde cells but 
that after the interval CDe/cDE scored much less than cDE/cde. This 
change was reflected also in titrations with cde/cde, CDe/cde, and CDe/CDe. 
The scores are shown in table 5; they are strictiy comparable wherever 
the ratio symbol is used. It is now clear that either the saline anti-e was 
deteriorating, perhaps because of storage, exposing an anti-f component 
or that the successive immunizations of Mr. Br. were stimulating anti-/ 
while making the anti-e incomplete. 

We were somewhat disturbed by the frequent presence of anti-f in these 
sera until we remembered that it was not until anti-C was discovered that 
so many anti-D sera were realized to be anti-D + C. 


TABLE 5 
CHANGE IN THE BEHAVIOR OF THE ANTI-€ SERUM, BR.: SCORES FOR COMPARATIVE 
TITRATIONS IN SALINE 

cDE/cde CDe/cDE 
Before October, 1948 742 864 
205 $3 
Ff FF 

cde/cde CDe/cde CDe/CDe 
Before October, 1948 85 70 : 64 
From May to December, 1949 506 : 288 171 
ff Ff FF 


From May to December, 1949 


It is probably because our anti-e and anti-c sera contain anti-f that we 
have not been able to identify cDE{/cDEF or CDef/CDeF; cDE{/cDEF 
would presumably be agglutinated by the anti-f in our anti-e sera and 
would therefore be called cDE/cde, and similarly CDef;CDeF would be 
labeled CDe/cde if the anti-c serum contained anti-f. This occasional 
wrong diagnosing of cDE/cDE and CDe/CDe would not show as dis- 
crepancies in family studies if the same anti-e and anti-c sera had been 


used on all members of the family. 

It is clear that anti-e and anti-c sera are now exposed to investigation 
from a new angle. 

The most important consequence of the recognition of f is that it must 
be taken into account in any comprehensive consideration of the RA groups. 

Summary.—The reactions of a ‘‘new’’ RA antibody, anti-f/, are described ; 
they can most simply and conservatively be explained by the posulation 
of a fourth allelic series of antigens, f and F. The usual chromosome 
combinations, in white people, are: cdef, cDef, cD“ef, CDeF, cDEF, C’DeF, 
CdeF, cdEF, CD“eF, cD“EF, CDEF, CD“EF, C*deF, CdEF, cDE*“F, 
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C“Del’, and CYdEF. Certain less usual combinetions are expected: of 
these cdeF has been found, cDeF, CDef, and cDF/ are yet to be found. 
Anti-/ is commonly present in anti-e and anti-c sera. 
Anti-F is presumed to exist but it will be somewhat difficult to identify. 
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FURTHER STUDIES ON THE EFFECTS OF CARBON DIOXIDE 
AND OX YGEN ON THE FREQUENCY OF X-RAY INDUCED 
CHROMOSOME ABERRATIONS IN TRADESCA NTIA* 


By Howarp A, SCHNEIDERMANT AND Epwarp D. KING 
THe BroLoGicAL LABORATORIES, HARVARD UNIVERSITY 
Communicated by Karl Sax, Tune 8, 1953 


Previous studies have shown that exposure to carbon dioxide during ir- 
radiation markedly increases the radiosensitivity of cells.'~° It has also 
been demonstrated that the concentration of oxygen present during irra- 
diation has a profound effect on the radiosensitivity (e. g., Giles and Riley®). 
The experiments to be reported examine the interaction between these 
two gases and further define the mode of action of carbon dioxide in increas- 
ing radiosensitivity. 

Matertal and Methods.—-Tradescantia paludosa Anderson and Woodson 
was used in all experiments. Cytological observations of chromosomal 
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aberration types were made 5 days after irradiation. The principal tynes 
of aberrations analyzed at this period are double interchanges (dicentrics 
and centric rings). 

The experimental methods were essentially the same as those described 
previously. The gases employed were obtained in commercial cylinders 
and assayed as follows: Nitrogen (Airco), 99.5% Ne plus 0.5% O2; Car- 
bon Dioxide (Liquid Carbonic Co.), 99.8% COs, 0.05% Ne, 0.009% On, 
0.07% HO. The inflorescences were placed inside 3'/»-liter lucite cham- 
bers initially filled with air. In one series of experiments the chambers 
were first evacuated and then refilled with specific gases at specific pressures. 
In a second series, the chambers were not evacuated; specific positive 
pressures of the gas under consideration were superimposed on the atmos- 
phere of air initially present in the chamber. 

TABLE 1 
FREQUENCIES OF DICENTRICS AND CENTRIC RINGS INDUCED IN TRADESCANTIA MICRO- 
SPORES EVACUATED AND THEN X-RAYED IN VACUO OR IN MIXTURES OF AIR AND CARBON 
Droxipe. ALL DosaGEs ADMINISTERED AT A CONSTANT INTENSITY OF 40 R/MIN. FoR 
10 MINUTES 
EXPERIMENTAL CONTROL 
GAS MIXTURE DICENTRICS GAS MIXTURE DICENTRICS EXPERIMENTAL 
COn at wie pHERES), AIR OF CONTROL 
0) 0 0.9 1 6 14 
0.33 0 2.0 0 117° 
0.33 0 1.9 0 1062 
I 0 1.0 0 x ti be 
1° 1 23.8 1 


] 
| 
0.33° 1 13.4 | 7. 184 
] 
7 


7 309 
“ Differences of this order of magnitude between experimental and control with 
aberration frequencies shown here are not statistically significant for the number of 


chromosomes counted (i.e., 3000). 
+ These preparations were not subjected to preliminary evacuation 


After the desired gas mixture had been established within the chamber, 
the latter was sealed and was maintained at 25°C. in the dark for a speci- 
fied period. It was then placed in the x-ray machine alongside an iden- 
tical chamber containing the dosimeter. Irradiation was performed at 
room temperature in subdued light. This procedure in each experiment 
was carried out twice—once with the experimental gas (usually CO.) and 
once with the appropriate control gas (usually air). Within 30 minutes 
after irradiation the chambers were opened and the influorescences re- 
turned to air. 

Ten to 15 inflorescences were exposed in each chamber. Five days fol- 
lowing irradiation as many slides as possible were prepared. Three 
thousand chromosomes were counted in each control and experimental 
group. The scoring was performed by one investigator (E. D. K.). 
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Results.—In an effort to examine the interaction between carbon dioxide 
and oxygen in conditioning radiosensitivity, the following experiments 
were carried out. The chambers containing the influorescences were evac- 
uated for 3'/, minutes and then either maintained in the evacuated state 
or refilled to a desired pressure with a particular gas. Fifteen minutes 
after the initial evacuation the chambers containing inflorescences were 
irradiated for 10 minutes at a constant rate of 40 r/min., after which the 
chambers were slowly decompressed and the inflorescences returned to 
air. The results summarized in table | permit the following conclusions: 

1. Exposure to a vacuum during irradiation markedly decreases the 
aberration frequency below that encountered in control preparations ir- 
radiated in air. This result is in agreement with numerous observations 
by other workers.® 

2. Exposure to '/; of an atmosphere of carbon dioxide plus 1 atmosphere 
of air greatly increases the aberration frequency above that of controls ir- 
radiated in air. 


TABLE 2 
FREQUENCIES OF DICENTRICS AND CENTRIC RINGS INDUCED IN TRADESCANTIA MICRO- 
SPORES EXPOSED TO X-RAYS IN AIR AND THEN COMPRESSED WITHIN 18 SECONDS WITH 
CarBON Droxipk. ALL DosaGes ADMINISTERED AT A CONSTANT INTENSITY OF 200 
R/MIN. FoR 2 MINUTES 
GAS MIXTURE DICENTRICS EXPERIMENTAL 
(ATMOSPHERES) AND RINGS, AS PER CENT 
COs AIR / OF CONTROL 
Experimental 0.33 1 0) 962 
Control l 


7 
7.3 


“ Differences of this order of magnitude between experimental and control with aberra- 
tion frequencies shown here are not statistically significant for the number of chromo- 


somes counted (i.e., 3000). 


3. irradiation in the presence of as much as | atmosphere of carbon 
dioxide in the absence of air is equivalent to irradiation in a vacuum. 

Thus it appears that carbon dioxide acts synergistically with oxygen in 
increasing radiosensititivity. But the effect of carbon dioxide on radio- 
sensitivity disappears when oxygen is absent. In short, it appears that 
carbon dioxide acts by modifying radiochemical reactions that involve oxy- 
gen; it affects the ‘indirect action”’ of x-rays and not the “direct action” 
(cf. Sparrow and Rubin’). 

In a second series of experiments we sought to determine whether carbon 
dioxide acts at the instant of irradiation or subsequent to radiation. To 
this end, the inflorescences were first irradiated at 200 r/min. for 2 minutes. 
Then, within 18 seconds after irradiation, they were compressed with '/; 
of an atmosphere of carbon dioxide. After 60 minutes of exposure to the 
carbon dioxide, the chamber was decompressed and the inflorescences re- 
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turned to air. A control group was treated in identical fashion save for 
omission of the compression with carbon dioxide. Table 2 shows that post- 
irradiation treatment with carbon dioxide had no detectable effect. Con- 
sequently, we conclude that carbon dioxide exerts its action at the time of 
irradiation. Moreover, since carbon dioxide and oxygen act synergisti- 
cally, it follows that oxygen acts at the time of irradiation. 

Figure | records the effects of mixtures of one atmosphere of air plus 
varying pressures of carbon dioxide. The results of the present study are 
plotted along with those previously reported.® In all these experiments the 
air-filled tanks and the inflorescences therein were compressed with carbon 
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PRESSURE OF ADDED GAS (ATMOSPHERES) 
FIGURE 1 


Curve showing the effect of mixtures of one atmosphere of air plus varying pressures 
of carbon dioxide or nitrogen on the incidence of x-ray-induced dicentrics and centric 
rings in Tradescantia microspores 


dioxide (or nitrogen in the case of controls) for | hour prior to irradiation. 
The inflorescences were then irradiated for 10 minutes at a constant rate 
of 40 r/min. Within 30 minutes after irradiation the inflorescences were 
returned to air. As previously noted® the results show that the greater the 
carbon dioxide concentration (up to at least 1'/; atmospheres of carbon di- 
oxide), the greater the increase in aberration frequency. Moreover, car- 
bon dioxide concentrations of less than | per cent of an atmosphere appear 
to cause a significant increase in aberration frequency. 

The relation of carbon dioxide pressure to radiosensitivity in Trade- 
scantia (Fig. 1) differs strikingly from that observed in fern spores and 





838 GENETICS: SCHNEIDERMAN AND KING Proc. N. A. 5S. 


Drosophila eggs by Zirkle.2~* A complex ‘‘wave-form’’ relation between 
‘“arbon dioxide concentration and radiosensitivity was reported by this 
investigator: certain concentrations of CO, increased radiosensitivity 
and other concentrations decreased radiosensitivity. Thus, in the case of 
Drosophila eggs, carbon dioxide tensions between 16 and 21% and above 
64% of an atmosphere exerted a protective effect against radiation damage.’ 

Our experiments with Tradescantia reveal no such “‘wave-form’’ rela- 
tion between carbon dioxide concentration and radiosensitivity. To the 
contrary, radiosensitivity underwent a steady, albeit non-linear increase, 
as the carbon dioxide concentration increased. Thus carbon dioxide affects 
chromosome aberrations in Tradescantia in a manner strikingly different 
from its effects on the radiosensitivity of fern spores and Drosophila eggs. 
Indeed it appears that carbon dioxide may have widely different influences 
on different radiobiological systems. 

In conclusion we suggest that in Tradescantia carbon dioxide is affecting 
the “indirect action’ of x-radiation. Furthermore, carbon dioxide here 
appears to affect radiosensitivity only at the time of irradiation. Since 
‘arbon dioxide and oxygen act synergistically these observations are con- 
sistent with the view that oxygen exerts its action on radiosensitivity at 
the time of irradiation. Comparison of our data for Tradescantia with 
Zirkle’s data for Drosophila eggs and fern spores leads us to conclude that 
in different radiobiological systems carbon dioxide may act in very differ- 


ent ways. The information at hand permits no conclusion as to the pre- 
cise biochemical nature of the carbon dioxide effects. 


* This work was financed by aid from the Office of Naval Research, Contract No. 
N8onr-73100, and by an institutional grant to Harvard University from the American 
Cancer Society, Inc. We wish to thank Professor Karl Sax for his helpful suggestion 
and Professors Carroll M. Williams and Raymond E. Zirkle for reading the manuscript. 
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GENE ACTION IN RELATION TO GROWTH AND DEVELOPMENT. 
I, PHENOTYPIC VARIABILITY 


By F. W. WENT 
CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA 
Read before the Academy November 10, 1952; communicated May 25, 1953 


Among the most fundamental problems in biology belongs the action of 
the gene. The gene, however, is elusive since no direct approach is pos- 
sible: the gene can only be recognized by its effect. These effects may 
be used either in their final expression—-the heritable characters—and 
then we are in the domain of classical genetics. Or we can try to infer from 
the final and intermediate steps of gene action what lies at the beginning of 
this chain. This is the approach of both the biochemist, the embryologist, 
and the physiologist. The successes of the biochemist in approaching the 
action of the gene are spectacular, such as in the analysis of biosynthesis 
in Neurospora. 

The physiologist or embryologist cannot hope to contribute much in the 
direction of the chemistry of gene action, since they are dealing with proc- 
esses rather than specific biochemical reactions. But they can contribute 
in several other directions. On the one hand they are dealing with growth 
and development, which cannot be expressed in chemical terms as yet but 
which are of such importance in inheritance. On the other hand they 
stress the quantitative side of biological processes more than the biochem- 
ist, who, in his analysis of the action of the gene, largely stresses the path- 
ways of gene action, 

The following discussion will deal primarily with physiological observa- 
tions which have a bearing on the quantitative side of gene action. This 
problem has many different aspects. The problem of practically absolute 
uniformity in gene reproduction will not be dealt with. 

In gene action, we find absolute reproducibility in ontogeny, in which 
the sequence of events and of organs meets only very seldom with a slight 
deviation in timing, but the number and order of, for example, flower parts 
in the pea flower are practically without variability. 

Not only the qualitative aspects, but also the quantitative side of growth 
is genetically controlled. Here we meet with a much greater variability: 
Even within a pure line or clone there is considerable variation in size or 
weight or any other quantitative character.' Two factors contributing 
to this phenotypic variability have been recognized. In the first place, 


the environment has a great effect on the expression of quantitative char- 


acters. Hormonal control of growth insures proportionality between the 
different parts of the organism: when for some reason one part is small, 
most other organs are proportionately reduced in size. The absolute di 
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mensions of the organism are therefore less directly gene-controlled. 
The proportions, however, are under gene control. Yet the size of an 
organism under a specific set of growing conditions, and particularly its 
maximal size, is definitely an inheritable character. 

In the second place, all cellular reactions occur on a microscopic scale, 
where the molecular structure of matter starts to interfere with the regular- 
ity of classical statistical physics. This has been pointed out by many 
physicists,?~* and was very clearly expressed by P. Jordan in 1936: “‘Wide 
applicability of macrophysical laws to organisms is unquestioned; yet 
there are good grounds to believe that the basic life processes do not belong 
in the realm of macrophysics any more. The macroscopical processes 
occurring in the larger parts of the body are governed by considerably more 
delicate processes. The ultimately controlling reactions are generally of 
atomic dimension.'’ This leads Jordan to his ‘‘amplifier’’ theory of or- 
ganisms, in which he stresses the fact that in organisms the basic controlling 
“reactions of atomic dimensions” such as gene reactions, are amplified to 


macroscopic observability by cellular processes. ‘‘Genetics provides the 
most convincing and widest basis for the thesis, that organisms are not 
macrophysical, but microphysical systems.” 

Continuing his argument, Jordan then states that unstability, or in- 
constancy, is a fundamental characteristic of quantum physics and that 
actually inconstancy occurs in biology as well, as the statistical nature of 


the Mendelian laws shows. From this, one might conclude that when 
genes control development, one would expect a fundamental variability in 
organisms. This thesis was actually elaborated by Otto Rahn° who cal- 
culated from the observed variability in biological material (in his case 
bacteria) the numbers of genes involved. He obtained figures ranging 
from 100 to 1000, which agree with estimates made of the number of genes 
occurring in organisms using other premises. 

Probably few other biologists rationalized their acceptance of pheno- 
typic variability as due to control of growth and development by atomistic 
reactions in the sense of Jordan (see, e.g., the criticism of Biinning). Yet 
we find that not only geneticists, soil scientists and agriculturists, who have 
to work with variable material of which they cannot reduce the variability, 
but also biologists, have adopted statistical analysis as the general method 
of coping with variability in biological material. This indicates to what 
extent phenotypic variability was accepted as a basic property of living 
matter by biologists as a whole. If phenotypic variability had been con- 
sidered as being largely due to environment, more serious efforts would 
have been made to control the external environment of growing plants. 

Experience gained in the Earhart Plant Research Laboratory shows that 
environment rather than the atomistic nature of biological reactions is 
responsible for phenotypic variability. It was generally observed that plant 
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material grown under the controlled conditions of this laboratory was less 
variable than similar material grown in an ordinary greenhouse. This 
was of practical importance in the design of experiments, in which fewer ob- 
jects could be used to obtain the same significance, even to the extent that 
only four plants are now used per treatment in many experiments. Jus- 
tification for this reduction in size of sample in experiments is found in 
figure |. This shows the effects of an increase in uniformity of growing 
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FIGURE 1 
Coefficient of variability of dry weight of tomato plants grown in 
either an ordinary greenhouse (/eft) or in an air-conditioned green- 
house (right). The latter plants were kept for 2, 3, 9, or 10 days 
(abscissa), prior to weighing, in controlled light rooms with even more 
uniform environment 


conditions on variability in weight of young tomato plants, when they are 
grown with different degrees of environmental control. In most plant 
material we find that the coefficient of variability-—-when the standard de- 
viation of the mean is expressed in per cent of the average—amounts to 
anywhere between 20 and 40 per cent. In the Earhart greenhouses this 
variability is reduced, usuaily to about 15 per cent. Yet, when similar 
plants are grown in temperature- and light-controlled rooms, their coef- 
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ficient of variability is decreased to about 8 per cent. The details of these 
determinations will be published later. 

At present, the question cannot be answered how much of the remaining 
8 per cent variability is innate, and how much is imposed by inconstancy 
of growing conditions. In the following observations the coefficient of vari- 
ability in pea plants decreases to below 4 per cent (table 1) and it is very 
likely that further work will reduce this still more. 

We can view this problem of variability from many other angles. If 
growth were controlled by processes whose rates are subject to appreciable 
random fluctuations, the variability of a batch of plants should remain the 
same or perhaps even increase as they continue to grow. This is not ac- 
tually the case as field observations will show: Even though germination 
may be variable in a wheat or corn field and the seedling stand may consist 
of fairly uneven plants, a well-tended field of mature wheat or corn is 


TABLE 1 


LENGTH IN MM. OF PEA PLANTS, GROWN IN A 16-Hr. PHoTopeRIop aT 17°C. IN 
DIFFERENT INTENSITIES. EACH Group CONSISTED OF 9 PLANTS; IN EACH CASE THE 
STANDARD ERROR OF THE MEAN Is GIVEN AND THE COEFFICIENT OF VARIABILITY 
vinco, 1500 FT.-c vinco, 1000 Fr.-c vinco, 500 rT.-c KRONBERG, 1500 rT.-c 

F COEFF COEFF COEFF COEFF 


oF OF OF 
DATE LENGTH VARIA. LENGTH VARIA LENGTH . LENGTH VARIA 


23 VII 34.3 + 1 11.5 35.8 d i 37.5 EF d 32.4 e 2. 
25 VII 51.8 + 2 10.9 49.3 p é 53.3 2.2 ae 45.9 16. 
31 VII 107.2 + 2 6.8 107.8 2.3 5.9 104.% p B 80.9 18 
7 Vill 189.443 4.5 191.3 é 163.! : 142.4 15.3 
+ 3 é 
+ 3 


2 
2 
9 
5 
a 
7 


12 VIII 235 3.6 248.9 5 205 f 
13 VIII 249.6 3.4 2 223. - §.% 194.9 
Dry wt. in 
mg 870 2% 5.2 ‘ 415 ¢ j 762 
Actual increase 
in length 
from 31 VII 
to 13 VIII 142.3 2.6 ‘ ° ceepee ee ws 1414.0 + 3.6 


surprisingly uniform. The following set of measurements substantiates 
this field observation. Groups of ‘‘Vinco” and ‘‘Kronberg’’ peas were 
grown in three different light intensities: 500, 1000, and 1500 ft.-c., and 
as they grew, the distance to the artificial lights was adjusted so that their 
tops continued to receive the same intensities. 

To start with, groups of nine plants were selected to be as uniform as 
possible in respect to size and form. They were measured from time 
to time, and table 1 and figure 2 show the results. The variability 
was small to start with in all three groups; in absolute figures the vari- 
ability increased least in the peas growing in the highest light intensity. 
But the relative variability decreased very much—to less than one-half 
of what it originally was. Because the variability in these peas was so 
small to start with, the causes for its decrease are not so clear. In the 
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Kronberg variety which was more variable at the start, variability remains 
much greater, but not because the longest plants continue to grow fastest. 
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FIGURE 2 


Length of Vinco pea plants grown in a 16-hour photoperiod at 
17°C. in different light intensities. Variability of each group of 
plants is indicated by a 2 * standard error 


For it is seen that the variability of the length increase during the last 
fourteen days is less than the previously existing variability (the largest 


TABLE 2 
GROWTH OF PEA PLANTS DEVELOPING FROM SEEDS OF DIFFERENT WEIGHT AVERAGE 
OF 25-35 PLANTS PER GROUP 


TOTAL 
LENGTH 
auG. 18-21 AUG. 21-26 AUG. 26 
180 130-199 ye 36.5 73.5 134.2 
230 201-246 22.4 37.5 7238 137.5 


” 
330 298-377 y 31.2 84.0 136.0 


WEIGHT OF DRY SEED, MG GROWTH IN MM 
AVERAGE RANGE 


differences on July 31 were 54 mm.; on August 13, 69 mm., whereas for 
the length increments they were only 44 mm.). 
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These measurements show that the growth rate of these peas does not 
depend on the preceding conditions, and that differences within a well- 
selected strain are neither genetically nor phenotypically conditioned. 
Differences between individuals do not become magnified, but decrease 
relatively. Therefore we have to consider the possibility of a self-regula- 
tory mechanism which adjusts the growth rate of these plants. This also 
follows from the observation that originally the plants from bigger pea 
seeds grow slower than those from smaller ones; after some time they 
catch up with the plants from the smaller seeds and reach the same over- 
all size, as shown in table 2. Similar results were obtained when plants 
were arranged according to their size and growth rate after one week’s 
growth. Those growing fastest reduce their growth rate in the following 
eleven days to below that of the originally slower growing peas to such an 
extent that their final size was within 0.7 per cent of each other, whereas 
originally this group was 17 per cent larger (see Fig. 3). Even the very 
smallest peas (only 45 per cent of the largest) caught up with the others. 


TABLE 3 
COEFFICIENT OF VARIABILITY OF TOMATO PLANTS WHEN GROWN IN SUBOPTIMAL, 
OPTIMAL AND SATURATING LIGHT INTENSITIES AT 17°C. IN A 16-Hr. PHOTOPERIOD. 
THE VARIABILITY IS BASED ON THE Dry WEIGHTS OF 20 CONTAINERS WITH 20 PLANTS 
EACH, AS IN FIGURE | 
LIGHT INTENSITIES, FT.-C COEFFICIENT OF VARIABILITY AVERAGE, “% 
500-700 12.0; 12.5; 14.2 12.9 
900 8.4; 8.9; 8.9 8.7 
1100-1500 9.8; 10.2; 10.7; 11.6; 12.0; 13.3 11.3 


This was especially evident when a correction was made for the physiological 
age of the small plants—which is about two days less than that of the other 
peas—then their growth rate fell within the range of the originally taller 
plants. 

In connection with variability of experimental material another obser- 
vation has to be recorded. Under optimal growing conditions, when rate 
of elongation is greatest, the coefficient of variability is least (see Fig. 4). 
This was found in young pea plants grown at different temperatures. It 
is interesting to note that in the youngest stage, the growth rate is unaf- 
fected by temperature over the 15.5° to 26°C. range. Only below 15° is it 
controlled by a process with a high temperature coefficient. And it is 


just where growth is presumably limited by a chemical process that varia- 


bility increases. 

In tomato experiments it was also found that under most nearly optimal 
growing conditions the variability of the plants was the least; therefore, 
also at the greatest rates of growth. For instance, the average dry weights 
of tomato plants grown like those in figure 1 showed a coefficient of varia- 
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FIGURE 4 
Growth rates (left ordinate, mm./day) of Vinco peas at different 
temperatures (abscissa) for the period of 5-13 days after planting 
(crosses and solid line) and 13-20 days after planting (plusses and 
dotted line). Coefficient of variability (right ordinate) of the 
growth rate from 5 to 13 days shown in circles with broken line. 
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bility of 12.9 pet cent when grown at suboptimal light intensity, and of 
11.3 per cent at supraoptimai intensities, but of 8.7 per cent at the optimal 
intensity (table 3). 

We can draw a number of conclusions from these observations: 

(1) Variability in the material described is smaller than ordinarily sup- 
posed to be typical for genetically uniform plants. The selected examples 
are not exceptional, however, and variability of plants in the Earhart 
Laboratory is generally much less than when similar plants are grown un- 
der uncontrolled conditions. This is due to two independent effects: 
(a) through proper air distribution in the greenhouses and controlled light 
rooms, the individual plants are more nearly all under the same growing 
conditions, which they are not, e.g., on an ordinary greenhouse bench 
where the plants in the center have much less air circulation and therefore 
are subjected to higher temperatures during the day, and (4) under optimal 
growing conditions variability is also smaller (Fig. 4); with the great vari- 
ety of growing conditions in the Earhart Laboratory plants in general can 
be grown under more nearly optimal conditions and consequently will be 


less variable. 
(2) From figure | it can be seen that in genetically uniform material, 
the upper limit of any innate variability is around 5 per cent. The addi- 


tional variability usually observed is imposed largely by uncontrolled 
growing conditions. There are two explanations for this low variability. 

(3) On the one hand, it is possible that low variability is an innate 
characteristic of biological material. This seems to be the view of both 
Bohr and Schrédinger. Bohr’ states: ‘The explanation of the prop- 
erties of living organisms certainly lies in their specific organization, in 
which typical atomic traits are interwoven with mechanistical traits to an 
extent which has no counterpart in inorganic nature.”’ 

(4) On the other hand, we could assume that growth is variable like 
any physicochemical process, but that superimposed on it is an auto- 
regulatory process, which corrects any deviation imposed on it by chance. 
This would seem to be substantiated by the data of figures 2 and 3, from 
which it follows that the originally slower growing plants grow faster after- 
ward, but in later development this was less pronounced: It took the peas 
(shown in Fig. 3) which were growing slowest from August 22 to 26 three 
days to catch up with the fastest growing ones and even then they did not 
surpass them (originally 23.8 and 27.4 mm./day; next three-day period, 
24.0 and 28.7 mm./day; last four-and-a-half days, 28.7 and 28.7 mm./ day). 
Therefore, although there is a tendency for growth of the slower growing 
plants to pick up, there is no immediate response in all cases, and the self- 
regulatory mechanism is not very intimately tied up with the growth 
process itself. 

(5) Is it possible to imagine a process which controls growth which 
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could explain the phenomena described in this paper? From figure 4 we 
can conclude that at the opt‘mal growth rate a diffusion process, largely 
independent of temperature, controls growth in young peas. Along a very 
different line of reasoning,’ it was concluded that diffusion within a cell be- 
tween the nucleus and cytoplasm is the limiting process in growth. 

(6) Such a diffusion process is slower when the cells or organs are larger. 
Thus when faster growth has resuited in larger cells in the growing point, 
this would automatically reduce the growth rate. This is just the opposite 
of control by chemicals, where a larger original structure would tend to in- 
crease the supply of hormones or other growth factors and would result in 
faster growth. We would expect, then, the larger individuals to be both 
the result of and to result in a greater growth factor supply. Therefore 
through chemical control of the growth rate, original small differences in 
size would tend to become magnified, whereas just the opposite occurs 
(Figs. 2 and 3 and table 2). By the simple assumption that growth is 
controlled by an intra- or intercellular diffusion process, we can explain a 
decrease in growth rate of the growing points which have become larger 
than the average and vice versa. 

(7) How can we tie the observed phenomena in with genes? The gene 
theory was based upon the discontinuity of the hereditary process. The 
distribution of characters over the offspring is a typical quantized phenom- 
enon, explained by the presence in cells of a few discrete units or genes. 
The same genic concept can account for effects of ionizing radiation on 
organisms (see, e.g., Timofeeff-Ressovsky and Zimmer’); through these 
radiation effects the approximate size of the gene can even be calculated. 

There is, of course, no doubt that growth and the growth rate in plants is 
dependent upon the hereditary mechanism. For each pea variety, al- 
though very uniform within itself, differs from most other varieties in its 
growth rate. Yet, since the growth rate is so remarkably uniform from 
plant to plant, we must conclude that either the growth rate is only very 
indirectly dependent upon genes, or that growth is controlled by hereditary 
units which are present, not in twofold, but in very large numbers within 
each cell. 

The former possibility resolves itself into two mechanisms: Either the 
laws of quantum mechanics do not hold for the reproduction and physi- 
ological effects of genes, as Schrédinger* suggests, or the gene becomes 
multiplied in each cell before it exerts its effect. In practice, the latter 
mechanism becomes almost identical with the assumption of polygenes 
controlling quantitative processes inside the cell. But whereas formerly 
the nature of polygenes could only be inferred from the distribution of the 
quantitative characters over the offspring, it should now become possible 
to study this problem by the quantitative effect of the hereditary mecha- 


nism upon the growth process within an individual. In subsequent papers 


this will be analyzed in more detail. 
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The physiologist and embryologist are now in the uncornfortable position 
where they will have to explain why genes, which were just proved by the 
geneticist and biophysicist to be following quantum mechanics, do not im- 
part the degree of variability upon growth and development which one 
would expect from quantum mechanics. We can view this difficulty also 
from another standpoint, however. If the uncertainty principle of quan- 
tum mechanics did hold for gene reproduction and gene effect, organic life 
in its present form would be impossible, because there would be only statis- 
tical reproducibility of organisms, and very soon life would have reverted 
to the more stable and probable inorganic state. Only because develop- 
ment does not follow the laws of quantum mechanics can the less probable 
condition of life be continued indefinitely. The reduced variability of 
genetically uniform plant material when grown under controlled conditions 
is therefore nothing exceptional; it is the logical consequence of the rigid 
laws of organic development which fall outside the laws of statistical or 
quantum mechanics. This has been expressed before by Schrédinger: 
“The unfolding of events in the life cycle of an organism exhibits an ad- 
mirable regularity and orderliness unrivaled by anything we meet in in- 
animate matter.” 

Summary.—The major part of the commonly encountered variability 
of genetically uniform material (phenotypic variability) is not due to statis- 
tical fluctuations of the numbers of molecules on which development de- 
pends, but is largely caused by inconstancy of and irregularities in the ex- 
ternal environment. The genetic homogeneity of pure lines can be trans- 
lated under completely controlled conditions into phenotypic uniformity 
of plant material which exceeds anything observed thus far. The more 
optimal the growing conditions are, the smaller is phenotypic variability. 
When the standard deviation is expressed as per cent of the mean, varia- 
bility in peas and tomatoes is regularly reduced to 2 per cent, and in the 
most favorable material it was as little as 1.2 per cent. Due to an auto- 
regulatory mechanism phenotypic variability decreases in the course of 
development. All these facts lead to the conclusion that the laws of or- 
ganic development fall outside the laws of statistical or quantum mechanics. 

! Johannsen, W., Elemente der exakten Erblichkeitslehre, Jena, Gustav Fischer, 1909. 

2 Bohr, N., ‘Licht und Leben,” Naturwiss., 21, 245-250 (1933). 

§ Jordan, P., Anschauliche Quantentheorie, Springer, Berlin, 1936. 

4 Schrédinger, E., ‘What is life. The physical aspect of the living cell,”’ Univ. Press, 
Cambridge, 1945 

®’ Rahn, O., “Die durch chemische Gesetze bedingten Variationen der Lebewesen,” 
Bioch. Zeit., 284, 40-62 (1936) 

® Binning, E., ‘‘Physikalisch-chemische Grundlagen der biologischen Vorgange,”’ 
Fortschr. Bot., 8, 154-165 (1939). 

7 Went, F. W., “Physical Factors Affecting Growth in Plants” (in press), 1958. 

’ Timofeeff-Ressovsky, N. W., and Zimmer, K. G., Biophysik I. Das Trefferprinzip 
in der Biologie, Hirzel, Leipzig, 1947. 
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TRAVEL TIMES OF LONGITUDINAL WAVES FROM 
, SURFACE FOCI 


By B. GUTENBERG 
CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA * 
Communicated June 8, 19538 


Tables and graphs for travel times of seismic waves propagated through 
the earth are based on earthquake records. This requires accurate knowl- 
edge of the longitude and latitude of the corresponding epicenters, focal 
depths, and origin times. It is highly desirable to compare the resulting 
travel time curves and tables with results derived from experiments in 
which all four fundamental quantities are accurately known from independ 
ent data.! Sets:nologists are indebted to scientists and organizations for 
measurements of records and other corresponding information from which 
accurate travel times of longitudinal waves originating at known surface 
sources in the Pacific area could be found. 

Residuals for 33 observed direct longitudinal waves, P, were calculated 
relative to travel times of Jeffreys and Bullen’ for average continental sur- 
face foci, as well as relative to similar tables by Gutenberg and Richter.* 
The following average residuals (with standard errors) from the Jeffreys- 
Bullen tables’ are found in four ranges of distances containing most of the 
observations: 


Distances, deg 35-45 50-65 68-79 80-91 
Residuals, sec ~-1.7+0.6 -3.9+0.4 -1.32+0.6 -1.1+0.7 


It has been suspected previously! that travel times from surface sources in 
the Pacific basin are about 2 seconds shorter than the corresponding times 
from continental surface foci due to differences in velocity in the surface 


TABLE 1 
CORRECTIONS TO THE TRAVEL TIMES OF JEFFREYS-BULLEN? (JB) AND GUTENBERG- 
RIcHTER® (GR) FOUND FROM PACIFIC SURFACE SouRCES (.S, REDUCED TO CONTINENTAL 
STRUCTURE) AND THE EARTHQUAKE OF JULY 21, 1952 (£); FoR Dertats, Ske TEXT 
A = ANGULAR DISTANCE 


30 10) 50 60 70 & 90 100 deg 
4] 0 0 +] +1 sec 
+] 1 0 + ’ 0) sec 

l ‘ 1'/ L sec 


3 2 
0 21/» 3 


-2 
3 
1 


; 
-2 2 2 sec. 
layers. Consequently, 2 seconds were added to all residuals from Pacific 
surfaces. Averages based on graphs of the individual residuals similar to 
figure | are given in the lines marked 5S in table 1. 
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Records of the Kern County, California, earthquake of July 21, 1952, 
gave an opportunity to calculate travel times from a well-located conti- 
nental source. Over 150 stations all over the world responded to a request 
to mail originals or copies of their records to Pasadena. Readings of arrival 
times from 40 additional stations are available at the time of writing. 
Careful determination of the elements of this shock (to be described else- 
where) resulted in an origin time of 11"52"14.°3 with a standard error of 
about 0.2 second and an epicenter at 35°00’ North Longitude, 119°02’ 
West Latitude, probably correct within 3 km. 

For all measured and reported arrival times of the direct longitudinal 
wave P travel times were calculated as well as the corresponding geocentric 
distances. Residuals from the travel times of Jeffreys-Bullen’ for their 
focal depth ‘*h = 0.00" (supposed to correspond to 33 km.) and from the 
travel times of Gutenberg-Richter‘ for a focal depth A of 25+ km. were cal- 
culated for each reading and plotted as a function of distance in the re- 
spective two sections. Half the travel time differences between the curves 
of Jeffreys-Bullen for ““h = 0.00” and ‘‘surface focus’’ are also indicated in 
the lower part of figure 1. Since the earthquake of July 21, 1952, probably 
had a focal depth of about 16 km., the deviation of the observed times from 
this curve corresponds approximately to the corrections to be applied to 
the Jeffreys-Bullen tables. Similarly, the differences in travel time be- 
tween h = 25 and h = 16 km. is indicated in the upper part of the figure. 
It has been found previously® that the zero point of the Gutenberg-Richter 
travel times is about | second too early. To adjust for this difference 
1 second would have to be added to all residuals in the lines GR of table 1. 

Table 1 shows that the assumptions made, especially those for the dif- 
ference between Pacific and continental velocities in the earth’s crust, are 
good approximations. 

From table | and the corresponding graphs (including figure 1) correc- 
tions to the /B and the GR curves were estimated for every degree and ap- 
plied to the corresponding tables for surface focus. Table 2 contains a 
smoothed revised travel time curve for continental surface focus on the 
basis of these results. 

The data for the Pacific surface sources contained five instances of P” 
(longitudinal waves through the outer and inner core) observed at distances 
between 110 and 140 degrees. Three of them showed two phases between 
5 and 15 seconds apart. The earlier phases probably correspond to early 
P” waves observed in earthquake records; ‘“‘these early readings are prob- 
ably not due to longitudinal waves approximately following the laws of geo- 


metrical optics; a phenomenon of diffraction is probably involved.’ 
For the five later impulses the residuals based on the tables of Gutenberg 
and Richter’ for surface focus are very small if the hypothetical effect of 
Pacific structure is considered. 
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SKP (along one path through the mantle longitudinal, along the other 
transverse; along the path through the outer core longitudinal) was re- 
ported by one station a few seconds earlier than calculated from the 
tables* * One P’P’ (longitudinal wave which has passed twice through 
the whole earth, travel time about 40 minutes) had a residual of —8 seconds 
from the Gutenberg-Richter tables" (corresponding to —4 seconds of 
P’) and one PKKP (reflected inside the core from the core boundary) 
observed near a distance of 110° had a travel time exactly in agreement 
with the tables of Gutenberg-Richter."! 

The errors of the average travel time curves for direct longitudinal 
waves are now probably smaller than the effects of the local structure near 
the source and the station. The corresponding differences in travel time 
depend mainly on the depth of the Mohoroviti¢ discontinuity at the 
source (in shallow earthquakes) and the station and may amount to 3 
seconds or even more at distances of a few degrees and at least 2 seconds at 
greater distances. 

* Division of the Geological Sciences, Contribution No. 637. 

' Gutenberg, B., and Richter, C. F., Trans. Am. Geophys. Union, 27, 776 (1946). 

2 Jeffreys, H., Geaphys. Suppl. Mon. Not. R. Astron. Soc., 4, 522 (1939). 

4 Gutenberg, B., and Richter, C. F., Gerlands Beitr. Geophys., 54, 97 (1939). 

‘ Gutenberg, B., and Richter, C. F., [bid., 43, 82 (1934). 

5 Gutenberg, B., Bull. Seism. Soc. Am., 41, 145-146 (1951) 

6 Gutenberg, B., and Richter, C. F., Gerlands Beitr. Geophys., 54, 113 (1939). 

7 Tbid., p. 115, table 19, column C. 

8 [bid., p. 126. 

9 Jeffreys, H., Geophys. Suppl. Mon. Not. R. Astron. Soc., 4, 609 (1939). 

© Gutenberg, B., and Richter, C. F., Gerlands Beitr. Geophys., 54, 131 (1939) 

11 Thid., p. 127 


APPROXIMATE METHODS IN MATHEMATICS* 
By M. HERZBERGER 
RESEARCH LABORATORIES, EASTMAN KODAK COMPANY 


Communicated by C. E. K. Mees, June 15, 1953 


The following simple ideas have proved advantageous in solving sys- 
tematically certain problems in geometrical optics and appear to be of 
general mathematical interest. The method is not restricted to power 
series approximations, but is applicable to any system of approximating 
functions, P,, which have the quality that P,P, can be approximately 
represented as a linear combination of polynomials P, with r smaller than 
a constant ». However, in this note we shall restrict ourselves to poly- 


nomials. 
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In approximating functions by polynomials of one, two, or more 
variables, it is customary to neglect terms of degree greater than a given 
number n, so that the product of two polynomials of degree n again be- 
comes a polynomial of degree n. In this case, the polynomials of degree 
n form an algebra with the unity as a module and the quantities x’y*z 
as independent basis elements. The system is, in the language of algebra, 
a commutative, nilpotent system with unity added, since each of the other 
basis elements is a root of zero. 

Power Series in One Variable.—Let us consider the function 


f(x) = yo + Ax + Qy;x? + Pe ee x”. (1) 


Taking the case n = 5, for example, and considering 1, x, x”, ... x° as the 
basis elements, we obtain ihe multiplication table (assuming x” = 0 for 
m > 5) 




















which shows that the corresponding algebra is commutative. The theory 
of hypercomplex numbers' shows that the matrix: 


a0 
a) 
4 
ann 


#1 4111 





inn bl BOS! 





gives a regular representation of the function space, i.e., the product of 


i) 5 5 
two functions f(x) = Yaw, g(x) = bx, is a function h(x) = Vex, 
0 0 0 
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such that the matrix (3) formed from the coefficients is the product of 
the matrices (3) formed from the a; and 6,, respectively. 

The reader may notice that here, and in the following examples, only 
the first column of a matrix of the form (3) has to be computed; the others 
are then already determined. 

Power Series in Two Variables. Let us take as the next example the 
polynomial with two variables 
f = dy + AX + ayx? + ainx® + Aux? + V(d2 + AyyX + Ay2Xx* + 

Qy12X°) + V?(d22 + Ay2X + Ayi2X*) + V*(do22 + AiewXx) + Vo. 


The multiplication table becomes 























neil 


awe 


in which the blank spaces in the matrices stand for zeros, and, according 


to the theory of hypercomplex numbers, the regular representation is given 
by the group matrix 


} + 
“yin | 411 
*112 | 412 
4122 | #22 
4222 

} + 
®yynn) 4an 


} #1012) 4112 
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Three Variables.—Finally we shall investigate three variables to show 
the typical form for the parastrophic matrix as well as the group matrix. 


f(x, ¥, 2) = do + (aux + dey + ag2) + (Aux? + Gyxy + Ay2x2 + day? + 
o3V2 + 13327) + (dyyx* + Ay12X"y + dy3x*z + Ayx2X'y? + Aysxy2 + 
y33X2" + Aomy*® + 23,9" + A233V2" + 3332"). (7) 


We find as multiplication table 
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with the group matrix 














We have proved, for any polynomial function, the existence of a corre- 
sponding triangular matrix (containing in its first column the coefficients 
of the polynomials) such that the sum or product of two functions is ob- 
tained from the sum or product of the two corresponding matrices. This 
makes it possible to transform any equation between polynomials into 
equations between triangular matrices and thus to solve by an iteration 
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procedure any problem where the polynomial approximation is permissible 
in the neighborhood of the origin. 

An important application of the method is in the approximate solution 
of partial differential equations. 
of a os os , for the function (7) of three variables, keeping in 
Ox Oy Ox’? Oxdy 


We write down the group matrices 


mind that each differentiation decreases the approximation by one degree. 


re) O 
We find the matrices, A, for” and J as 
- 


et 








24122 
4123 
3ag22 


24223 











4233 


L 





nd going to the fourth power, we find for —, and 


24122 

aris 9129 
; } ' 

Nagy 4101 | 3ayyp2/2aq42 
| 





$4y122)24122 24112 

Sayyig/24i13 24912 S4yn 24432919123 
1 

2ay199] 


341222 





241223) 








"1233 | 23 ay | 


J 
(11) 

Let us give some examples which 
ow how easily partial differential equations are solved by this method. 





; ee of 
To find a function f(x, y) so that p = ~ and q = 


of ; 
obey the equation 
Oy 


Pr+¢? = (12) 
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we calculate with the matrices (10) and obtain the successive equations 
a,? + a? = 1 
40,;0;; + 2a, = O 
2ajdy. + 4Acd2n = O 
6a,Ai, + 4d? + 2Zdedy2 + Ay? = O 
4ajdyyo + 4a 42 + 4dody22 + 4dydy = O 
2a;A120 + dy” + 6222 + 4d? = 0. 


The reader will notice that, if we are given a), @), dy, our coefficients, with 
the exception of a2, can be computed linearly. This is a common occur- 
rence in this method. Having chosen for a, one of the roots of the first 
equation, we find 


Qoljy, = —2a\ay 

027d = 4,0 
A539 = — (3a \d27A44, + 204,") 
Ay4Ay29 = 0, (3a4027a111 + 4a)”) 


2°99 = — ,7(4097044) + 2d"). 


> 


As another example, we choose a differential equation which plays a 
certain role in optics and is hard to attack by classical methods 


(p? + @’)y = g. (15) 


The calculation is simplified by concluding that qg vanishes for y = 0. 
This means that the coefficients 


Q, = dy = Ay = due = Guy = 0. (16) 
We find from the matrix calculation the recursion formulae 


2d = a,’ 244122 = AYP ATUATT | + fay," 
2d) = 44,0) Sdi22 = O (17) 
3d22 z= () 409209 = 24422 + 4d" 


which determine all coefficients, if a), @, Gin, Gin are given. Up to the 
fourth order we find thus 


ay = 0 
Aye = 2a,ay 
Qu = 3d,dy1 + 2a4)° 
('/q)ay* + ay?ay,. 
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The method described here gives the development of the solution of the 
partial differential equation for the neighborhood of x = y = 0. It can 
break down only ‘f the zero point is a singular point. 

Of course there is, in general, no difficulty in obtaining the approximate 
solution for a series around an arbitrary point. In boundary problems, 
there remains the question of fitting the undetermined coefficients (q,, 
(ay), ...) to the boundary values. 

The method having been established, the reader will find no difficulty 
in applying it to the more common way of writing a polynomial, 


f( x,y = by +- bx -{- bey + (; 9) (b,x? + 2dyoxy + boy?) + (* 6) (Oy x3 + 
BO ox*y + Bdi20xy? + Doooy*) . ae (19) 


or to the problem of developing functions of several variables, instead of 
into a power series, into a series of other functions, having the quality that 


rove 


A,:A, = > aj), Ay, (20) 


where every product is a linear combination of the functions with index 
no greater than that of the highest of the original functions. 


* Communication No. 1589 from the Kodak Research Laboratories. 

' Cf. the matrices S; appearing in Frobenius, G., Theorie der hyperkomplexen Grossen, 
Berl. Sitzungsber., 1903, and Herzberger, M., Ueber Systems hyperkomplexer Grossen, 
[naugural Dissertation, Berlin, 1923 


ERGODIC THEORY OF MARKOV CHAINS ADMITTING AN 
INFINITE INVARIANT MEASURE* 
By T. E. Harrist AND HERBERT ROBBINS 
COLUMBIA UNIVERSITY AND THE INSTITUTE FOR ADVANCED STUDY 
Communicated by H. Whitney, June 3, 1953 


1. Introduction. —We consider real valued Markov chains 
Vo, X1, with a stationary transition probability function 


H(u, B) = Prob (x,+1¢€ Bix, = u). 
Under certain conditions there exists a probability measure m defined on 
the Borel sets B and satisfying 
r(B) = fr(du)H(u, B). (1) 
There is then a stationary Markov chain ... x1, Xo, %1, ... with the transi- 


tion function /7 such that Prob (x, « B) = 2(B),n = 0 + 1,....  Birk- 
hoff's ergodic theorem can then be applied to this chain (cf. ref. 2). 
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It may be that no probability measure m exists satisfying (1). (Such is 

the case, for example, if for every bounded Borel set B, lim //,(u, B) = 0, 
k—« 

where //, denotes the k-step transition function, //; = //.) There may, 
however, be an infinite measure + which satisfies (1). The purpose of 
this note is to show how the ergodic theorem of Hopf* can then be applied. 
A similar treatment will be given elsewhere for certain stochastic processes 
involving a continuous time parameter (cf. ref. 5). The present note 1s 
closely related to recent work of Robbins’ and of Kallianpur and Robbins.‘ 

2. Definition of the Measure m.—We shall assume henceforth that (1) 


‘ 


has an ‘“‘admtissible’’ solution x; 1.e., 7 iS a measure on the real Borel sets 
which satisfies (1), does not vanish identically, and is finite for bounded 
Borel sets. (Current unpublished work of C. Derman shows that this is so 
in a variety of circumstances.) 

Let Fy be the class of real Borel sets, 2 the space of sequences x of real 
numbers, x = (..., X-1, Xo, 1, ...), and F the Borel extension of the 
cylinder sets in & If A e F is determined by the coordinates x, X¢+1, 
..., then Q(A|x, = u) will denote the probability of A relative to the 
Markov chain starting with x, = u, as specified by /7. 

A measure m may be set up on F as follows. If A is a cylinder set 
determined by x,, ..., x,, let 


m(A) = J r(du)Q(A|x; = u), (2) 


where (1) ensures that any value of 7 < k may be used indifferently in (2). 
It can be verified that m may be extended to F, there being no essential 
change in the proof of the extension theorem of Kolmogorov. The 
extended measure m satisfies (2) whenever A ¢ F depends on xy, Xx41, 
and j < k. Moreover, m(A) = 0 if and only if Q(A|x, = u) = 0 

for a.e. (7) u; thus m(A) = 0 has a meaning in terms of ordinary Markov 
chain probabilities. From (1) it follows that if A « F and if 7 is the shift 
transformation, (7'x); = Xi, then m(A) = m(TA). 

3. The Dissipative Part of Q.—For any B ¢ Fy, let Ry be the event that 
x, € B infinitely often for m = 1, 2, (eo ae y, 

AssuMPTION |. Jf B ¢€ Fy then Q(Rg|\ xo = u) = 1 for a.e. (x) win B. 

THEOREM |. Jf Assumption 1 holds then the dissipative part of Q has m 
measure 0 (see ref. 3 for terminology). 

Sketch of Proof: Let A be a cylinder set defined by, say, (xo, x1, 
x,) € Beas, a (R + 1)-dimensional Borel set, and let 


U = tus < Q(A\x aan 
P Cr geagery eer ak 


Equation (2) ensures that A differs from }° A+(xp ¢ U’,) by a set of 
n=) 
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m-measure 0, and hence for a.e. (m)x in A, 7"x € A i.o., so that 


m(A+* >> 7-"A) = m(A). Next, let W be a ‘‘wandering’’ set with m(W) 
l 


< © (it suffices to consider this case) and let A be a cylinder set such 
that m(W — A) + m(A — W) <e. Then (all sums 7 = 1 to ~-) 


0 = m(WSST"W) > m(AsOT"W) — « 
= Sim|(T-"A)W] — « > m(W:>3T-"A) -— «> 
m(A*)“T-"A) — 2« = m(A) — 2, 


and hence 
m(W) < m(A) + € < 3e, m(W) = 0. 


4. Metric Transitivity of 7.--The following assumption is stronger 
than Assumption 1. 

AssumpTION 2. Jf B € Fy and 4(B) > O then Q(Rg\xo = u) = 1 for 
a.e. (mr) u in Q. 

THEOREM 2. Jf Assumption 2? holds then the shift transformation T is 
metrically transitive. 

Sketch of Proof: First, let A « F be such that for every k = 1, 2, ..., 
A depends only on xy, X¢+1, ... (call such a set a “‘terminal’’ set), and 
suppose A = 7A. Assumption 2 and the invariance of A can be shown 
to imply that Q(A|x» = wu) 1s constant for a.e. (w)u, and the 0-1 law 
shows that the constant is 0 or 1. It follows from (2) that m(A) = 0 or 
m(Q2 — A) = 0. 

Now let g(x) be a fixed positive m-summable function which depends on 
a finite number of coordinates of x. Hopf’s theorem implies that if f(x) 
is m-summable then 

> f(T’x) 


lim = u(f; x) 


n 
n> © ry. 
> g(7%x) 
0 
exists for a.e. (m)x. If f(x) depends on only a finite number of coordinates 


it follows from the result above on ‘‘terminal’’ invariant sets that 


Er S f(x)m(dx) 
S g(x)m (dx) 


To prove metric transitivity it is sufficient to show that (3) holds for 


u(f; x) = const. (3) 


every m-summable f. 
Now let f(x) be any m-summable function and let f,(x) be a sequence of 
functions of a finite number of coordinates such that 


lim J’ f(x) — f(x)| m(dx) = 0. 
ko 
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Then 


S fi(x)m(dx) 


(f> 2) = nlf — fox) : . 
uJ wf — Si JS g(x)m(dx) 


whence 


“f(x)m(dx) sae 
iu(f; x) Jf = lim |u(f — fis x) 


7 S g(x)m(dx) k 


Since for any m-summable h(x), S ulh; x)g(x)m(dx) = Sh(x)m(dx), 
it follows that 


lim S \ulf — fi; x)| g(x)m/(dx) < lim S wl f — fails x)g(x)m(dx) = 
ko « k— = 


lim J) if — fx| m(dx) = 0, 
k—>0 
and since g(x) > 0 the desired result follows. 

5. Applications.-In what follows h(u) and k(u) will be #-summable 
functions of a real variable u with J °k(u)x(du) # 0. Theorems | and 2 
in conjunction with Hopf’s theorem then imply 

THEOREM 3.° Under Assumption 2, for a.e. (m)Xxo, 


f,... A(xo) + + h(x») S.. h(u) (du) 
sibel lim = ; 


n—o R(X) + _ + R(x,) iF: ” k(u)a(du)S 


Coro.tiary. Under Assumption 2, if m is another admissible solution 
of (1) such that x(B) = O implies 1,(B) = O, then for some constant c > 0, 
Tw, = CT. 

Next, consider a sequence of independent random variables y, ys, 
with a common distribution function F(y). 

ASSUMPTION 3. The sums y, + . + y, are ‘“‘interval-recurrent’’; 1.é., 
for any interval I, Prob (y, + ... + y, €1 1.0.) = 1 (see ref. 1). 

We state without proof 

LemMMA |. Under Assumption 3, if B € Fy has positive Lebesgue measure 
then for a.e.a, Probla + y+... + ¥,€ Bio.) = 1. This ts true for every 
a if some convolution F"*(y) has an absolutely continuous component. 

Set H(u, B) = SydF(y — u); then x = Lebesgue measure satisfies 
(1), and from Theorem 3 and Lemma | we obtain 

THEOREM 4. Under Assumption 3, for a.e. (Lebesgue measure) a, 


d h(a + + FI) Of" hiw)de 


Prob ( lim ” =! my | 
‘ii p kia + vm + J . R(ujdu 
y=1 
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If some convolution F”’*(y) has an absolutely continuous component then 
(4) holds for every a. 

From Theorem 4 we can deduce the following 

COROLLARY. Under Assumption 3 (even if no F"*(y) has an absolutely 
continuous component) for every a the limit relation in (4) holds with proba- 
bility | simultaneously for every two functions h(u), k(u) which are Riemann- 
integrable in a finite interval and vanish elsewhere, the integral of k(u) assumed 
to be #0. 

Thus the prob | density of the sequence of sums y; + ... + ¥, implies 
the prob | eqguidistribution of these sums. This establishes a conjecture 
made in ref. 4. 

If the y,; have only integer values then the function 7(B) = number of 
integer points in B, is an admissible solution of (1). Theorem 3 then 
implies that under Assumption 2 (which holds whenever every integer 
visa ‘recurrent’ value of the sums y, + ... + y,), for every two functions 
of integers h(v), k(v) such that 


> |h(v)| < ©, Yl k(v)| < ©, Yk) #0, 


> AC. + +y) D> hv) 


Prob «lim - =- 
= 2, k(n 4+ y,) > kr) 
j= a 

* Work sponsored by the Office of Scientific Research of the Air Force. 

t On leave from the RAND Corporation. 

' Chung, K. L, and Fuchs, W. H. J., “On the Distribution of Values of Sums of 
Random Variables,’ Mem. Am. Math. Soc., 6, (1951). 

2 Doob, J. L., Stochastic Processes, Wiley & Son, New York, 1953, Chap. 10. 

5’ Hopf, E., Ergodentheorie, reprinted by Chelsea, New York, 1948, pp. 46-53. 

‘ Kallianpur, G., and Robbins, H., ‘On the Sequence of Partial Sums of Independent 
Random Variables,”’ to appear. 

® Kallianpur, G., and Robbins, H., ‘Ergodic Property of the Brownian Motion 
Process’; to appear in Proc. NATL. ACAD. Sct. 

® Kolmogorov, A., Grundbegriffe der Warscheinlichkeitsrechnung, reprinted by Chelsea, 
New York, 1946, pp. 24-30. 

7 Robbins, H., “On the Equidistribution of Sums of Independent Random Variables”’; 
to appear in Proc. Am. Math. Soc. 

* Chung has obtained this result independently and by a different method for the case 
in which the Y; have only a denumerable set of values. His results will appear in 
Trans. Am. Math. Soc. 
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ON COHOMOLOG ¥Y GROUPS OF COMPACT ANALYTIC VARIETIES 
WITH COEFFICIENTS IN SOME ANALYTIC FAISCEAUX* 


By K. KopDAIRA 
DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY 
Communicated by S. Lefschetz, June 9, 1953 


1. Introduction-Let V be a compact complex analytic variety of 
complex dimension » and let F be an analytic fibre bundle over V’ whose 
fibre is a complex line, i.e., the complex number field C and whose structure 
group is the multiplicative group C* of complex numbers acting on C. 
Moreover let {2’(F) be the faisceau' over |’ of germs of holomorphic p- 
forms with coefficients in F. The faisceau (?(/) introduced recently by 
D. C. Spencer and, independently by J.-P. Serre turned out to be of 
importance to applications of faisceaux to the theory of compact analytic 
varieties. However, for those applications, we need a basic theorem to 
the effect that the cohomology groups //*(V; Q°(F)) of V with coefficients 
in 2?(F) have finite dimension. The purpose of the present short note is 
to give an outline of a proof of this basic theorem. 

2. Harmonic Forms with Coefficients in Complex Line Bundles.—Take 
a sufficiently fine finite covering % = {U,{ of V. Then, letting be the 
projection of F onto V, r~!(U,) is a direct product: m-'(U;) = CX U, 
and ¢; X eC X U;j is identical with ¢, K ge C X U, if and only if ¢; = 
fix(2&, where f(z) is a non-vanishing holomorphic function defined in 
U, n U,. Under these circumstances we say that the bundle F 1s defined 
by the system of functions }\ fix}. 

Setting ay = |fyx|%, we infer that the system of functions {a} defines 
a principal bundle over V of class C“ whose structure group is the multi- 
plicative group of positive real numbers. Obviously this bundle is topo- 
logically trivial and therefore we can find a system {a,} of positive real 
functions a, of class C” defined, respectively, in UL’, such that 


lf? = @,/dk, inU, n Ux. 


By a form ¢ of type (p, q) of class C“ with coefficients in F we shall mean 
a system }y,} of exterior differential forms yg, of type (p, q) of class C” 
defined, respectively, in U, such that 


1 = fix’ Ge inU, n U,. 


We denote by I’ “(F) the linear space consisting of all forms ¢ of type 
(p, q) of class C” with coefficients in F. We introduce an arbitrary but 
fixed positive definite Hermitian metric ds? = 25° g,3 (dz“d2’) on V and 
denote by *¢, the adjoint form of g, with respect to this metric ds’. 
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Now, for any ¢ = fy)! Pak "(F), we define 0g = {(d¢),! and dg = | (d¢),! 
by 


(Og); = O¢,, 


l 
(dg); = ~ a2 -¢,) 
aj 


respectively, where O or O denotes the exterior derivation with respect to 
the variables 2“ or 2“. As one readily infers, Og and d¢ have invariant 
meaning, 1.e., we have (Og), = fix*(Og), and (dy), = fxe(de),. Therefore, 
Og « 1% **'(F) and dg ¢€ 1 "~'(F) for any g ¢ I’ “F). Next, for any 
pair of forms ¢, ¥ « Tl “(F), we define the inner product (¢, ¥) by 


; a 
(y, ¥) = I ej A*¥;, 
Vv a; 


1 1 
— OAV, = oe AM 
a; ay 


where we note that 


in U’, nm Ly. The linear operator d is the formal adjoint of 0 in the sense 
that 


(y, db) = (d¢, W) (1) 


for arbitrary form g ¢ 1 “(F), pe 1 **'(F). In fact, this formula (1) 
follows from the identity 


| | ] 
a( eA -t,) =— 09, Ath; — ~ 9 A *(dp);. 
ay (lj ay 


Finally we introduce the complex Laplace-Beltrami operator? 
( ] = d0 + Od. 
Then we have, for arbitrary forms ¢, y « T “(F), 
(le. ¥) = (¢ LW) = (Og, OF) + (d¢, dy). (2) 
Moreover we infer readily that the principal part of the differential operator 


3 ys, a8 


|} is — 2g that 1s, 


oO” 
02%02”” 
0? 
02°02" 


}= -De" 


where denotes the terms containing at most the first partial derivatives. 
Thus | | is an elliptic partial differential operator of the second order. 
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The theory of harmonic forms—that is the theory of differential cpera- 
tors, d, 6 and A = éd + dé—can be easily generalized to the differential 
operators, 0, } and | | = d0 + Od. In particular, letting //” “(F) be the 
subspace of I? “(F) consisting of all solutions ¢ of the Laplace-Beltrami 
equation | |g = 0, we have the decomposition formula*® 


I (F) = Or -\(F) @ dr “+"(F) @ HF), (3) 


where @ denotes a direct sum of mutually orthogonal linear subspaces. 
Since the equation | }y = O is equivalent to 0g = d¢ = 0, as (2) shows, 
I]” “(F) ts the linear space consisting of all forms ¢ of type (p, q) of class 
C” with coefficients in V satisfying Og = d¢ = 0. We call such ¢ a harmon 
form of type (p, q) with coefficients in F. Now, since |_| is a partial differ- 


o 
02702" 
defined over the compact variety V, we infer that the number of linearly 
independent solutions ¢ of the equation [ |g = O is finite. Thus we obtain 


y aB 


ential operator of e//iptic type with the typical principal part — Lg 


dim []” “(F) < + o, (4) 


3. Cohomology Groups with Coefficients in Q?(F).—Let Q’(F) be the 
faisceau over V of germs of holomorphic p-forms with coefficients in F 
and let //*(V; Q?(F)) be the gth cohomology group of V with coefficients 
in Q2?(F), In the special case in which F is the trivial bundle C X V, P. 
Dolbeault* has proved that //(V; 2?(C X V)) is isomorphic to the 0- 
cohomology of forms of type (p, g) on V. This result of Dolbeault can 
be easily generalized to our general case, as was pointed out by J.-P. 
Serre. Namely we have the isomorphism 

HV; Q?(F)) = Z” “(F)/or? *-"(F), (5) 
where Z” “(F) is the subspace of I” “(F) consisting of all ¢ « I” “(F) 
satisfying 0g = 0. 

Now we infer from the decomposition formula (3) that 

ZF) = OF? *-"(F) @ H® “(F). 

This proves the isomorphism 
Z” (F)/OT? 9 F) = H" (PF), 
Combining this with (5), we obtain 
HW" V; Q(F)) = HP). 
Therefore, by (4), the cohomology group //‘(V; Q?(F)) has a finite dimen 
sion. Thus we have proved the following 


THEOREM. The qth cohomology group Il V; Q?(F)) of a compact analytic 
variety V with coefficients in the analytic faisceau Q?(F) is isomorphic to the 
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linear space I”: *(F) of all harmonic forms of type (p, q) on V with coefficients 
in F. In particular, the dimension of H*( V; Q?(F)) is always finite. 


* This note was prepared while the author was partially supported by a research proj- 
ect at Princeton University sponsored by the Office of Ordnance Research, U.S. Army. 

! For the notion of analytic faisceau, see Séminaire, H. Cartan, 1951-1952. 

2 In the case where the bundle F is trivial, the operator {_] has been studied in detail 
by Garabedian, P. R., and Spencer, D. C., ‘‘A Complex Tensor Calculus for Kahler 
Manifolds,”’ Technical Report 17, Stanford University (1951). 

8 Cf. Kodaira, K., ‘‘Harmonic Fields in Riemannian Manifolds,”’ Ann. Math., 50, 
587-4565 (1949), Sections 4 and 10. 

‘ Dolbeault, P., ‘‘Sur la cohomologie des variétés analytiques complexes,’ Compt 
rend., Paris, 236, 175-177 (1958). 


GROUPS OF COMPLEX LINE BUNDLES OVER COMPACT 
KAHLER VARIETIES* 
By K. KopaAIRA AND D. C. SPENCER 
DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY 
Communicated by S. Lefschetz, May 29, 1953 


1. Introduction.—This is a preliminary report of a generalization of 
divisor class groups on algebraic varieties to groups of complex line bundles 


over compact Kahler varieties. In the case of algebraic varieties, a divisor 
class defines a complex line bundle whose structure group is the multi- 
plicative group of complex numbers. Therefore, on Kahler varieties, 
we replace the notion of divisor classes by that of complex line bundles, 
and determine the structure of the group of such complex line bundles. 
In this way we are led to the definition of Picard varieties attached to 


Kahler varieties. 

2. Groups of Complex Line Bundles.-We denote by C the complex 
line (i.e., the field of all complex numbers) and by C* the multiplicative 
group of complex numbers acting on C. Let V be a compact Kahler 
variety of complex dimension m and let F be an analytic fibre bundle 
over V’ with the fibre C and the group C*. In what follows we call such 
a bundle F simply a complex line bundle. Letting U = {U,} be a suffi- 
ciently fine finite covering of V’ and p be the projection of F onto V, we 
have p-'(U,;) = CX Uj, and ¢; X ze C X U;, is identical with ¢, *K 2€C X 
l’, if and only if ¢; = fix(z)&, f(z) being a non-vanishing holomorphic 
function defined in U; n Uy. We say that the bundle F 1s defined by the 
system of functions f,. Clearly we have 


Sir Sis = |, fihiehy = |, (1) 


Two bundles F, F’ defined, respectively, by f;,, fj, are said to be analytically 
equivalent if there exist non-vanishing holomorphic functions f, defined, 
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respectively, in U; such that fy = fifi'fn. Then we identify F’ with F 
and consider F" the same as F. 

Now, for any pair of bundles F, G determined, respectively, by fix, 
gx, we define F + G to be the bundle determined by fxg %'. Then the 
set § = |F} of all complex line bundles can be regarded as an additive 
group. The zero element of the additive group § is the trivial bundle 
Cx ¥. 

Given a divisor D on V defined in each LU’, by a local meromorphic fune- 
tion RD) = R,(z; D), the system of functions f(D) = R,(D)/R(D) 
defines a complex line bundle over V which will be denoted by {PD}. It is 
obvious that |D} coincides with |D’} if and only if D is linearly equivalent 
to D’. Thus each {D} may be identified with the divisor class of D and 
therefore the divisor class group can be regarded as a subgroup of &. 

3. Characteristic Classes.-It follows from (1) that, for each triple 
(j, k, 1), 


log fix + log fir + log fi, = 2rtcjpni (2) 


is a constant in U; n U, n U, and that cy; is a rational integer. The 
system of these integers c,, defines a 2-cocycle on the nerve N of the 
covering Ul and therefore determines an element cy of the 2nd cohomology 
group //*(N; Z) with integer coefficients Z. Moreover, if WU’ < U in the 
sense that each member U; Ul’ is contained in some U; ¢€ Ul, then the 
canonical homomorphism //*(N; Z) — IHP?(N’; Z) maps cy onto Cy’; 


thus c = lim cy is a well-defined element of the 2nd cohomology group 
\ 


IT?(V; Z) of the variety V, where lim denotes the direct limit with respect 
N 
to the partial order <. This ¢ is called the characteristic class' of the 
bundle F. We denote the characteristic class ¢ of F by c(F). Clearly 
the mapping F — c(F) is a homomorphism of the group § into //*(V; Z). 
Now we shall show that, for given ce II?(V, Z), there exists a complex 
line bundle F with «(F) = c tf and only tf the harmonic part Hc of the co- 
homology class c is of type (1.1). Let A” ” be the faisceau over V of germs 
of forms of class C” of type (q, p), =” ’ the subfaisceau of A“ ’ consisting 
of germs of forms closed under 0. Again, let B’ be the faisceau over V 
of germs of r-forms of Class C” and #’ be the subfaisceau of germs of r- 
forms closed under d. Clearly 2” = =” ° is the faisceau of germs of holo- 
morphic g-forms and C = ° is the faisceau of complex constants. We 
have the diagram 


> BP 
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where the two horizontal lines are exact, where II = Ifo, » is the canonical 
projection of B’ onto A® ”, and i denotes the inclusion map. We obviously 
have commutativity, that is, 7IIo , = Io, ,t, Ollo, » = Io,»41¢. In par- 
ticular, for p = 0, we have 


1 d 
()— eer “a4 eennbinienioi B° —— gi —— 


i | II | II | 
| | | 
+ ¥ 4 
m) 
0 >? ——+ Ao 0-4 E014) 
From the corresponding exact cohomology sequences we obtain the 
diagram 
HV; Z) 


Y 


> H*(V; C) 


ie 
| 


M 
HI” *(V) > FLAC SP), 


where //’(\’) is the space of harmonic p-forms on V, H® ?(V) is the sub- 
space of harmonic form of type (0, p), A, # are isomorphisms onto, and 
ull = i*\. The isomorphism ) is the classical one defined originally by 
Hodge, and the isomorphism y» is due to Dolbeault.* Applying this to 
the case p = 2, we see that, for given ce H?(V; Z), (1*)2c(e H2(V; 2°)) 
is equal to zero if and only if IA~%*c is zero, while, since A~4*c = He, 
A> 't#c = Ip effec. On the other hand, we infer from (2) that there 
exists F with c(F) = c if and only if (1*)*c = 0. Thus we see that there 
exists F with c(F) = c if and only if IIo ./7c = 0. Since c, and hence Hc, 
is real, Ilo .//c = O implies Il, o//c = 0. Consequently there exists F 
with c(F) = c¢ if and only if the harmonic part Hc of ¢ is of type (1, 1), 
q. e. d. 

Thus the homomorphism F — c(F) maps & onto the subgroup H}, (V; Z) 
of H?(V; Z) consisting of all c such that Hc = 1h, Ic. Letting 2B be the 
kernel of this homomorphism, we obtain 


B/B = Hi, (V; Z) (3) 
4. Picard Varieties~—Consider a bundle F « 8 defined by f,. Then, 
choosing rational integers c;, suitably, we infer that the holomorphic func- 
tions hy = log fy + 2 mic, satisfy 
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Ny thee thy = 0, inU,;n Uzn U;. 


Thus the system {/,,} defines an element /y of J/'(.N; Q) for each covering 
U. Moreover, if Ul < WW’, the canonical homomorphism //'(N’; Q) — 
T1'(.N; Q) maps hy onto hy. Thus h = lim hy is a well-defined element 


of 17'(V; 2%). The diagram 


uu 
HT (V) ICES P&P) 


shows that @ = wo'h « HI” '(V) ¢ H'(V) is the conjugate of a simple 
differential a of the first kind on V and h = 7*\@&. In this manner, each 
bundle Fe ¥ is determined by an eiement &@ « //° '(V’) and, moreover, 
the correspondence & — F is a homomorphism from //° '(V’) onto B. It 
is clear that the bundle F is trivial if and only if h = (1*)*(2 mic), c « H'(V; 
Z), while the relation h = (1*)*(2 mic) is equivalent to a = wh) = 
2 willy Jc. Thus we see that the bundle F is trivial if and only if the 
corresponding & has the form & = 2 milo ,J/c, ce H'(V; Z). Let 2g be 
the first Betti number of V. Then //° '(V’) is a complex vector space of 
dimension g and 
Y= fwlw = 2wrello fc, celitV: 2) (4) 

is a discrete subgroup of //° '(\V’) generated by 2g vectors which are 
linearly independent with respect to real coefficients. Hence the factor 
group //® '(V)/ 3 is a complex torus of dimension g and the isomorphism 


B= H*V)/3 (5) 


holds. 

In order to determine the explicit form of the correspondence & — F, 
take a point 2(7) in each U’,e U. Then the element hy e //'(N; 2) 
approximating h = 71*\& is given by the system of constants 


*2(k) = 
hy = e a7) Oe 
Hence the bundle F corresponding to & is (equivalent te the bundle) 


defined hy the system of constants 


*2(k) 
fix = OCXD Sutjy &. 
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Let a = >> 1,4, where a, are a base for the simple differentials of the 


y= 
first kind on V. Then the above formula shows that the bundle F ¢« 3% 


depends “‘analytically”’ on the parameters ft), f2, ..., t,. Thus it is reason- 
able to introduce a complex analytic structure on $2 by identifying ‘8B with 
the complex torus //® '(V)/¥. Wecall B = H® '(V)/Q the Picard variety 
attached to the Kahler variety V. This is a generalization of the Picard 
varieties attached to algebraic varieties over the complex number field.‘ 


» 


Thus we have proved the following 

THrorem. The group ¥ of all complex line bundles over a Kahler variety 
V contains the Picard variety B = H% '(V)/¥ attached to V as a subgroup 
and the factor group &/% is isomorphic to the subgroup Hi, ,(V; Z) of the 
2nd cohomology group IT’?(V; Z) with integral coefficients consisting of all 
ce ll(V; Z) whose harmonic parts Hc are of type (1, 1). 

* The authors were partially supported by a research project at Princeton University 
sponsored by the Office of Ordnance Research, U. S. Army Ordnance, while this 
note was being prepared. 

' Chern, S., ‘Characteristic Classes of Hermitian Manifolds,’’ Ann. Math., 47, 85 
121 (1946) 

2 Séminaire H. Cartan, 1951-52 

* Dolbeault, P., “Sur la cohomologie des variétés analytiques complexes,’’ Compt 
rend. (Paris), 236, 175-177 (1953). 

* Chow, W. L., ‘On Picard Varieties,’"’ Am. J. Math., 74, 895-909 (1952); Igusa, J., 
“On the Picard Varieties Attached to Algebraic Varieties,” Jbid., 1-22; Wiel, A., ‘On 
Picard Varieties,’’ [bid., 865-894. 


DIVISOR CLASS GROUPS ON ALGEBRAIC VARIETIES 
By K. Koparra AND D. C. SPENCER 
DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY 
Communicated by S. Lefschetz, May 29, 1953 


In a previous paper! we have determined the structure of the group of 
complex line bundles over a compact Kahler variety. In this paper we 
show that for an algebraic variety this group is isomorphic to the divisor 
class group, and we obtain, in particular, a new proof of the Lefschetz- 
Hodge theorem concerning algebraic cycles. Moreover, we prove Igusa’s 
first and second duality theorems. We continue with the notation of the 
previous paper. 

Let V be a non-singular algebraic variety of dimension n imbedded in 
a projective space. As was shown in Section 2 of the previous paper, 
every divisor D on V determines the corresponding complex line bundle 
{D} over V which may be identified with the divisor class of D. In fact, 
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letting Ul = {U} be a sufficiently fine finite covering of V and R,(D) = 
R,(z; D) be local meromorphic functions defining D in L’;, the bundle 
|D} is determined by the system of non-vanishing holomorphic functions 
fx(D) = R\D)/R(D). Clearly |D + D’| = {D} + {D’} and therefore 
the divisor class group may be regarded as a subgroup of the group ¥§ of 
complex line bundles. 

Let Q(F) be the faisceau over |’ of germs of analytic sections of the 
bundle F. An element of Q(F) is therefore a mapping 2 —~ ¢)(z) X z ofa 
small domain U’, ¢ V into p-'(U’,;) = C X U, ¢ F, where £,(z) is a holo 
morphic function of z. In case F = {D}, the mapping ¢)(2) K 2 => 
6(2)/Rj(z; D) maps Q(}D}) isomorphically onto the faisceau 2(D) of 
germs of meromorphic functions which are multiples of —D. Therefore 
we have the isomorphisms 


HV, QU D})) = HYV; O(D)). 
In particular we have 
dim 19(V; Q({D})) = dim | D} + 1. (1) 


Let S be a general hyperplane section of V and let 2”"(.S) be the faisceau 
of germs of meromorphic n-forms which are multiples of — S. 

LEMMA: The cohomology group I14(V, 2"(S)) vanishes for q 2 1. 

This can be proved in the following manner. Let 2" be the faisceau 
over V of holomorphic n-forms and Q% ' the faisceau over S of holomorphic 
(n — 1)-forms. Then we have the exact sequence’ 

i 4 ' 


0—— 0” a(S) — > OE > 0, 


where P denotes the Poincaré residue operator. The corresponding exact 
cohomology sequence is 


6 
> HWS, QO" ') > 11 V7, 2") > 114 V,.2"(S)) 


Moreover, we have the commutative diagram 


6* 
He-(S, 25") > HV, 2") 
| 


4 
H 
al fe) >H™YV), 


where yu is the isomorphism due to Dolbeault and where // is the mapping 
g > Ile, gell"~' *(S), Hg being the harmonic part of g considered 
as a current of type (m,qg) on V. Considering the linear spaces //"~% "(V), 
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H" ©" S) which are dual to 77" “(V), 77" *~'(S) in a canonical manner, 
we infer that the restriction map 


een) + H*-© %S) 


in the adjoint of the map //. Now, by virtue of a theorem of Lefschetz, 
the map r is an isomorphism into or onto according as g = | or g 2 2. 
Consequently //, and hence 6*, is a mapping onto or an isomorphism onto 


according as gq = lorqg2 2. Therefore we infer from the exact cohomology 
sequence (2) that //4(V, 2"(S)) = Oforg2 1,q.e.d. 

Letting AK be the canonical divisor on V, we have the isomorphism 
(DD) = o2"(D — K). Combining this with the above lemma, we obtain 
the following result :* /f the complete linear system | D) is sufficiently ample 
in the sense that |S| = |D — K| is a system of hyperplane sections of V in 
one of its ambient spaces, then the cohomology groups H%(V, Q(D)) vanish 
forq2 1. 

THEOREM |. For every complex line bundle F over an algebraic variety 
V, there exists a divisor D on V such that {D\ = F. 

Proof: We prove the theorem by induction on the dimension n of the 
variety V. First we consider the general case in which nm 2 2. Let S 
be a general hyperplane section of V which is fixed and let S,, be a general 
hypersurface section of V of order m. We denote by FS the restriction 
of F to S and by Q(F-S) the faisceau on S of germs of analytic sections of 
FS. We have the exact sequence 


. 


U r 
0) ——> Q(F — {S}) ——> a(F) ——> AF-S) —— 0 (3) 
where 7 is the inclusion map defined by 
i(2z) X 2 —— Ris; S)&(s) X 3, 


and where yr is the restriction. 
Now, let F be a given bundle. Setting 


Fn = F + [Snj 
we have, by (3), the exact sequence 
Ul 
0 2(F 4,1) > *OQ(F,,.) —— 9 OFF 5« 5) 


and therefore the corresponding exact cohomology sequence 


:* r* 
0 > FHV, QF »-1)) — HV, Q(F,)) ——a HS, Q( Fy: S)) 
5* a* “al 


> ID'(V, QO Fn s)) >ID(V, QUE ,)) - > 
H(S, O(F,,:5S)) -————> ... (4) 
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By the inductive hypothesis, there exists a divisor D on S such that’ 
1D! = F-Son S. We have therefore 


F-S = F-§ + {S,}S = {D + S,,-S}, on S, 


where S,,-.S is the hypersurface section of S cut out by S,,. Obviously 
the complete linear system |D + S,,-S) is sufficiently ample if m is large. 
Hence, by the above result, there exists an integer my such that 


ID(S; Q(F,,S)) = H(S, 2(D + S,,-S)) = 0, for m = mp. 


Since the cohomology groups //*(V, Q(F,,)) have finite dimensions,‘ we 
therefore obtain from (4) 


dim 77°(S, Q(F,,:.S)) = dim H°(V, Q(F,,)) — dim H°(V, Q(F»—1)) + 
dim //'(V, Q(F,,-1)) — dim H7'(V, Q(F,,)), (m > my) 


Adding both sides of this formula on m, we get 


dim /1°(V, Q(F,,)) = >> dim /1°(S, Q(F,,-.S)) — do, (5) 
m = mo + | 
where dy = dim //"(V, Q(F,,)). Now, since F,,-S = |D + S,,-S}, we 
have, by (1), 


dim 17°(S, Q(F,,:S)) = dim |D + S,,-S) + 1—~+ +o for m > +o. 
We therefore infer from (5) that 
dim /7°(V, Q(F,,)) ~ + © form—~> +o, 


This shows that, for large m, the bundle F,, has a global analytic section 
¢j(s) X 2 different from 0 X z. Letting D,, be the divisor on V defined 
by ¢)(z) = Oin each U,, we infer readily that F,, = |D,,{. Hence, setting 
D = Dan — Sm, we obtain F = {D}. 

In the case n = 1, H°(S, Q(F,,:S)) = QUF,,:S) = Q(S) is a linear space 
of dimension s 2 1, s being the degree of the divisor S, and J/'(S, Q(F,,-.S)) 
vanishes for any m. Hence, applying the above argument to this case, 
we infer the existence of D with |D| = F. This completes our inductive 
proof. 

Letting //,,.(V;Z) be the (2n — 2)th homology group of VY with 
integral coefficients Z, we have the canonical isomorphism //7*(V; Z) = 
H»,»(V; Z) which maps each cohomology class ¢ ¢ //°(\’; Z) onto its dual 
*¢ € Fy, o(V; Z). Since H(*c) = Ic, the mapping c ~ *c¢ gives the 


isomorphism 


H?, ,) (V; Z) & HS 2 (V; Z), (6) 


where //$)) (V; Z) is the subgroup of H:,2(V; Z) consisting of all 


homology classes whose harmonic parts are of type (1, 1). Each divisor 
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D is an integral (2n — 2)-cycle on V. Moreover, it can be shown that 
the homology class of D is equal to the dual of the characteristic class of 
the bundle {D}: 
De*c({D}). (7) 
In fact, we choose a suitable covering U = {l/;| and a suitable corre- 
sponding simplicial decomposition of Vin which each vertex 2( 7) lies in 
U’, and we denote by s(jk/) the simplex with the vertices 2(j), 2(k), 2(/) 
where lL’, 9 LU, A U, is not empty. Let e; be the dual 2n-cell of 2(j) 
and let D be defined by Rj(z; D). Then we have 
2 wil(D, s(jkl)) = A;log Rj(z; D) + A, log R,(z; D) + A, log Riz; D) = 
log fix(D) + log f(D) + log f(D) = 2 TIC jp 1. 


where A, log R,(z; D) is the variation of log R,(z; D) around the boundary 
of e, 9 s(jkl) and where the branches of log f;,(D), ... are chosen in the 
obvious manner. This proves the formula (7) above. Formula (7) 
shows that every divisor D belongs to the group H$}_'} (V; Z). 


Now, let G be the additive group of all divisors D on V, G, the subgroup 
of G consisting of all D which are homologous to zero with respect to 
integral coefficients, and let G, be the subgroup of G, of all D which are 
linearly equivalent to zero. Then it follows from Theorem | that 

G/G. = 8 (8) 
In view of (7), D belongs to G, if and only if {D} «8. Hence we get 
THEOREM 2. (Jgusa's first duality theorem.)® We have 


Ga/Gi = B (9) 


where YS = H%'(V)/X ts the Picard variety attached to V. Combining (6), 

(8), and (9) with the theorem in the previous paper, we get §/G, & H3)_'} 

(V; Z), while (7) shows that this isomorphism coincides with the natural 
i 

inclusion map G/S, — > HS!) (V; Z). Hence we obtain the following 


TueoreM 3. The factor group G/G, coincides with H)},"} (V; Z): 


G/Ga = Hyy-'2 (V; Z) (10) 


This leads immediately to the following criterion for algebraic (2n — 2)- 
cycles due to Lefschetz and Hodge :® 

THEOREM 4, An integral (2n — 2)-cycle T on an algebraic vartety V 
of dimension n is homologous to a divisor if and only if the harmonic part 
HT of T ts of type (1, 1). 

Let G, be the subgroup of G consisting of all divisors which are homolo- 
gous to zero with respect to rational coefficients, and let 7,(1’) be the 
gth torsion group of V. Since 7», .(V) ¢ Hi}: . (V; Z), we infer from 


2n 
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(10) that S,/Ga = Ton-2(V), while 7>,2(V) is dual to 7\(V). Conse- 
quently we obtain the following 

THEOREM 5. (/gusa’s second duality theorem.)' The factor group 
G,/ Ga ts dual to the first torsion group T\(V) of V. 


' Kodaira, K., and Spencer, D. C., these PROCEEDINGS, 39, 868-872 (1953). 

? Kodaira, K., and Spencer, D. C., “On Arithmetic Genera of Algebraic Varieties,”’ 
Ibid., 39, 641-649 (1953). 

§ This result has been pointed out by J. P. Serre as a consequence of his duality 
theorem. 

* Kodaira, K., ‘On Cohomology Groups of Compact Analytic Varieties with 
Coefficients in Some Analytic Faisceaux,’’ these PROCEEDINGS, 39, 865-868 (1953). 

* Igusa, J., ‘On the Picard Varieties Attached to Algebraic Varieties,’’ Am. J. Math., 
74, 1-22 (1952). Weil, A., ‘On Picard Varieties,’ /bid., 865-894 

® Hodge, W. V. D., The Theory and Applications of Harmonic Integrals, Cambridge 
University Press, 1941; p. 214-216. 

’ Igusa, J., loc. cit 


THE GROWTH OF AMBLYSTOMA CHIMERAS 
By GILBERT CHURCH 
DEPARTMENT OF BIOLOGICAL SCIENCES, STANFORD UNIVERSITY 
Communicated by V. C. Twitty, June 1, 1953 


The production of chimeras by combining anterior and posterior halves 
of Amblystoma tigrinum and A. punctatum embryos at the tail bud stages 
was shown to be possible by Stone in 1934.' Although normally A. 
tigrinum far surpasses A. punctatum in growth rate and size, Stone reported 
that the anterior component of the chimeras regulated the growth of the 
combination, at least up to metamorphosis. According to his findings, the 
“anterior half was always species dominant in its effect.’’ That is, the 
anterior punctatum halves reduced, and the anterior tigrinum halves in- 
creased, the intrinsic growth rate of the posterior halves. Although not 
specifically stated in Stone’s report of his work, it seemed to be implied 
that this species dominance of the anterior half resulted in the formation of 
harmoniously proportioned chimeras. 

These results, as Twitty later pointed out in his review of size-controlling 
factors,? were in contrast with the results obtained in heteroplastic grafts of 
eyes and limbs between the two species. In these two instances the trans- 


planted organs maintained their normal growth rates and achieved normal 
size.** There are, however, several instances in which size effects are ex- 
erted between structures intimately related in their development, such as 
the development of harmoniously proportioned eyes when the optic vesicle 
and lens epidermis are interchanged between A. tigrinum and A. puncta- 
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tum.’ Since many organ systems were involved in the growth of Stone’s 
chimeras, Twitty suggested that the total regulation might represent the 
cumulative effect of several contributing factors. But he pointed out that 
before a final interpretation could be given Stone’s results, further experi- 
mentation was needed. 

Accordingly, the present writer has repeated and extended Stone’s ex- 
periments employing A. ftigrinum and A. punctatum embryos both from 
Tennessee and from the Chicago area. Because it seemed conceivable 
that the dominance of the anterior part might vary with the proportions in 
which the two components were united, the fusions were made at levels 
ranging from in front of the gills to behind the anus. Most of the animals 
and controls were fed maximally, but some cases were purposely not. 

The results differ significantly from those achieved by Stone. Except 
for slight regulation at the line of fusion, necessary for the connection of the 
two unequally proportioned parts, in no instance has the anterior compo- 
nent been observed to exert a significant effect on the size of the posterior 
part. (See figures.) Each part has consistently maintained up to metamor- 
phosis a size an hereditary pattern characteristic of its species irrespective 
of the level at which the union was made. Although some evidence has 
been accumulated to show that partial size adjustment in the posterior 
part does occur six to eight months after metamorphosis, the findings up to 
metamorphosis are in agreement with those obtained from heteroplastic 
grafts of eyes and limbs between these two species, and lend themselves to 
similar analytical study and interpretation. 

Results of further anatomical and histological studies will be forthcoming. 

' Stone, L. S., Proc. Soc. Exptl. Biol. Med., 31, 1084 (1934) 

2 Twitty, V.C., Growth Suppl., 108-120 (1940) 

8 Harrison, R. G., Proc. Nati. Acap. Sctr., 10, 69 (1924). 

‘ Twitty, V. C., and Schwind, J. L., J. Exptl. Zool., 59, 61 (1981). 

5’ Harrison, R. G., Arch. Entwmech., 120, 1 (1929). 


FIGURE 1 


(A) A. tigrinum head combined with A. punctatum body. (B) A. punctatum control 
(C) A. punctatum anterior part combined with A. tigrinum posterior part. 


FIGURE 2 


(A) A. punctatum control. (B) Lateral view of chimera shown in figure 1 (C). 
(C) A. tigrinum control. Brackets will assist the reader in comparing the size of the 
chimeric and control components 
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DEVELOPMENT OF THE FINE STRUCTURE OF THE MYELIN 
SHEATH IN SCIATIC NERVES OF CHICK EMBRYOS* 


By Betty BEN GEREN AND JOSEPHINE RASKIND 


CHILDREN’S CANCER RESEARCH FOUNDATION, THE DIVISION OF LABORATORIES AND 

RESEARCH OF THE CHILDREN’S MEDICAL CENTER, THE DEPARTMENT OF PATHOLOGY, 

HARVARD MEDICAL SCHOOL; AND THE BIOLOGY DEPARTMENT, MASSACHUSETTS 
INSTITUTE OF TECHNOLOGY 


Communicated by Francis O. Schmitt, May 25, 1953 


Polarization optical studies" * of fresh nerve fibers produced evidence of 
the radial orientation of lipid molecules and the tangential orientation of 
protein molecules in the myelin sheath. The low-angle x-ray diffraction 
analyses* * of fresh nerves have defined a repeat period in the radial 
direction of the myelin sheath of about 172 A in fresh amphibian nerve and 
about 185 A in fresh mammalian nerve. In 1950, Fern4ndez-Moran® 
described concentric lamellar structures in frozen sections of osmium 
tetroxide fixed nerve fibers and reported the unit spacings to be about 70 A 
in rat sciatic and about 60 A in frog sciatic nerves. Sjéstrand® recently 
measured the unit spacing in sections of methacrylate-imbedded, osmium 
tetroxide fixed mouse nerve fibers and found an average spacing of 119 A. 
Our own measurements’ of osmium-tetroxide fixed, methacrylate-imbedded 
adult frog nerve spacings average about 85-95 A. 

It has been of interest here to distinguish between two possible modes 
of formation of the fine structure of the myelin sheath in the embryo. 
One possibility would be the ‘‘crystallizing out’ of the characteristic fine 
structure in accumulated lipid and protein at some time related perhaps 
to cellular activity or axonal function. An alternative possibility would be 
the deposition of lipid and protein as oriented layers from the beginning 
of the myelination process. The evidence presented here gives strong sup- 
port to the latter view. 

Methods.—-Fertile chicken eggs were incu vated for 18 days at 39°C. 
Embryos were removed, the sciatic nerves exposed and flooded with fixa- 
tive at room temperature. After preliminary blackening, the fibers were 
dissected and allowed to fix for a total period of 1 to 1'/: hours. The fixa- 
tive used was 1% osmium tetroxide in Hank's balanced salt solution.’ 
After washing, alcohol dehydration and infiltration, specimens were im- 
bedded at 45°C. in n-butyl methacrylate’ catalyzed with benzoyl peroxide. 
Sections were cut on the Minot-type microtome" with glass knives."! 
The electron microscope is an RCA type EMU fitted with an externally 
movable objective aperture.'” 

Results —The concentric lamellar structure of the myelin sheath is 
seen in cross section of the embryonic nerve fibers. The unit spacing 
(the distance from the beginning of one dense band to the beginning of the 
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next) was measured in undisrupted areas of 22 fibers (1056 total spacings 
measured). Figure 1 shows a distribution curve of these measurements, 
each position representing the average measurement for 3 to 14 spacings. 
The mean is 97 A with an average deviation of 9 A. Figure 2 (a) illustrates 
a region in a cross-section of a sciatic nerve, and figures 2 (>) and 2 (c) are 
higher magnifications of indicated areas in the myelin sheath. The thick 

ness of the dense lines in the spacing varies according to focus, thickness of 
section, orientation of the fiber in the plane of section, and degree of astig- 
matism caused by the objective aperture. In our best in-focus pictures 


15 


Total No. Fibers ae 
Total No. Positions 148 
Total No. Spacings 1056 
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FIGURE 1 


of through-focus series the width of this dense line varies between 20 and 
50 A. No regularly arranged structures were observed in the less dense 
band of the unit repeat. 

Within a single sciatic nerve are fibers with varying numbers of layers in 
the myelin sheath. Preliminary study of earlier ages of development 
(14 and 15 days) reveals many fibers with 2 to 12 concentric lamellae 
whereas in the 18-day material described here most fibers have 24 or more 
lamellae, although a few are seen with lesser numbers of layers. It is thus 
apparent that the increase in myelination of nerve fibers in the embryo 
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proceeds by the addition, in an oriented manner, of lipids and proteins to 
the surface of the uxonal projection within the cytoplasm of the Schwann 
cell. The Schwann cells in the early (14-day) embryonic nerve trunk are 
extremely large compared to the size of the individual myelinating fibers 
which appear as large Schwann cell cytoplasmic inclusions. This is in 
marked contrast to the satellite type of appearance of the Schwann cell 
around the adult fiber. Detailed studies of the earlier (14- and 15-day) 
stages will be reported later. 

Discussion.—\t is apparent that the myelin sheath has well-ordered 
fine structure in the 18-day chick embryo. Our observations indicate that 
from the beginning, the deposition of lipid and protein about the axon oc- 
curs in well-ordered layers and the thickness of the myelin sheath is in- 
creased by an increase in the number of layers. 

The unit spacing measures 97 + 9 A. One might expect that a unit 
structure approximately one-half the large period shown by x-ray diffrac- 
tion® * would be based on repeating lipid double layers arranged in some 
alternating fashion with the predicted tangential protein lamellae, assuggested 
by Schmitt, Bear, and Palmer. However, several factors make detailed 
studies of the period and correlation with the x-ray data still premature. 
Swelling and shrinking may occur during fixation with undefined changes 
of water content. One observes obvious areas of distortion of the myelin 
sheath as in figure 2 (c). The swelling and gross distortion of tissues im- 
bedded in methacrylate is easily recognized and such specimens are dis- 
carded, but the evaluation of small changes in volume due to imbedding is 
difficult. Fernandez-Moran’s data on frozen sections of fixed material are 
of significance here, but more data are needed for detailed comparisons. 
Very slight distortion of the specimen from heating in the electron beam 
can further prevent accurate measurements. Measurements of fibers cut 
somewhat diagonally with respect to the fiber axis would give greater than 
actual values for the unit spacing in the radial direction. In our specimens, 


the myelin sheath and its spacings appear to be compressed in the direction 


of cutting. 

Our lack of knowledge regarding the nature of the reactions of osmium 
tetroxide with tissue components further limits detailed analysis of the 
tissue fine structure as revealed by the present methods. 

Correlation of the fine structure of developing nerve fibers with the time 
of appearance of conduction cannot be attempted yet, but further advances 
in techniques and in interpretation of results may make such correlation 
possible. Investigations of this type may also permit analyses of factors 
influencing myelin formation at levels of resolution that allow interpreta- 
tion of the lipid and protein molecular arrangement. 

Conclusions.—The forming myelin sheath of embryonic nerve fiber has 
a repeat spacing similar to that of the myelin sheath of mature fibers. 
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FIGURE 2 
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The forming myelin consists of concentric lamellar structures arranged as 
a Sheath around the axon as it courses in the cytoplasm of the Schwann 
cell. The thickness of the forming myelin sheath is dependent on the num- 
ber of concentric lamellae and the younger the fiber the smaller the number 
of unit layers. 
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ACOUSTIC ORIENTATION IN THE OIL BIRD, STEATORNIS 
By DonaLp R. GRIFFIN* 
DEPARTMENT OF ZOOLOGY, CORNELL UNIVERSITY 
Communicated by H. A. Bethe, June 15, 1958 


The discovery that bats use high frequency sounds to orient themselves 
during flight has naturally given rise to speculation and inquiry concerning 
the possible existence of similar types of orientation in nocturnal birds. 
The limited studies of owls that have been made to date indicate that these 
birds rely primarily on vision, although their dark-adapted thresholds may 
be slightly lower than that of the human eye.' There are a few species of 
birds that habitually roost in caves, and some build their nests well inside 
caverns of various types where the light is dim at best. For example, the 
“edible bird’s nest swiftlets’’ (genus Collocalia) of southeast Asia are re- 
ported to nest in completely dark parts of caves.” As far as I know the 
orientation of these swiftlets has not been studied in any detail, but I 
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hope that the observations described below may serve to stimulate such 
investigations in the future. 

In northern South America there is another bird which roosts during the 
day almost exculsively in caves, and flies out to feed only during the night. 
This bird, Steatornis caripensts, the oil bird or guacharo, was described by 
von Humboldt in 1800, and it attracted wide interest among zoologists 
during the 19th century because it has no close relatives and is consequently 
placed in a family by itself. The bird was given its name because of the 
high fat content of the nestlings which were gathered in large numbers to 
yield edible fat. Steatornis is somewhat similar in structure to the family 
Caprimulgidae (nighthawks, whip-poor-wills, goatsuckers, and the like), 
but has a hawk-like bill and several anatomical pecularities such as the most 
extreme form of the paired, bronchial syrinx to be found in birds.* 

Over the past several years I have attempted to collect all available in- 
formation about the roosting places of these birds in caves, in particular 
the degree of darkness prevailing where they roost and fly. Out of several 
eye-witness accounts only a few stated with any assurance that the birds 
penetrate beyond the twilight zone, and it is certain that many of the nests, 
probably a majority, are located where some light is available. A few 
reports, however, gave strong indication that at least some of the birds 
fly into regions of total darkness. One of the most convincing of these came 
from William H. Phelps, Jr., on the basis of his visits to the Cave of the Gua- 
charos at Caripe, State of Monagas, Venezuela—the same cavern where 
von Humboldt first discovered the birds in 1800. As a result of our cor- 
respondence Mr. Phelps arranged for E. de Bellard in 1951 to expose some 
film for 4 or 5 minutes in the innermost chamber of the cave that is occu- 
pied by the oil birds. This roll of film appeared after development to be 
virtually undarkened, although a second, unexposed strip was very slightly 
more transparent. 

A still more convincing report was later received from O. P. Pearson in a 
letter describing a visit to las Cuevas de Tulini, Province of Tolima, Colom- 
bia, July 27-28, 1950. The following pertinent passage may well be quoted: 
“The ceiling (of the cave) is in places 75 or more feet high and there 
are additional holes and cracks containing many bats. . At dusk and 
after, two oil birds flew for perhaps a half hour near the ceiling of the 
cave...sometimes being chased by large bats (probably Phyllostomus 
discolor).... These birds made a loud snapping sound as well as a rau- 
cous rattling cry.... I looked through the cave at some time in the mid- 
dle of the night.... The oil birds were flying in the cave before dawn 
even with my eyes dark adapted I could not see my hand in front of my 
face. The birds were making a snapping sound as they flew, like a wooden 
ratchet, 5 to 10 clacks per second, perhaps an echolocation mechanism.” 

On March 14 and 15, 1953, I was able to accompany Mr. Phelps to the 
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Cave of the Guacharos at Caripe and to study this problem directly. Some 
of the apparatus used to study the high frequency sounds of bats was 
taken to the cave, and a portable, high fidelity tape recorder was also used 
to record the sounds made by the oil birds while flying through the cave. 
The other members of the party, all of whom rendered assistance in these 
studies were: Mrs. W. H. Phelps, Jr., Dr. and Mrs. W. MeCrudy, Dr. 
(;. Zuloaga and his son, G. Zuloaga, Jr. The trip was an extension of the 
author's studies of acoustic orientation in tropical bats supported in the 
Canal Zone by the Off-ce of Naval Research. But all expenses of the trip 
from Panama to Caripe were generously contributed by Mr. Phelps, and 
without his support these studies would not have been possible. 

Our first step on arrival at the cave was to walk into the deepest chamber 
where the oil birds roost and nest (650 meters from the entrance) and to 
remain there without lights of any sort for 25 minutes. This period was 
ample to attain complete dark adaptation, especially since it followed 15 
or 20 minutes during which we walked into the cave with only the light of 
flashlights and gasoline lanterns. Even after this complete dark adapta- 
tion none of us could detect any signs of light. On the following day I 
repeated this observation both during the morning and afternoon, and in 
addition exposed several frames of Eastman Super XX film for varying 
periods of time by removing the lens mount of a camera and pointing the 
emulsion toward the entrance to this chamber of the cave. Frames ex- 
posed for 5 and 9 minutes (the longest exposures made) showed no detect- 
able darkening when developed. Adjacent frames had not been exposed 
at all, and even though they underwent the same development procedure 
as the test frames, no line of demarcation could be seen in the finished roll 
between exposed and control areas. This chamber is 75 to 100 ft. high, and 
the birds roost on ledges near the ceiling. But the passage leading into 
this room is lower, and at the position where we dark-adapted ourselves, 
and where the film was exposed, there was more opportunity for light to 
penetrate than could possibly have reached the upper space where the birds 
were flying. It thus seems clear that these birds habitually fly in total 
darkness where vision cannot play any role in their orientation. 

Our next step was to make tape recordings of the sounds emitted by the 
birds in this innermost occupied chamber while we remained as quiet as 
possible without lights on the floor of the cave 50 to 75 ft. below them. The 
equipment used for this purpose was a Western Electric 640AA condenser 
microphone, a cathode follower, a Ballantine model 220 decade amplifier, 
and an Amplifier Corporation of America model 610E tape recorder, all 
battery operated. The protective grid was used over the microphone, but 
no horn was employed. The frequency response of this system was in all 
probability uniform from 50 to 15,000 c. p. s. within +6 db. 

All eye-witness accounts of oil bird colonies in caves, beginning with von 
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Humboldt’s‘ have stressed the clamor made by the birds when disturbed 
by visitors to the cave. Even to one thus forewarned the noise was hair- 
raising, so that one could feel only sympathy fer Humboldt’s Indian guides 
who refused to advance beyond a certain point in the cave because they 
believed the oil birds were in fact ghosts of their ancestors and that the in- 
ner parts of the cavern could only be entered by the dead. The sounds 
that we heard and recorded were of several sorts. Some were obviously 
call notes from birds perched on the ledges, others were squawks and 
screeches made as the birds flew back and forth above us, and still others 
consisted of the rapidly repeated ‘‘clucks” or ‘‘clicks’’ which Pearson had 
evidently heard in Colombia. The latter were especially noticeable when 
a bird flew down closer to where we stood. 

Of special interest were the sounds made by the oil birds as they spon 
taneously flew out through the cave in the evening. To record them the 
microphone was held at the top of a large rock 50-60 ft. below the ceiling 
of the cave and within 100 yd. of the entrance. In addition to the apparatus 
used in the inner room, the sounds picked up by the condenser microphone 
were also studied by means of a DuMont model 304H cathode-ray oscillo 
graph, powered by storage battery and vibrator. The oscillograph screen 
was photographed with a Bell-Howell Model 70DA 16-mm. motion picture 
‘amera in order to obtain graphic records such as figure 1. For some rec- 
ords and tape recordings a battery-operated Spencer-Kennedy model 
301 variable electronic filter was used between the amplifier and oscillo- 
graph; but unfiltered records and tape recordings were also obtained un- 
der all conditions of interest. These unfiltered records showed no appre- 
ciable components below 2000 c. p. s. and consequently the filter was used 
at 2000 c. p. s. high pass to improve the signal-to-noise ratio in records 
such as figures | and 2. 

The most striking fact about the sounds we heard and recorded as the 
birds flew out from the cave was the almost complete absence of calls and 
screeches such as those that predominated in the roosting chambers during 
the day. Instead there were almost nothing but loud, sharp clicks, 
repeated rapidly and almost continuously as the birds flew past. Only 
at long intervals would one of the hundreds of birds emit a longer duration 
sound like those that had bombarded our ears in the inner chamber. 

A typical oscillograph record of these clicks is shown in figure |; this 
was photographed directly from the oscillograph screen at the cave with 
the 2000 c. p. s. high pass filter. This picture is typical of hundreds of 
frames of the 16-mm. film, and while the wave-form varies to some extent 
from one record to another, the duration, frequency, and interval between 
pulses are approximately constant. 

As will be explained below, such records as figure | show that most of 
the energy in the clicks lay well within the frequency range of the tape 
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recorder. Hence the tape recordings provide a continuous record of ade- 
quate fidelity for most purposes. Figure 2 is a cathode-ray oscillograph 
record made from one of these tapes using a camera which moved the film 
continuously and which had a shutter that remained open throughout the 
period of film exposure. Study of many similar records, and watching 
the oscillograph while the tapes were being played back, show that figure 2 
is typical of the sounds recorded on the tape when only one bird was within 
range of the microphone. 

The frequency of these clicks emitted by the oil birds lies in most cases 
between 6000 and 10,000 c. p.s. No accurate time standard was available 
for the studies at the cave, so that the time scale for the direct cathode- 
ray oscillograph records such as figure | is probably less accurate than those 
like figure 2 made from the tapes. The average of 33 frequency measure- 
ments in records similar to figure 2 is 7300 c. p. s., the extremes of this 
series being 6100 and 8750. In pulses of such short duration the very con- 
cept of frequency begins to lose some of its usual meaning, and many of 
the records suggest that even within a single pulse the waves may differ 
slightly in their period or wave-length. Harmonics can be noted in the 
direct oscillograph records such as figure 1, and their presence interferes 
with precise measurements of frequency. The true amplitude of the har- 
monics is somewhat greater than the photographs indicate, because the 
sensitivity of the 640AA microphone declines at approximately 12 db per 
octave at frequencies above 12 ke. But even when allowance is made for 
this factor it is still clear that the principal sound energy is contained in the 
fundamental waves that lie well within the range of human hearing. While 
we have no direct knowledge of the frequency range of hearing in Steator- 
nis, birds in general are insensitive to frequencies above 20,000 c. p. s.° 

It should be emphasized in this connection that the apparatus which 
yielded figure | is the same as that which I have previously used for studies 
of the high frequency sounds of bats.’ Hence any strong components 
above 20 ke. would certainly be visible in the records, provided they 
reached the microphone. To be sure, the distance of about 50 ft. from bird 
to microphone would have caused considerable attenuation of ultrasonic 
frequencies; for instance, at 50,000 c. p. s. the attenuation in the air would 
be approximately 30 db.° It would thus be possible to contend that during 
flight out through the cave the birds emitted sounds having strong com- 
ponents at ultrasonic frequencies, but that these components were so 
greatly attenuated by 50 ft. of air that they are not apparent in the records 
such as figure |. But oscillograph records of the clicks emitted by one of 
the freshly caught birds te be described below revealed no signs of any 
such high frequency components, at least none having amplitudes approach- 
ing those of the fundamental frequencies below 10,000 c. p. s. Yet the 
sounds emitted by the only bird studied at close range by means of the os- 
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FIGURE 1 


Cathode-ray oscillograph record of two clicks emitted by an oil bird 
as it flew out through the entrance of the cave. 


Each sweep occupied 
5 msec. (0.005 sec.). 


The camera shutter was open for 16 msec. 
The frequency response of the apparatus was nearly uniform from 
2000 to 10,000 ¢, p.s. and declined above that frequency at approxi 
mately 12 db per octave. 


ee 
FIGURE 2 


Oscillograph record made from tape recording. In this picture the 


film moved vertically while the sweep circuit moved the oscillograph 


trace from left to right, each sweep requiring 16.7 msec. The fre 
quency response of the apparatus used was approximately uniform 
from 2000 to 15,000 ¢.p. s. 
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cillograph were probably only a limited part of its vocal repertoire, and 
further studies would be necessary before one could completely exclude 
the possibility that high frequency components are important under some 
conditions. 

The clicks emitted by the oil birds do not occur at a regular rate, but 
rather come in short bursts such as the one illustrated in figure 2. The 
number of clicks in such bursts varied from as few as two up to six or more; 
and the longer bursts seem to intergrade with certain of the call notes, 
although their true relation to other sounds made by the birds cannot be 
determined without further study. Within the bursts, however, the in- 
terval between clicks did not vary widely. Out of 43 intervals measured 
the average is 2.6 milliseconds, the extremes of this series being 1.7 and 4.4 
msec. 

The amplitude or sound pressure of these sounds could not be measured 
since the apparatus was not accurately calibrated and because the distance 
from bird to microphone and the angle of incidence were beyond our con- 
trol. But the clicks were readily audible even when the birds were 200 
yd. away; and all of our records were obtained as they flew at least 50 ft. 
above the microphone. The duration of the individual clicks is approx- 
imately | msec. Occasionally there are very short pulses containing only 
two or three waves (such as the first pulse of figure 2); but the longest 
pulse discovered among hundreds examined to date had a duration of al- 
most exactly 1.5 msec. 

It thus seemed clear that the oil birds fly through the virtually complete 
darkness of caves, and that while doing so they emit pulses of sound very 
similar to those used by many species of bats, but having a lower fre- 
quency—well within the range of human hearing and that of most birds. 
This evidence alone is not conclusive, although it certainly supports the idea 
that these birds guide their flight through dark caves by the type of acous- 
tic orientation I have called echolocation.* We therefore carried out a 
few simple experiments of the same type as those which had established 
the auditory basis of obstacle avoidance in bats.* 

For this purpose Mr. Phelps caught four live oil birds by means of a tram- 
mel net. Observations in the cave indicated that any significant experi- 
ments required a dark room, for in the light of our flashlights the birds 
greatly reduced the number of clicks and sometimes stopped them altogether. 
It therefore seemed best to carry the birds to Caripito where a darkroom 
was improvised the following evening through the cooperation of Dr. G. 
Zuloaga and the Creole Oil Company. This room was about 10 or 12 ft. 
square and § ft. high; and the birds flew about inside it reasonably well 
despite the fact that they had been in captivity for 24 hours and had not 
eaten any appreciable amount of food since the night before their capture. 
Either in light or in darkness they avoided any direct collisions with the 
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walls, although they did brush their wing tips from time to time as they 
turned in the room, which was only about four times their wingspread in 
its dimensions. In the dark their movements could easily be followed by 
the noise of their wingbeats as well as by the clicks which sounded approx- 
imately like those emitted while flying out through the cave. None of the 
birds was able to avoid frequent collisions with an electric light cord sus- 
pended from the ceiling; but it is quite possible that an object of this 
size could have been detected and avoided by oil birds flying under more 
natural conditions. 

The following tests were carried out with each of the three strongest 
birds, one at a time. After being allowed to fly about the room two or 
three times (in the dark) to show that it was able to avoid the walls success- 
fully, each bird was captured and its ears plugged with absorbent cotton 
and Duco cement. The cotton was packed fairly tightly into the internal 
‘ar canal but only as far as the first bend. The cement was then applied 
in a double layer over the cotton and the immediately adjoining skin. 
The bird was held quietly for two or three minutes to allow partial drying 
of the cement and then released in the dark. After two to five such flights 
the plugs were removed and the bird’s flight again observed in the dark 
with the ear canals unencumbered. 

In every case the bird with its ears plugged flew directly against the first 
wall it came to, and noisy collisions were the immediate result of all at- 
tempts to fly. Yet all three birds recovered their previous skill once the 
plugs had been removed. The general weakness of the birds made it im- 
possible to carry out more than a few flight tests in each condition— intact, 
-ars plugged, and intact again after removal of the plugs—but the results 
were uniform and consistent. The strongest bird gave the most impressive 
performance; for it continued to fly vigorously even after the ears were 
stopped, and began to show signs of weakness only after several collisions 
with the walls. Yet it recovered its skill and vigor after removal of the 
plugs, and at the end of the observations seemed as agile as the other birds 
had been at the beginning of their tests. 

With the room light turned on the untreated birds flew skilfully, but 
emitted many fewer clicks than when flying in the dark. One bird was 
flown in the light with its ears plugged. It emitted very few clicks, and had 
no difficulty avoiding the walls of the room. Evidently when light is avail- 
able the oil birds orient themselves by vision, as one might expect from 
the fact that their eyes are roughly as large as those of other birds of similar 
size. When flying in the dark with ears plugged, however, all three birds 
clicked vigorously, and almost without interruption. 

While this evidence should ideally be supplemented by more thorough 
studies of orientation under a wider variety of conditions, it already seems 
clear that sound plays an essential role in guiding the flight of Steatornis 
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through the dark. Not only can this bird now be added to the list of an- 
imals known to orient themselves by echolocation, but the sound frequen- 
cies employed by Steatornis are far lower than those used by any of the 
bats yet studied. Yet the oil bird, like most bats, employs pulsed sounds 
having a duration so short that the whole wave train in a single click ex- 
tends over less than a meter as the sound travels through the air. Short 
pulses have the advantage that echoes return even from objects at close 
range during a silent period when they can more easily be heard. While 
Mohres'” (1953) has recently described a special type of echolocation in 
bats of the family Rhinolophidae (the horseshoe bats of the Old World), 
in which pulses lasting roughly one-tenth second are utilized, all other bats 
yet investigated use pulses lasting only a few milliseconds, 

The fact that the oil bird can orient itself successfully during flight 
through dark caves by echolocation based on sounds having a frequency of 
roughly 7000 c. p. s. reinforces the idea that similar types of orientation 
are quite possible in our own species. Blind men are known to make ex- 
tensive use of sound in their orientation, and recent experiments show that 
hearing is the most important sense used by the blind to detect objects at 
a distance.'' The fact that bats can guide their flight by echolocation 
based on ultrasonic sounds has perhaps seemed to have only an academic 
interest for blind people. Perhaps the success of the oil bird will prove 
more encouraging. I even venture to hope it may stimulate efforts to 
improve the methods of acoustic orientation now used by the blind. One 


important aspect of such attempts should surely be a study of various types 
of sound that might be emitted in order to produce useful echoes, and a 
promising first step would be the use of very short pulses, or clicks, like 
those of bats and oil birds. 

Summary. 1. The oil bird, Steatornis caripensts Humboldt roosts dur- 
ing the doytime in caves, and many of the birds choose nesting sites in total 


darkness where vision could not be used for their orientation. This was 
established for the innermost occupied chambers of the Cave of the Gua- 
cheros at Caripe, Venezuela, by (1) noting that even a’ter 25 minutes of 
dark adaptation no light could be seen as we stood under the birds’ nesting 
ledges, and (2) by exposure film for 9 minutes without causing any detect- 
able darkening of the emulsion. 

2. When flying about inside the cave the oil birds make a variety of 
sounds, but when flying out for the night the predominant type is a sharp 
click, repeated in short bursts while the bird is on the wing. When light 
is available the clicks are much less frequent or cease altogether. 

3. The individual click consists of a very few sound waves having a 
frequency of about 7000 c. p. s. (average of 33 measurements 7300, ex- 
tremes 6100 and 8750 ¢. p. s.). The duration of each click is about | 
msec. the shortest measured to date being 0.3 msec., and the longest 1.5 





eee ek atone hee aa 


eet is eg ec ats et eee 


Vor. 39, 1958 ZOOLOGY; D. R. GRIFFIN 893 


msec. The interval between clicks within the short bursts varied between 
1.7 and 4.4 msec. in a series of 43 measurements, the average being 2.6 msec. 
While no measurements of intensity were feasible, the clicks were loud 
enough to be easily audible at a distance of 200 yd. within the cave. Ex- 
cept for their lower frequency these pulses of sound are very similar to 
those used by most genera of insectivorous bats for their acoustic orienta- 
tion. 

4. The external ear canals of three captive oil birds were plugged with 
absorbent cotton covered with Duco cement, and all three birds became 
disoriented wher: flying in the dark with their ears plugged. Before, and 
immediately after this treatment, they flew about within a small dark 
room avoiding collisions with the walls. But when their ears were plugged 
they never turned away from a wall and collided with every object en- 


countered. 
5. It thus seems clear that Steatornis guides its flight through dark 
caves by a type of acoustic orientation similar to that used by bats, but 


the oil bird employs for this purpose short pulses of sound that lie well 
within the frequency range of human hearing. 


* Present address: Biological Laboratories, Harvard University 
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